Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 843,196 T→G 16.0% F289V (TTT→GTT)  rlmF → 23S rRNA m(6)A1618 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,1960TG16.0% 1182.3 / 60.1 492F289V (TTT→GTT) rlmF23S rRNA m(6)A1618 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (127/273);  new base G (90/2);  total (217/275)
Fisher's exact test for biased strand distribution p-value = 4.46e-35
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 504 total aligned reads displayed.

AGAGCAAAGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACG  >  NC_000913/843047‑843345
                                                                                                                                                     |                                                                                                                                                       
agagCAAAGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCt                                                                                                                                                         <  1:4135913/150‑1 (MQ=255)
    cAAAGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATTGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTAt                                                                                                                                                    >  1:1091678/1‑150 (MQ=255)
    cAAAGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTGGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTAt                                                                                                                                                    >  2:6404365/1‑150 (MQ=255)
       aGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGc                                                                                                                                                 <  2:7465432/150‑1 (MQ=255)
        ggCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGcc                                                                                                                                                <  1:5008456/150‑1 (MQ=255)
         gCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGccg                                                                                                                                               >  1:6691127/1‑149 (MQ=255)
         gCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGccg                                                                                                                                               >  1:5154824/1‑149 (MQ=255)
           ttCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTgg                                                                                                                                             >  2:5681679/1‑150 (MQ=255)
             cgcgAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGAc                                                                                                                                           <  1:8309416/150‑1 (MQ=255)
             cgcgAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGAc                                                                                                                                           <  1:10024559/150‑1 (MQ=255)
               cgAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCt                                                                                                                                         <  1:164312/150‑1 (MQ=255)
               cgAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCGGGAcca                                                                                                                                         >  2:1463043/1‑149 (MQ=255)
               cgAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTAATGACTGGAcca                                                                                                                                         >  2:7569574/1‑149 (MQ=255)
                 aaGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTGTTGGCGGGACCttt                                                                                                                                       >  2:5578055/1‑150 (MQ=255)
                  aGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGGTATTGCCGGGACCTTTa                                                                                                                                      >  2:5202867/1‑150 (MQ=255)
                   gCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTAt                                                                                                                                     <  1:643573/150‑1 (MQ=255)
                   gCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCGGGGCCATTAt                                                                                                                                     >  2:8867320/1‑150 (MQ=255)
                     aGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGa                                                                                                                                   <  2:421640/150‑1 (MQ=255)
                      ggTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGaa                                                                                                                                  <  2:6613143/150‑1 (MQ=255)
                             tgGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGag                                                                                                                           <  1:2950072/150‑1 (MQ=255)
                             tgGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGag                                                                                                                           <  2:5608662/150‑1 (MQ=255)
                             tgGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCAGGAGCGTTATTAacgtagag                                                                                                                           >  1:6355718/1‑145 (MQ=255)
                               gTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTTTTTTCTGGACCGTTATTATagacgatga                                                                                                                         >  1:319303/1‑147 (MQ=255)
                               gTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCGGGACCTTTATGTACGACGagca                                                                                                                         >  1:7913776/1‑150 (MQ=255)
                               gTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGagca                                                                                                                         <  1:3717344/150‑1 (MQ=255)
                               gTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTGTATTTCCTGGGCCGTTATTAAAGAACagca                                                                                                                         >  1:1922843/1‑150 (MQ=255)
                               gTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCGGGACCATTATTAACGACGagcc                                                                                                                         >  2:4899096/1‑149 (MQ=255)
                                tttACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGagcag                                                                                                                        <  2:2012950/150‑1 (MQ=255)
                                      tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGGCCGTTATGAACCACGAGCCCcgccgc                                                                                                                 >  1:7347107/1‑150 (MQ=255)
                                      tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCAGGACCTTTATGAAAGACGAGCAGcgccgc                                                                                                                 >  2:8244020/1‑150 (MQ=255)
                                      tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCAGGACCTTTATGAACGAAGAGCAGcgccgc                                                                                                                 >  2:620349/1‑150 (MQ=255)
                                       cACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCgccgcc                                                                                                                <  1:5520877/150‑1 (MQ=255)
                                         cTGGTATCTCGTGGTGAAAACTTACCGCCGGTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCAGGACCATTATTAAACTCCAGCAGCgccgccct                                                                                                              >  2:5153041/1‑148 (MQ=255)
                                           ggTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGttt                                                                                                            <  2:3699676/150‑1 (MQ=255)
                                           ggTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGttt                                                                                                            >  2:5846069/1‑150 (MQ=255)
                                           ggTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGttt                                                                                                            <  2:9590519/150‑1 (MQ=255)
                                            gTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAATAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGGCGGGGCCGTTATGTACAGAGAGGAGCGCCGCCGattt                                                                                                           >  1:8806055/1‑150 (MQ=255)
                                            gTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGtttt                                                                                                           <  2:339989/150‑1 (MQ=255)
                                            gTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGAAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGtttt                                                                                                           <  1:4363979/150‑1 (MQ=255)
                                              aTCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGt                                                                                                         <  1:10175869/150‑1 (MQ=255)
                                              aTCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGt                                                                                                         <  1:2215863/150‑1 (MQ=255)
                                               tctcGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTc                                                                                                        <  1:1193886/150‑1 (MQ=255)
                                                ctcGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCa                                                                                                       <  1:2566289/150‑1 (MQ=255)
                                                 tcGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTGCATGGCCGATATTAACGACGAGCCGCGCCCCCCGtttttaaa                                                                                                       >  1:5286968/1‑146 (MQ=255)
                                                 tcGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTAATTCCGGGCCCTTTATTAAAGACGAGCAGCGCCCCCGGTTTGTCaa                                                                                                      >  2:8642674/1‑150 (MQ=255)
                                                   gtggtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATc                                                                                                    <  2:7154309/150‑1 (MQ=255)
                                                       tgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGACGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTAGTGCCGGGACCGTTATGAACGACGAGCAGCGCCGCCGGTGTGTCAAACGccc                                                                                                >  2:7746889/1‑149 (MQ=255)
                                                       tgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCa                                                                                                <  1:2474382/150‑1 (MQ=255)
                                                         aaaaCTTACCGCCGGTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTAATTCCTCGACCTATATTAAAGACCACCACCGCACCCCGTTAGTCAATCCCCCGc                                                                                              >  2:8046699/1‑150 (MQ=255)
                                                              ttACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCAGGACCGTTATTAAACAAGATCAGCGCCCCCCGTTTTTAAATAGCCacccgtta                                                                                         >  2:6908533/1‑143 (MQ=255)
                                                               tACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTaaa                                                                                        <  1:2272011/150‑1 (MQ=255)
                                                               tACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTaaa                                                                                        <  2:4005500/150‑1 (MQ=255)
                                                                 ccgccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGc                                                                                      <  2:5864661/150‑1 (MQ=255)
                                                                   gccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGt                                                                                    <  2:4518997/150‑1 (MQ=255)
                                                                   gccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTGCGGGACCATTATGAACGACGTGCAGCGCCCCCGTTTTTTCAATCGCCAGCGGTAAAACGt                                                                                    >  1:8551017/1‑150 (MQ=255)
                                                                    ccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGtt                                                                                   <  2:4338691/150‑1 (MQ=255)
                                                                    ccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCGGGACCATTATTAAACAACAGCAGCCCCCCCCGTTTTTCAAACCCCCCCGTTAAAACggt                                                                                   >  2:1745078/1‑148 (MQ=255)
                                                                     cgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTg                                                                                  <  2:4578625/150‑1 (MQ=255)
                                                                           aTCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTccc                                                                            <  2:7028835/150‑1 (MQ=255)
                                                                              gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTGCAGGAACATTTAGAAAGACAAGCACCCCCCCCCGTTTTTGAAAAGCCACCCTTAAAACCTTGGTTcccccc                                                                         >  2:978442/1‑149 (MQ=255)
                                                                                gCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCcggc                                                                       <  2:10146708/150‑1 (MQ=255)
                                                                                  ccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCcggcgg                                                                     <  2:8605381/150‑1 (MQ=255)
                                                                                  ccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCTTTATTTAAGACGAGCAGCGCCCCCCTTTTGTCAATCGCCAGCGGTAAAGCGTTGGTTCCCCccgcgg                                                                     >  1:1349757/1‑150 (MQ=255)
                                                                                  ccTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGGTATTTCCTGGGCCATTATTAAAGACGAGCAGCGCCCCCGGTTTGTGAAAAGCCCGCCGTAAAGCGGTGGGTGCCCccgcgg                                                                     >  1:8901132/1‑150 (MQ=255)
                                                                                     gacggacgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTGCGGGGCCATTATTAAACAACAGCACCCCCCCCCGTTTGTCAAACGCCACCCTTAAAGCGGTGGTTCCCCCCGCCGTaa                                                                  >  2:2744993/1‑150 (MQ=255)
                                                                                         gacgTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGt                                                              <  1:9943250/150‑1 (MQ=255)
                                                                                            gTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCGGGACCGTTATGTAAGACGAGCAGCGCCGCCGGTTTGTGAATAGCCACCCGTAAAAAGGTTGTTGCCCCCGCGGGAAAGGGTCa                                                           >  1:3411813/1‑150 (MQ=255)
                                                                                                 cgcgGTGATTTTGTTTAAAAAATAGCTGGCCCAGTGGCAAAAGAAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGGATCGCCGCCGTTTTGTCAATCGCCAGAGTTAAAGCGTTGGTTACCCCTGCGTTTACGGTTCAGCCgg                                                      <  1:6022462/150‑1 (MQ=255)
                                                                                                    ggTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCATTATGAAAGACGAGCAGCGCCGCCGGTTTGTCAATAGCCACCGTTAAAGCGTTGGTTCCCCCCGCGGTAACCGTTCAGCCCGGGt                                                   >  2:10088609/1‑150 (MQ=255)
                                                                                                        aaGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTcgcg                                               <  1:8022813/150‑1 (MQ=255)
                                                                                                         aGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCgcgc                                              <  1:2836847/150‑1 (MQ=255)
                                                                                                               ttAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAg                                        <  2:791089/150‑1 (MQ=255)
                                                                                                                 aaaaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGACCGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgca                                      <  2:9905043/150‑1 (MQ=255)
                                                                                                                 aaaaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgca                                      <  1:3989765/150‑1 (MQ=255)
                                                                                                                 aaaaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgca                                      <  2:4711859/150‑1 (MQ=255)
                                                                                                                   aaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCTGGACCATTATGAACCACGAGCAGCGCCGCCCTTTTGTCAAAAGCCACCCGTAAAACGTTGGGTCCCCCCGCGGGAACCGGTCAGCCCGGGGCGCGCCCCGCgcagg                                    >  2:1029954/1‑150 (MQ=255)
                                                                                                                     aaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTGCGGGACCGTTATGAACGCCCAGCACCGCCCCCGGGTTGTCAATAGCCACCGGTAAAGCGGTGGGTCCCCCCGCGGGAACGGGTCAGCCCGGGGCGCGCCCCGAgccggcc                                  >  2:2956195/1‑149 (MQ=255)
                                                                                                                         gaTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAc                              <  1:5782201/150‑1 (MQ=255)
                                                                                                                         gaTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAc                              <  1:8483115/150‑1 (MQ=255)
                                                                                                                          aTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCGTGACCTTTATGAAACACCAGCAGCGCCCCCGTTTTTTCAATACCCAGCGTTAAAGCGTGTGTTACCCCGGCGGTAACCGGTCAGCCCGGGTCGCGCCCGGAGCAGGCAGCAcc                             >  2:3247401/1‑150 (MQ=255)
                                                                                                                          aTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAcc                             <  2:6445374/150‑1 (MQ=255)
                                                                                                                          aTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCAcc                             <  2:9244384/150‑1 (MQ=255)
                                                                                                                          aTGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCATTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATAGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGGAACGGTTCAGCCCGGGTCGCGCCCGGAGCAGGCAGcccc                             >  2:2553908/1‑150 (MQ=255)
                                                                                                                              cccAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGAt                         <  1:750298/150‑1 (MQ=255)
                                                                                                                               ccAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATa                        <  2:7309190/150‑1 (MQ=255)
                                                                                                                                 aGGGGCAAAAGCAGAGTCGCGTTATTGCCAGGACCTTTATGAAACACGAGCAGCGCCCCCCGTTTGTCAATAGCCAGCGTTAAATCCTTGGGTCCCCCCGCCGGACCCGTTCAGCCCGGGGCCCGCCCCGAGCACGCCGCACCAGATAAg                      >  2:6147552/1‑150 (MQ=255)
                                                                                                                                  ggggCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGt                     <  1:4659885/150‑1 (MQ=255)
                                                                                                                                    ggCAAAAGCAGAGTCGCTTTATTTCCGGGCCCTTTATTAAACACGAGCAGCCCCCCCCGTTTTTCAAACGCCCGCGGTTAAGAGTGGGTTCCCCCCGCGGGAACCGTTCACCCCGGGTCGCCCCCCGCGCCCGCAGCACCAGATCAGTCt                   >  2:9686721/1‑150 (MQ=255)
                                                                                                                                    ggCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCt                   <  2:5011443/150‑1 (MQ=255)
                                                                                                                                     gCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTg                  <  1:9377266/150‑1 (MQ=255)
                                                                                                                                     gCAAAAGCAGAGTCGCGTTATTTCCGGGACCTTTATTAACGACGAGCAGCGCCGCCGGTTTTTCAAACGCCAGCGGTAAAGCGTTGGGTGCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTg                  >  2:2411421/1‑150 (MQ=255)
                                                                                                                                       aaaaGCAGAGTCGCTTTATTTCCGGGACCATTATTAAACACGAGCAGCGCCCCCGTTTTTTCAAACGCCCGCGTTAAAGCGTTGGTTCCCCCCGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCACGCAGCACCAGATAAGTCAGCa                >  2:3035120/1‑150 (MQ=255)
                                                                                                                                         aaGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACt              <  2:8769851/150‑1 (MQ=255)
                                                                                                                                         aaGCAGAGTCGCTTTAGTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACt              <  1:9092738/150‑1 (MQ=255)
                                                                                                                                               agTCGCTCTTTTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTATTGTCAATCGCCAGCGATAAAGCGTTGGGTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgcg        <  1:6932463/150‑1 (MQ=255)
                                                                                                                                                gTCGCTTTAGTGCCTGGACCTTTCTGACCGCCGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGgcgcgc       <  2:8281282/150‑1 (MQ=255)
                                                                                                                                                 tCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGcgcgcg      >  2:6014464/1‑150 (MQ=255)
                                                                                                                                                   gCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa    <  2:6849317/150‑1 (MQ=255)
                                                                                                                                                   gCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa    >  2:5369809/1‑150 (MQ=255)
                                                                                                                                                   gCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGaa    >  1:5093829/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
AGAGCAAAGGCTTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTGCACTGGCGCGCGAACG  >  NC_000913/843047‑843345

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