Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,327 A→C 5.7% intergenic (‑54/‑265) ygcW ← / → yqcE putative deoxygluconate dehydrogenase/putative transport protein YqcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3270AC 5.7% 426.0 / 12.2 137intergenic (‑54/‑265)ygcW/yqcEputative deoxygluconate dehydrogenase/putative transport protein YqcE
Reads supporting (aligned to +/- strand):  ref base A (61/68);  new base C (5/2);  total (66/71)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 140 total aligned reads displayed.

TGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAACAC  >  NC_000913/2900179‑2900476
                                                                                                                                                     |                                                                                                                                                     
tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                      <  1:1386576/150‑1 (MQ=255)
 ggCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                                                                     >  2:1331069/1‑150 (MQ=255)
   ccTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAgg                                                                                                                                                   >  1:94845/1‑150 (MQ=255)
    taaaaCCGCAATTGTTACCGGTGGGAATAGCGGTTTTAACTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg                                                                                                                                                  >  1:525091/3‑150 (MQ=38)
     tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                                                                 >  1:1021414/1‑150 (MQ=255)
     tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                                                                 >  2:370208/1‑150 (MQ=255)
     tAAACCGCTATTCCAACCGGTAACAATTGCGGTTTTACATTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAggg                                                                                                                                                  >  2:1258136/1‑150 (MQ=255)
      aaaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                                <  1:1065839/150‑1 (MQ=255)
       aaCGGCTATTCCCACCGGTAACATTTGGGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGTACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAg                                                                                                                                               <  2:1069835/150‑1 (MQ=255)
        aCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGt                                                                                                                                              <  1:909259/150‑1 (MQ=255)
           gCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAcc                                                                                                                                          >  1:19283/1‑150 (MQ=255)
             tATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCt                                                                                                                                         >  1:63673/1‑150 (MQ=255)
                tCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtc                                                                                                                                      <  2:342263/150‑1 (MQ=255)
                tCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtc                                                                                                                                      <  1:370208/150‑1 (MQ=255)
                 cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtct                                                                                                                                     >  2:429664/1‑150 (MQ=255)
                  ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCTATTGACATAGTTAGCTACAATTTCCGTATGGAAAAAAAAAAGTTAGGTAACTTCCCTCTCGGGGGTTAAGTTAGCTTtctt                                                                                                                                    <  2:260173/150‑2 (MQ=255)
                     ccGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTGCCTTTCTCAtt                                                                                                                                 >  2:956557/1‑150 (MQ=255)
                     ccGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCAtt                                                                                                                                 >  2:237189/1‑150 (MQ=255)
                      cGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTa                                                                                                                                >  1:1615636/1‑150 (MQ=255)
                       ggTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAt                                                                                                                               <  1:1315518/150‑1 (MQ=255)
                       ggTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAt                                                                                                                               >  1:1542020/1‑150 (MQ=255)
                       ggTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTAATTAt                                                                                                                               >  1:598343/1‑150 (MQ=255)
                          aaCAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTAtctc                                                                                                                            >  2:292200/1‑150 (MQ=255)
                           aCAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATctct                                                                                                                           <  1:909238/150‑1 (MQ=255)
                            cAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTg                                                                                                                          <  1:1247001/150‑1 (MQ=255)
                               ttGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTa                                                                                                                       >  1:1575968/1‑150 (MQ=255)
                                tGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTaa                                                                                                                      <  2:677966/150‑1 (MQ=255)
                                 gCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAAc                                                                                                                     >  2:788852/1‑150 (MQ=255)
                                  cGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACt                                                                                                                    >  1:238818/1‑150 (MQ=255)
                                                  ggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGc                                                                                                    <  1:1631300/150‑1 (MQ=255)
                                                  ggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGc                                                                                                    <  1:1366550/150‑1 (MQ=255)
                                                  ggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCAGGTAACTGGAGATAATTTTTAGc                                                                                                    >  2:1124033/1‑150 (MQ=255)
                                                                  gagaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGc                                                                                    >  2:197606/1‑150 (MQ=255)
                                                                  gagaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCACGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTga                                                                                    >  1:744086/1‑149 (MQ=255)
                                                                     aaCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTaa                                                                                 <  2:1093716/150‑1 (MQ=255)
                                                                     aaCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTaa                                                                                 <  1:392066/150‑1 (MQ=255)
                                                                       cGCATTGAGATATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtg                                                                               <  1:426497/150‑1 (MQ=255)
                                                                           ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGAt                                                                           >  1:1125598/1‑150 (MQ=255)
                                                                            tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                          >  1:1269809/1‑150 (MQ=255)
                                                                            tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCg                                                                         <  2:19283/150‑1 (MQ=255)
                                                                             gagagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTgg                                                                                                                         >  1:1556118/1‑102 (MQ=255)
                                                                              agagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGt                                                                        >  2:1001990/1‑150 (MQ=255)
                                                                                agaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTaa                                                                      <  1:1527530/150‑1 (MQ=255)
                                                                                 gaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAAt                                                                     >  1:386525/1‑150 (MQ=255)
                                                                                   ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATca                                                                   >  2:354992/1‑150 (MQ=255)
                                                                                          gACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGtgtg                                                            >  1:545648/1‑150 (MQ=255)
                                                                                          gACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGtgtg                                                            <  2:1514863/150‑1 (MQ=255)
                                                                                          gACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGtgtg                                                            <  1:225531/150‑1 (MQ=255)
                                                                                                 ttAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCg                                                     >  2:1360294/1‑150 (MQ=255)
                                                                                                  tAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGt                                                    <  1:560401/150‑1 (MQ=255)
                                                                                                   aGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa                                                   <  2:238818/150‑1 (MQ=255)
                                                                                                   aGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa                                                   <  2:576560/150‑1 (MQ=255)
                                                                                                    gCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAt                                                  >  2:149416/1‑150 (MQ=255)
                                                                                                    gCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAt                                                  <  1:274224/150‑1 (MQ=255)
                                                                                                     cTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAtt                                                 <  2:1254167/150‑1 (MQ=255)
                                                                                                      tCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTg                                                <  2:801425/150‑1 (MQ=255)
                                                                                                       ccAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGt                                               >  2:1126793/1‑150 (MQ=255)
                                                                                                         aaTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCa                                             >  2:140410/1‑150 (MQ=255)
                                                                                                            ttCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaat                                          <  1:201652/150‑1 (MQ=255)
                                                                                                              ccGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAACGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaataa                                        >  1:18061/1‑150 (MQ=255)
                                                                                                               cGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaataat                                       >  1:78560/1‑150 (MQ=255)
                                                                                                               cGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaataat                                       <  1:1113499/150‑1 (MQ=255)
                                                                                                                gTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtt                                      <  1:819171/150‑1 (MQ=255)
                                                                                                                  aaGGCAAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAttt                                     <  2:421823/150‑1 (MQ=255)
                                                                                                                  aaGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtttt                                    <  2:1337054/150‑1 (MQ=255)
                                                                                                                   aGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTa                                   <  2:1021414/150‑1 (MQ=255)
                                                                                                                   aGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTa                                   <  2:156048/150‑1 (MQ=255)
                                                                                                                   aGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTa                                   <  2:1557264/150‑1 (MQ=255)
                                                                                                                    ggCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaa                                  >  1:1031691/1‑150 (MQ=255)
                                                                                                                    ggCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTAAAATAATTTTaa                                  >  2:1312310/1‑150 (MQ=255)
                                                                                                                      cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaa                                <  2:9780/150‑1 (MQ=255)
                                                                                                                      cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATATCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATATCTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaa                                >  2:223256/1‑150 (MQ=255)
                                                                                                                        aaaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCg                                                                         <  2:1184169/108‑1 (MQ=255)
                                                                                                                        aaaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCg                                                                         >  1:1184169/1‑108 (MQ=255)
                                                                                                                        aaaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaaa                               <  1:429664/150‑1 (MQ=255)
                                                                                                                         aaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAg                              <  1:807735/150‑1 (MQ=255)
                                                                                                                           aaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCa                            <  2:1120421/150‑1 (MQ=255)
                                                                                                                             aaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGa                          <  2:525091/150‑1 (MQ=255)
                                                                                                                               gagaaagGTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAg                        >  1:1505830/2‑150 (MQ=255)
                                                                                                                               aagaaagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAg                        <  1:1310761/150‑1 (MQ=255)
                                                                                                                                agaaagGTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGc                       >  1:1472228/1‑150 (MQ=255)
                                                                                                                                 gaaagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCg                      >  1:1202806/1‑150 (MQ=255)
                                                                                                                                  aaagGTAACTTCCCCCGCGAGAGTGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGt                     <  2:404234/150‑1 (MQ=255)
                                                                                                                                   aagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTg                    <  1:1001990/150‑1 (MQ=255)
                                                                                                                                    agGTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAGGCGTGa                   >  2:722870/1‑150 (MQ=39)
                                                                                                                                     ggTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  2:1437055/1‑150 (MQ=255)
                                                                                                                                     ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  1:1584877/1‑150 (MQ=255)
                                                                                                                                     ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  1:303626/1‑150 (MQ=255)
                                                                                                                                     ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  1:303625/1‑150 (MQ=255)
                                                                                                                                     ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  1:878124/1‑150 (MQ=255)
                                                                                                                                      gTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAtt                 >  1:1056855/1‑150 (MQ=255)
                                                                                                                                          cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActc             >  2:58544/1‑150 (MQ=255)
                                                                                                                                           ttCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActct            <  2:1269809/150‑1 (MQ=255)
                                                                                                                                            tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtttt                                    <  2:1222395/125‑1 (MQ=255)
                                                                                                                                            tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtttt                                    >  1:1222395/1‑125 (MQ=255)
                                                                                                                                            tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACtctc           <  2:598343/150‑1 (MQ=255)
                                                                                                                                            tCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACtctc           <  1:1124033/150‑1 (MQ=255)
                                                                                                                                               cccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcaca        >  2:51004/1‑150 (MQ=255)
                                                                                                                                                ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcacac       <  1:1331069/150‑1 (MQ=255)
                                                                                                                                                 cgcgAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcacaca      <  2:78560/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAACAC  >  NC_000913/2900179‑2900476

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