Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,498,215 A→G 7.0% H39R (CAC→CGC)  PP_2193 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,498,2150AG 7.0% 157.4 / 3.0 57H39R (CAC→CGC) PP_2193ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (26/27);  new base G (2/2);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

TGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGT  >  NC_002947/2498072‑2498359
                                                                                                                                               |                                                                                                                                                 
tGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTGTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAg                                                                                                                                             >  4:142878/1‑149 (MQ=255)
  aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                           >  1:124429/1‑149 (MQ=255)
       cAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAgggggccga                                                                                                                                     <  3:129193/149‑6 (MQ=255)
           ccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACg                                                                                                                                  >  8:15658/1‑149 (MQ=255)
            cATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACgg                                                                                                                                 <  6:60998/149‑1 (MQ=255)
            cATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACgg                                                                                                                                 <  7:15658/149‑1 (MQ=255)
              ttCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCt                                                                                                                               <  5:208303/149‑1 (MQ=255)
               tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGAc                                                                                                                                   >  5:8001/1‑144 (MQ=255)
               tctcATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGAc                                                                                                                                   <  6:8001/144‑1 (MQ=255)
                   aTTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCg                                                                                                                           >  5:257061/1‑148 (MQ=255)
                               tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGAt                                                                                                              >  6:66147/1‑149 (MQ=255)
                               tCAGGAAGCTCATTCGTGCCGCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGAt                                                                                                              <  1:94215/149‑1 (MQ=255)
                                 aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGc                                                                                                            <  8:270596/149‑1 (MQ=255)
                                  ggAAGCTCATTCGTGCCTGTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCa                                                                                                           <  3:142878/149‑1 (MQ=255)
                                    atgCTCATTCGTGCCGCTGATGCGCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTGACCATGCTTTTATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAgg                                                                                                         <  7:281379/147‑1 (MQ=255)
                                      gCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAACAAGATGCAGGcc                                                                                                       >  3:295689/1‑149 (MQ=255)
                                          aTTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                           <  7:184520/109‑1 (MQ=255)
                                          aTTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                           >  8:184520/1‑109 (MQ=255)
                                                   tctGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAg                                                                                           >  8:39026/1‑148 (MQ=255)
                                                     tGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTc                                                                                        >  2:147516/1‑149 (MQ=255)
                                                       atccccAACCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTc                                                                                        <  4:295689/142‑1 (MQ=255)
                                                        tGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAc                                                                                     <  2:124429/149‑1 (MQ=255)
                                                         gCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAca                                                                                    >  7:165867/1‑148 (MQ=255)
                                                           tCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTg                                                                                  >  6:259908/1‑149 (MQ=255)
                                                           tCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTg                                                                                  >  2:107748/1‑149 (MQ=255)
                                                                  cTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCAGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGcc                                                                            >  8:227435/1‑148 (MQ=255)
                                                                       aCCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCg                                                                      >  8:64308/1‑149 (MQ=255)
                                                                           ccgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                                  <  7:64308/149‑1 (MQ=255)
                                                                              tCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                                  >  2:81448/1‑146 (MQ=255)
                                                                              tCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTgg                                                                  <  1:81448/146‑1 (MQ=255)
                                                                              tCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAcc                                                               <  5:302372/149‑1 (MQ=255)
                                                                                gcgGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCa                                                              <  1:107748/148‑1 (MQ=255)
                                                                                gcgGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAg                                                             <  7:227435/149‑1 (MQ=255)
                                                                                  ggTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgc                                                           >  7:198122/1‑149 (MQ=255)
                                                                                   gTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgcg                                                          >  2:259340/1‑149 (MQ=255)
                                                                                         ttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGTCCAACCGCTGGACCAGGCGCTGAc                                                     <  2:316901/148‑1 (MQ=255)
                                                                                         ttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGa                                                      <  6:127313/147‑1 (MQ=255)
                                                                                               aTGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTc                                              <  1:147516/149‑1 (MQ=255)
                                                                                                  ccTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGc                                            <  8:165867/148‑1 (MQ=255)
                                                                                                      gATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAc                                        <  3:102183/148‑1 (MQ=255)
                                                                                                              gCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGCACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGCCCTGCCCGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggc                               >  5:243881/1‑149 (MQ=255)
                                                                                                                  tGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTg                            >  2:154223/1‑148 (MQ=255)
                                                                                                                  tGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGc                           <  7:39026/149‑1 (MQ=255)
                                                                                                                    ccGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCACCCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAc                         >  5:27398/1‑149 (MQ=255)
                                                                                                                           ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCg                                          >  2:239779/1‑125 (MQ=255)
                                                                                                                           ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCg                                          <  1:239779/125‑1 (MQ=255)
                                                                                                                           ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgctg                   >  4:99168/1‑148 (MQ=255)
                                                                                                                           ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGt                  <  3:93832/149‑1 (MQ=255)
                                                                                                                                  cTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcaccac           <  6:27398/149‑1 (MQ=255)
                                                                                                                                   tCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcacca            >  1:39519/1‑147 (MQ=255)
                                                                                                                                   tCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcacca            >  1:38564/1‑147 (MQ=255)
                                                                                                                                     cggccccTGCGTGCAGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgc                     >  2:189834/4‑137 (MQ=255)
                                                                                                                                     cggccccTGCGTGCAGGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgc                     <  1:189834/134‑1 (MQ=255)
                                                                                                                                      ggccccTGCGTGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGc       >  6:247237/3‑149 (MQ=255)
                                                                                                                                          ccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACct                       <  7:104846/129‑1 (MQ=255)
                                                                                                                                          ccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACct                       >  8:104846/1‑129 (MQ=255)
                                                                                                                                           cTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGt  <  2:39519/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGT  >  NC_002947/2498072‑2498359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: