Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,642,358 T→A 11.6% intergenic (+27/+22) dtpB → / ← rsmJ dipeptide/tripeptide:H(+) symporter DtpB/16S rRNA m(2)G1516 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,642,3580TA11.6% 319.9 / 28.0 120intergenic (+27/+22)dtpB/rsmJdipeptide/tripeptide:H(+) symporter DtpB/16S rRNA m(2)G1516 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (50/56);  new base A (5/9);  total (55/65)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 120 total aligned reads displayed.

TTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCC  >  NC_000913/3642212‑3642503
                                                                                                                                                  |                                                                                                                                                 
ttGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATaa                                                                                                                                                <  1:1862519/150‑2 (MQ=255)
   cccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGc                                                                                                                                             >  2:379381/1‑150 (MQ=255)
    ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCt                                                                                                                                            <  1:1995895/150‑1 (MQ=255)
    ccGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCt                                                                                                                                            >  2:1053419/1‑150 (MQ=255)
      gTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAc                                                                                                                                          <  2:60076/150‑1 (MQ=255)
          ataCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTACCCTg                                                                                                                                      <  1:1788388/150‑1 (MQ=255)
          ataCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTg                                                                                                                                      >  2:1670824/1‑150 (MQ=255)
           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTa                                                                                                                                           >  1:271768/1‑144 (MQ=255)
           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTa                                                                                                                                           <  2:271768/144‑1 (MQ=255)
           taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTACCCTgc                                                                                                                                     <  1:1792249/150‑2 (MQ=255)
             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                   >  1:198142/1‑150 (MQ=255)
             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                   >  2:1881928/1‑150 (MQ=255)
             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttt                                                                                                                                   >  2:58469/1‑150 (MQ=255)
             ccAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgctt                                                                                                                                   <  1:121015/150‑1 (MQ=255)
               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttta                                                                                                                                 >  2:494558/1‑150 (MQ=255)
               aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttta                                                                                                                                 >  1:1775302/1‑150 (MQ=255)
                aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTACCCTgctttct                                                                                                                                <  1:379381/150‑3 (MQ=255)
                aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgctttct                                                                                                                                <  1:615867/150‑3 (MQ=255)
                 cGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgctttctg                                                                                                                               <  1:1670824/150‑4 (MQ=255)
                  gtgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtt                                                                                                                              >  1:514213/1‑150 (MQ=255)
                    gtTTGGTAAGATTGGTCTGGTCACGCTAGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttt                                                                                                                            <  1:1840401/150‑1 (MQ=255)
                     tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAACAGGGTAGCGATAACGCTACCCTgctttctgcaag                                                                                                                           <  2:1217421/150‑8 (MQ=39)
                       tGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTaca                                                                                                                         <  1:494558/150‑1 (MQ=255)
                        ggTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTacac                                                                                                                        >  2:1192898/1‑150 (MQ=255)
                         gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                       >  1:570379/1‑150 (MQ=255)
                         gTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACt                                                                                                                       >  1:1269678/1‑150 (MQ=255)
                          tAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAATCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTTTCGCTACCCTGTTTTAGTTTTACACTg                                                                                                                      <  2:789513/150‑1 (MQ=255)
                          tAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTg                                                                                                                      <  2:478399/150‑1 (MQ=255)
                            aGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGc                                                                                                                    >  2:746118/1‑150 (MQ=255)
                            aGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGc                                                                                                                    >  1:1561948/1‑150 (MQ=255)
                             gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCg                                                                                                                   <  1:1497101/150‑1 (MQ=255)
                             gATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCg                                                                                                                   <  2:1936947/150‑1 (MQ=255)
                                 ggtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTAcc                                                                                                               <  2:198142/150‑1 (MQ=255)
                                   tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTg                                                                                                             >  1:1697035/1‑150 (MQ=255)
                                    ctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGc                                                                                                            >  2:21330/1‑150 (MQ=255)
                                         cACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGa                                                                                                       >  1:1997540/1‑150 (MQ=255)
                                           cGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGata                                                                                                     >  1:433754/1‑150 (MQ=255)
                                            gCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGatat                                                                                                    >  2:1441238/1‑150 (MQ=255)
                                                     tGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGAt                                                                                           <  2:909588/150‑1 (MQ=255)
                                                     tGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGAt                                                                                           <  1:1161825/150‑1 (MQ=255)
                                                      gCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATg                                                                                          <  2:123138/150‑1 (MQ=255)
                                                       cAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTTCACTGGCGTACCTGCATAGATATCAAAGCGATGc                                                                                         >  1:1047501/1‑150 (MQ=255)
                                                           agtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCtt                                                                                     <  2:503526/150‑1 (MQ=255)
                                                            gtGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCttt                                                                                    >  1:557753/1‑150 (MQ=255)
                                                             tGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAg                                                                                               <  2:347723/138‑1 (MQ=255)
                                                             tGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAg                                                                                               >  1:347723/1‑138 (MQ=255)
                                                              gATGCTGTTGATGGTGGCGTGTCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                  <  1:1310585/102‑1 (MQ=255)
                                                              gATGCTGTTGATGGTGGCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                  >  2:1310585/1‑102 (MQ=255)
                                                              gATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAg                                                                                  <  1:1914367/150‑1 (MQ=255)
                                                                tGCTGTTGATGTTTCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTg                                                                                <  1:1980556/150‑1 (MQ=255)
                                                                tGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTg                                                                                >  1:981521/1‑150 (MQ=255)
                                                                tGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTg                                                                                <  1:21330/150‑1 (MQ=255)
                                                                 gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTgg                                                                               >  1:1114120/1‑149 (MQ=255)
                                                                  cTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGAc                                                                              <  2:1233053/150‑1 (MQ=255)
                                                                   tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACt                                                                             >  2:348496/1‑150 (MQ=255)
                                                                    gTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTa                                                                            >  1:701324/1‑150 (MQ=255)
                                                                      tGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTAcc                                                                          >  1:908561/1‑150 (MQ=255)
                                                                         tGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGGCTACcgcg                                                                       <  2:1114120/150‑1 (MQ=255)
                                                                          ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGt                                                                      <  1:970691/150‑1 (MQ=255)
                                                                             gCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtg                                                                   <  2:1225563/150‑1 (MQ=255)
                                                                             gCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtg                                                                   <  2:823293/150‑1 (MQ=255)
                                                                             gCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtg                                                                   >  1:1162712/1‑150 (MQ=255)
                                                                                   gCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCt                                                                                      <  2:1600766/125‑1 (MQ=255)
                                                                                    cTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggc                                                            >  1:1809097/1‑150 (MQ=255)
                                                                                    cTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtggc                                                            >  1:1816440/1‑150 (MQ=255)
                                                                                      gAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGa                                                          <  2:388087/150‑1 (MQ=255)
                                                                                      gAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGa                                                          <  1:1192898/150‑1 (MQ=255)
                                                                                         aCGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAt                                                       >  1:1666768/1‑150 (MQ=255)
                                                                                          cGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                      >  2:1756826/1‑150 (MQ=255)
                                                                                           gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                     <  2:1083791/150‑1 (MQ=255)
                                                                                           gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                     >  1:485987/1‑150 (MQ=255)
                                                                                            cATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCg                                                    <  1:1207622/150‑1 (MQ=255)
                                                                                            cATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCg                                                    <  1:590312/150‑1 (MQ=255)
                                                                                              tGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGcc                                                  >  1:1543392/1‑150 (MQ=255)
                                                                                                  tGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCTGGCGAATGTCGCCacgcc                                             <  1:1896294/150‑5 (MQ=255)
                                                                                                     gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                           <  2:433754/150‑1 (MQ=255)
                                                                                                     gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                           <  1:161132/150‑1 (MQ=255)
                                                                                                        gCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcg                                        <  2:197040/150‑1 (MQ=255)
                                                                                                         ccGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcgc                                       >  2:934254/1‑150 (MQ=255)
                                                                                                           gtaaaGGCATTAATTATTCTTGCAGAAAGCAGGGTAGGGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGt                                     <  1:867258/148‑1 (MQ=255)
                                                                                                           ggAAAGGGATTTATTTTTCTTGCAGAAAGCATGGTAGCGTTTTCGTTACCCTGTTTTAGTTTTACACTGGCGTACCTGGATAGATATCAATGCGATGCCCTTTAGTGACTTCCGCGTTTGGCGTGGCGACATTCGCCAGCTGTGGGGCGt                                     <  2:1561948/150‑1 (MQ=255)
                                                                                                               aGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCACCGCTGgcgaatgtcg                                 <  1:58469/150‑8 (MQ=255)
                                                                                                                 ccATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCg                               <  2:1256938/150‑1 (MQ=255)
                                                                                                                  cATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTTCACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCgg                              >  2:1721653/1‑150 (MQ=255)
                                                                                                                  cATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCgg                              >  1:1992692/1‑150 (MQ=255)
                                                                                                                   attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                             >  1:148893/1‑150 (MQ=255)
                                                                                                                         tattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCttt                       >  1:1856223/1‑150 (MQ=255)
                                                                                                                              ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTaccac                  <  1:1580188/150‑1 (MQ=255)
                                                                                                                               tGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTaccacc                 <  1:348496/150‑1 (MQ=255)
                                                                                                                                  aGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCAcg              <  2:1697035/150‑1 (MQ=255)
                                                                                                                                     aaGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGt           >  2:656334/1‑150 (MQ=255)
                                                                                                                                      aGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGtt          <  2:557753/150‑1 (MQ=255)
                                                                                                                                       gCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGttt         <  2:1775302/150‑1 (MQ=255)
                                                                                                                                       acAGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGttt         >  1:1005083/2‑150 (MQ=255)
                                                                                                                                        cAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTg        <  2:485987/150‑1 (MQ=255)
                                                                                                                                        cAGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTg        >  2:1772638/1‑150 (MQ=255)
                                                                                                                                         aGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTgg       <  2:138738/150‑1 (MQ=255)
                                                                                                                                           ggTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc     >  2:683729/1‑150 (MQ=255)
                                                                                                                                           ggTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc     <  2:1002343/150‑1 (MQ=255)
                                                                                                                                              aGCGTTATCGCTACCCTGTTTTAGTTTTTCACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGcc  <  1:1721653/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCC  >  NC_000913/3642212‑3642503

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: