Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 173,567 A→C 5.9% intergenic (+123/+35) fhuB → / ← hemL iron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913173,5670AC 5.9% 322.0 / 7.5 112intergenic (+123/+35)fhuB/hemLiron(3+)‑hydroxamate import ABC transporter permease/glutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base A (49/56);  new base C (4/3);  total (53/59)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 115 total aligned reads displayed.

TATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCT  >  NC_000913/173418‑173716
                                                                                                                                                     |                                                                                                                                                     
tATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATGCACGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATaa                                                                                                                                                       >  2:781849/1‑150 (MQ=255)
  attgttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGGACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCCGGCGCGTCTTATCGGGCCTACAAACGGCGCCGAATCGTAGGTCGGATAAgg                                                                                                                                                     <  2:949728/149‑1 (MQ=255)
   ttgttgAGAAAGCAGAGCCGTTAATATATCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGACCGATGCGACGCGGGCGCGTCTTATCCGGCCTACAAACGCCCCCCAATCGAAGGTCGGATAAGGc                                                                                                                                                    >  1:1247074/1‑150 (MQ=255)
     gttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATATAGCGTTTACGCGACAACTGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAATCTGCCCCGAATCGTAGGTCGGATAAGGCGt                                                                                                                                                  <  1:1639975/150‑1 (MQ=255)
     gttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTTTCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATACGGCGt                                                                                                                                                  >  2:1032707/1‑150 (MQ=255)
     gttgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGt                                                                                                                                                  <  1:1015268/150‑1 (MQ=255)
       tgAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAAGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGCAACGTAGGTAGGATAAGGCGtta                                                                                                                                                >  1:1721750/1‑149 (MQ=255)
          gaAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGGCGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcg                                                                                                                                             >  2:932764/1‑150 (MQ=255)
            aaGCCGAGCCTTTAATTTATCCATATACCCGCCCTGTACGCCTGGACGGATAAAGCCTTTACGCGACAACGGATGCCCGATGCGACGCTGCCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAcgcc                                                                                                                                           <  2:345376/150‑1 (MQ=38)
               cAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACACCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTCCAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCa                                                                                                                                        >  2:1187728/1‑150 (MQ=255)
                  agCCGTTAATTTTTCCATATAACCGCCCTGTACGCCTGGACGGATAAAGCGTGTCCGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  1:148069/150‑1 (MQ=38)
                  agCCGTTAATTTTTCCATAAAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGTCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                     <  1:1927925/150‑1 (MQ=255)
                   gCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                    <  2:185363/150‑1 (MQ=39)
                   gCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                    >  2:1235102/1‑150 (MQ=255)
                   gCCGTTAATTTATCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                    <  1:1943320/150‑1 (MQ=255)
                          aTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAcc                                                                                                                             <  1:1150410/150‑1 (MQ=255)
                          aTTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAcc                                                                                                                             >  2:1335497/1‑150 (MQ=255)
                             tttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCCCGCCGCATCCGACAACCCATg                                                                                                                          >  1:1481522/1‑150 (MQ=255)
                             tttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATg                                                                                                                          <  2:318286/150‑1 (MQ=255)
                             tttCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATg                                                                                                                          <  2:534188/150‑1 (MQ=255)
                                ccATATAACGGCCCTGTACGCATGCACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTATTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTg                                                                                                                       <  2:1769427/150‑1 (MQ=255)
                                 cATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTgg                                                                                                                      >  1:1723064/1‑150 (MQ=255)
                                 cATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAAACGTAGGTCGGATAAGGCGTTCACGCCGCATACGACAAACCAAGCTgg                                                                                                                      >  1:1709642/1‑150 (MQ=255)
                                   tataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACTCGACAACGGATGCCCGATGCGACGCGGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGaa                                                                                                                    <  2:1310268/150‑2 (MQ=255)
                                   tataACGGCCCTGTACGCCTGGACGGATAAAGCGATTACGCGACAACGGAAGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAt                                                                                                                    >  1:118303/1‑150 (MQ=255)
                                    ataACGGCCCTGTACGCCTGGACGTATAAAGCGATTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATc                                                                                                                   <  2:6836/150‑1 (MQ=255)
                                    ataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATc                                                                                                                   <  2:1636578/150‑1 (MQ=255)
                                    ataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATc                                                                                                                   <  2:423023/150‑1 (MQ=255)
                                    ataACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATc                                                                                                                   <  2:452691/150‑1 (MQ=255)
                                     taACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATca                                                                                                                  <  1:1766527/150‑1 (MQ=255)
                                     aaccGGCCCTGTACGCCTGGACGGATACAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATca                                                                                                                  <  2:1687032/147‑1 (MQ=255)
                                       aCGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATcaca                                                                                                                >  1:1178820/1‑150 (MQ=255)
                                        cGGCCCTGTACGCCTGGTCGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGCATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAAACATGCTGGATCACaa                                                                                                               >  1:355535/1‑150 (MQ=255)
                                               gTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCa                                                                                                        <  1:895629/150‑1 (MQ=255)
                                                 aCGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGAAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTCGCTCACAAATTAGCaaa                                                                                                      >  2:1812763/1‑150 (MQ=38)
                                                    ccTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAacac                                                                                                   <  1:1758101/150‑1 (MQ=255)
                                                     cTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACAcc                                                                                                  >  2:1511160/1‑150 (MQ=255)
                                                      tGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACAccc                                                                                                 >  1:815161/1‑150 (MQ=255)
                                                        gACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGAGTTCCCGCCGCATCCGACAACCCATGCTGGATCACAACTTCGCACACACCCgc                                                                                               >  1:1488031/1‑149 (MQ=255)
                                                         aCGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGAc                                                                                              <  1:1086463/150‑1 (MQ=255)
                                                          cGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTATGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACg                                                                                             >  1:1981696/1‑150 (MQ=255)
                                                            gATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCAATCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAACCACACGACGTg                                                                                           >  1:684333/1‑150 (MQ=255)
                                                             aTAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACTAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgc                                                                                          <  1:1174326/150‑1 (MQ=255)
                                                              tAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgca                                                                                         <  2:364837/150‑1 (MQ=255)
                                                              tAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTgca                                                                                         <  2:1291482/150‑1 (MQ=255)
                                                                   cGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATc                                                                                    >  2:655127/1‑150 (MQ=255)
                                                                   cGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGAACACAACTTCGCAAACACCCGACGTGCAGCATc                                                                                    <  1:171995/150‑1 (MQ=255)
                                                                    gTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGTGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCg                                                                                   <  2:1115376/150‑1 (MQ=255)
                                                                      ttACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGCCAACCCATGCTGGATCACAACTTCGCAAACACCCGACGAGCAGCATCGAt                                                                                 >  1:1323263/1‑150 (MQ=255)
                                                                       tACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATACGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATg                                                                                >  1:106098/1‑150 (MQ=255)
                                                                       tACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATg                                                                                <  1:43525/150‑1 (MQ=255)
                                                                        aCGCGACAACGGATGCCCGATGCGACGCTGGCGGGTCTTATCGGGCCTGCAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATgg                                                                               <  2:408969/150‑1 (MQ=255)
                                                                         cgcgACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCCCGCCGCATGCGCCAAACCATGCTGGATCACGACGTCGCAACCACGCGCCGTGCAGCATCGaggcg                                                                              >  1:344470/1‑146 (MQ=37)
                                                                          gcgACAACGGATGCCCGATGCGACGCTGGCGCGTGTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGgtg                                                                             <  1:1850145/150‑1 (MQ=255)
                                                                          gcgACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCCTCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGCCGTGCAGCATCGATGgtg                                                                             >  1:1329029/1‑150 (MQ=255)
                                                                           cgACAAAGGATGCCCGATGCGACGCTGGCGCGTCTTATCTGGCCTACCAACGGACCCGAATCGTAGGACGGATAAGGAGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCCACGTGCAGCATCGATGgtgt                                                                            >  1:146725/1‑150 (MQ=38)
                                                                               aaCGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTAtt                                                                        <  1:1119653/150‑1 (MQ=255)
                                                                                  ggATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGat                                                                     <  1:1189734/150‑1 (MQ=255)
                                                                                   gATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGata                                                                    <  2:1313694/150‑1 (MQ=255)
                                                                                       cccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCtt                                                                >  1:264282/1‑150 (MQ=255)
                                                                                       cccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCtt                                                                >  2:1873222/1‑150 (MQ=255)
                                                                                       cccGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGGGTTATTGATATCtt                                                                >  2:102580/1‑150 (MQ=255)
                                                                                         cGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTcc                                                              <  1:1973966/150‑1 (MQ=255)
                                                                                               gACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCtg                                                        >  2:189964/1‑150 (MQ=255)
                                                                                               gACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGGGCAGCATCGATGGTGTTATTGATATCTTCCctgctg                                                        >  1:1891007/1‑150 (MQ=255)
                                                                                                   cTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTgcg                                                    >  2:1471338/1‑150 (MQ=255)
                                                                                                      gcgcGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCa                                                 <  1:1074550/150‑1 (MQ=255)
                                                                                                       cgcgTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCAc                                                >  2:1201426/1‑150 (MQ=255)
                                                                                                         ggTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACgg                                              <  2:1950325/149‑1 (MQ=255)
                                                                                                           tCTTATCGGGCCTACAATCGGCCCCGAATCGTAGGTCTGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGAc                                            <  1:20480/150‑1 (MQ=255)
                                                                                                              tATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                         <  1:1667752/150‑1 (MQ=255)
                                                                                                              tATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                         >  1:802827/1‑150 (MQ=255)
                                                                                                              tATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                         >  1:910491/1‑150 (MQ=255)
                                                                                                              tATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGATTGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATa                                         <  2:1396569/150‑1 (MQ=255)
                                                                                                               aTCTGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCTCGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATaa                                        <  2:1839418/150‑1 (MQ=255)
                                                                                                                  gggCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTGCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGc                                     >  1:221408/1‑150 (MQ=255)
                                                                                                                     ccTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGAGCAGCAGCGATGGTGTTATAGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCg                                  <  1:441715/150‑1 (MQ=255)
                                                                                                                       tACAAACTGCCCCGAATCGTACGTCGGATAAGTCGTTCACGCCGCATCCGAAAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTTTGCGCCACTGACATAAAGCCCGCt                                <  2:13659/150‑1 (MQ=255)
                                                                                                                       tACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCt                                >  1:47747/1‑150 (MQ=255)
                                                                                                                         cAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTc                              <  2:1583055/150‑1 (MQ=255)
                                                                                                                          aaaCGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                             <  2:815161/150‑1 (MQ=255)
                                                                                                                          aaaCGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCa                             <  1:1634583/150‑1 (MQ=255)
                                                                                                                             cGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                          <  1:1108161/150‑1 (MQ=255)
                                                                                                                             cGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                          <  2:1267369/150‑1 (MQ=255)
                                                                                                                             cGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAAc                          <  1:1831609/150‑1 (MQ=255)
                                                                                                                               gCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATATAGCCCGCTTCAAACGc                        <  1:1715771/150‑1 (MQ=255)
                                                                                                                                ccccGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                       >  2:1116015/1‑150 (MQ=255)
                                                                                                                                ccccGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCt                       <  1:611315/150‑1 (MQ=255)
                                                                                                                                   cGAATCGTACGTAGTATAAGGCGATCACGCCGCATCAGACAAACCATGCTGGATCACAACTTCGCTATCACACGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACTGACATAAAGCCCGTTTCAAACGCTGAc                    <  1:90498/150‑1 (MQ=255)
                                                                                                                                     aaTCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCTTGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACgg                  <  2:136821/150‑1 (MQ=255)
                                                                                                                                      aTCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGt                 >  2:1527236/1‑150 (MQ=255)
                                                                                                                                       tCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAACGCCCGCTTCAAACGCTGACGGTg                >  2:1094485/1‑150 (MQ=255)
                                                                                                                                        cGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTTTTGATATCTTCAATGCTGTTCGCCACTGACATAAAGCCCGCTTCAAACGCTGCCGGTGc               >  2:104810/1‑150 (MQ=39)
                                                                                                                                           ttgTCGGATGCGGCGTGAACGCCTTATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg            >  1:1411209/3‑150 (MQ=255)
                                                                                                                                           tggTCGGATAAGGCGTTCACGCCGCATCCGAGAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTTTTATTGATATATTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAg            <  2:118303/149‑1 (MQ=255)
                                                                                                                                               cGGATGCGGCGTGAACGCCTTATCCGACAAACCATGCTGGATCACAACATCCCAAACCCCCGACGTGCAGCATCGATGGTGTTAGTGATATCTTCCATGCTGTGCGCCCCGGACATAAGGCCCGCTTCAAAAGCTGACGGTGCCCTGTaa        >  2:1256472/1‑150 (MQ=255)
                                                                                                                                                ggATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGAAATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTaaa       >  1:858733/1‑150 (MQ=255)
                                                                                                                                                 taTAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAac      <  1:396529/149‑1 (MQ=255)
                                                                                                                                                 gATAAGGCGTTCACGCCGCATCCGACAAACCAAGCTGGATCACAAATTCGCAAACACCCGACGAGCAGCATCGATGGTGTTAGTGATATCTTCCCTGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACTCTGCCGGTGCCAGGTacac      >  2:1210722/1‑150 (MQ=255)
                                                                                                                                                     aGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCt  >  1:186596/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAAACCATGCTGGATCACAACTTCGCAAACACCCGACGTGCAGCATCGATGGTGTTATTGATATCTTCCATGCTGTGCGCCACGGACATAAAGCCCGCTTCAAACGCTGACGGTGCCAGGTAAACACCT  >  NC_000913/173418‑173716

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