Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,156,071 C→G pseudogene (658/2179 nt) rlmKL → bifunctional 23S rRNA (guanine(2069)‑N(7))‑methyltransferase RlmK/23S rRNA (guanine(2445)‑N(2))‑methyltransferase RlmL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,156,0710CG99.8% 1830.8 / NA 462pseudogene (658/2179 nt)rlmKLbifunctional 23S rRNA (guanine(2069)‑N(7))‑methyltransferase RlmK/23S rRNA (guanine(2445)‑N(2))‑methyltransferase RlmL
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (217/244);  total (217/245)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 462 total aligned reads displayed.

GCAGGGCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTACG  >  NC_008463/2155923‑2156219
                                                                                                                                                    |                                                                                                                                                    
     gcgccgcgccGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATcgc                                                                                                                                               <  2:6867390/151‑1 (MQ=255)
         cgcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCg                                                                                                                                           >  2:7547412/1‑151 (MQ=255)
          gcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGa                                                                                                                                          >  2:4520577/1‑151 (MQ=255)
           cgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGaa                                                                                                                                         >  2:10941348/1‑151 (MQ=255)
           cgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGaa                                                                                                                                         <  2:3274809/151‑1 (MQ=255)
            gccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAc                                                                                                                                        <  1:9185788/151‑1 (MQ=255)
            gccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAc                                                                                                                                        >  2:10578838/1‑151 (MQ=255)
             ccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCTCCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAcc                                                                                                                                       >  2:6512470/1‑151 (MQ=255)
             ccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAcc                                                                                                                                       <  1:5453663/151‑1 (MQ=255)
               gcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTc                                                                                                                                     <  1:2903075/151‑1 (MQ=255)
                cTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCa                                                                                                                                    >  1:10561019/1‑151 (MQ=255)
                cTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCa                                                                                                                                    <  2:3453716/151‑1 (MQ=255)
                      gagaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGcgcg                                                                                                                              <  1:6627596/151‑1 (MQ=255)
                        gaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAg                                                                                                                            <  1:1901264/151‑1 (MQ=255)
                        gaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAg                                                                                                                            <  1:3840650/151‑1 (MQ=255)
                         aaCCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAgc                                                                                                                           >  2:4696130/1‑151 (MQ=255)
                             tGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCTGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTg                                                                                                                       <  2:4099033/151‑1 (MQ=255)
                             tGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTg                                                                                                                       <  1:4724310/151‑1 (MQ=255)
                               gcggcggcgGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTggg                                                                                                                     <  2:8737604/151‑1 (MQ=255)
                                cggcggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTgggg                                                                                                                    >  1:8862296/1‑151 (MQ=255)
                                cggcggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTgggg                                                                                                                    >  2:6329636/1‑151 (MQ=255)
                                    ggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTc                                                                                                                >  2:5375559/1‑151 (MQ=255)
                                     gcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCa                                                                                                               <  1:8650992/151‑1 (MQ=255)
                                      cggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAg                                                                                                              >  1:8357315/1‑151 (MQ=255)
                                         tgtTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCaa                                                                                                           >  1:9751675/1‑151 (MQ=255)
                                            tGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTg                                                                                                        <  1:4188199/151‑1 (MQ=255)
                                            tGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTg                                                                                                        <  1:7701315/151‑1 (MQ=255)
                                              aTCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGc                                                                                                      <  1:9398594/151‑1 (MQ=255)
                                               tCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCt                                                                                                     <  2:8588307/151‑1 (MQ=255)
                                                 cGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGTTATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTgg                                                                                                   <  1:11069014/151‑1 (MQ=255)
                                                 cGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTgg                                                                                                   <  1:570275/151‑1 (MQ=255)
                                                 cGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTgg                                                                                                   <  2:1632792/151‑1 (MQ=255)
                                                   tGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTgggg                                                                                                 <  2:154142/151‑1 (MQ=255)
                                                     ccGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCa                                                                                               >  2:7818170/1‑151 (MQ=255)
                                                     ccGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCa                                                                                               <  2:6665243/151‑1 (MQ=255)
                                                      cGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCAc                                                                                              <  1:7411422/151‑1 (MQ=255)
                                                       ggCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACg                                                                                             >  2:10069372/1‑151 (MQ=255)
                                                       ggCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACg                                                                                             <  2:381431/151‑1 (MQ=255)
                                                         cTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTg                                                                                           >  1:4933719/1‑151 (MQ=255)
                                                          tGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGc                                                                                          >  2:10806183/1‑151 (MQ=255)
                                                          tGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGc                                                                                          >  1:4342546/1‑151 (MQ=255)
                                                                aaGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGgcgc                                                                                    >  2:2773516/1‑151 (MQ=255)
                                                                 aGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCt                                                                                   <  1:6519180/151‑1 (MQ=255)
                                                                 aGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCt                                                                                   <  2:5671766/151‑1 (MQ=255)
                                                                 aGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCt                                                                                   >  1:2664534/1‑151 (MQ=255)
                                                                      gccgccGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACTTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGgc                                                                              <  2:1161112/151‑1 (MQ=255)
                                                                      gccgccGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGgc                                                                              >  1:3327307/1‑151 (MQ=255)
                                                                       ccgccgAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGgcg                                                                             <  2:9410559/151‑1 (MQ=255)
                                                                           cgAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAg                                                                         <  1:10162479/151‑1 (MQ=255)
                                                                            gAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGc                                                                        >  1:9118573/1‑151 (MQ=255)
                                                                              gggCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTc                                                                      <  1:9948812/151‑1 (MQ=255)
                                                                                gCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCa                                                                    <  2:7961174/151‑1 (MQ=255)
                                                                                 cGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCAc                                                                   <  1:6250985/151‑1 (MQ=255)
                                                                                  ggTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACg                                                                  >  1:7942877/1‑151 (MQ=255)
                                                                                  ggTGCCCTGGCCGACCCGATGGGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACg                                                                  >  2:10893363/1‑151 (MQ=255)
                                                                                   gTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGa                                                                 >  2:3991403/1‑151 (MQ=255)
                                                                                    tGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAg                                                                <  1:5237043/151‑1 (MQ=255)
                                                                                     gCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAgg                                                               <  1:1206263/151‑1 (MQ=255)
                                                                                     gCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAgg                                                               >  2:8863090/1‑151 (MQ=255)
                                                                                      cccTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGa                                                              <  1:7453166/151‑1 (MQ=255)
                                                                                       ccTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGaa                                                             <  2:6085866/151‑1 (MQ=255)
                                                                                       ccTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGaa                                                             <  2:8173354/151‑1 (MQ=255)
                                                                                         tGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGc                                                           >  1:8117301/1‑151 (MQ=255)
                                                                                         tGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGc                                                           <  1:1416454/151‑1 (MQ=255)
                                                                                           gCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGccc                                                         <  2:7265805/151‑1 (MQ=255)
                                                                                           gCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGccc                                                         >  2:2479812/1‑151 (MQ=255)
                                                                                            ccGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCca                                                        <  1:2536104/151‑1 (MQ=255)
                                                                                             cGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcag                                                       >  1:711911/1‑151 (MQ=255)
                                                                                             cGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcag                                                       >  1:1311501/1‑151 (MQ=255)
                                                                                                cccGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagcag                                                    >  2:1289596/1‑151 (MQ=255)
                                                                                                cccGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagcag                                                    >  1:5133391/1‑151 (MQ=255)
                                                                                                  cGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcg                                                  <  1:2321710/151‑1 (MQ=255)
                                                                                                   gATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgc                                                 >  1:10307827/1‑151 (MQ=255)
                                                                                                    aTGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgcg                                                >  1:817600/1‑151 (MQ=255)
                                                                                                     tgtgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgcgc                                               >  1:7159271/1‑151 (MQ=255)
                                                                                                      gtgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGcgcgcg                                              <  1:1801919/151‑1 (MQ=255)
                                                                                                      gtgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGcgcgcg                                              >  2:10435929/1‑151 (MQ=255)
                                                                                                        gtgGCGTCGGCACCTTCCTGGTCGTGGCGGCGCTGATGGCTGCCGATATCGCGCCGTACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCATGCTCCACGAGGAAGCCCAGCAGCGCgcggc                                            <  1:1702065/151‑1 (MQ=255)
                                                                                                        gtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcggc                                            <  2:9030164/151‑1 (MQ=255)
                                                                                                         tgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcggcg                                           <  2:2425872/151‑1 (MQ=255)
                                                                                                            cGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggcc                                        <  2:9483123/151‑1 (MQ=255)
                                                                                                              tCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggccgg                                      <  2:8016932/151‑1 (MQ=255)
                                                                                                                ggCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCggccggcc                                    <  2:9670612/151‑1 (MQ=255)
                                                                                                                 gCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCt                                   <  2:4904518/151‑1 (MQ=255)
                                                                                                                  cACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTg                                  >  1:6142476/1‑151 (MQ=255)
                                                                                                                    ccTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGc                                >  2:256623/1‑151 (MQ=255)
                                                                                                                     cTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGcc                               <  1:10020212/151‑1 (MQ=255)
                                                                                                                      ttCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCa                              >  1:10868038/1‑151 (MQ=255)
                                                                                                                             tCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgcc                       >  2:522493/1‑149 (MQ=255)
                                                                                                                                aGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGtt                   <  1:7852434/151‑1 (MQ=255)
                                                                                                                                    ggcgCTGATGGCGGCCGATATCGCGCCGCACCTCAGGCGCGAGCGCTGGGGCGTCCGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTgg               <  1:7602808/151‑1 (MQ=255)
                                                                                                                                         tGATGGCTGCCGATATCGCGCCGACCCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGGGGCGCAAGCTCCTCGAGGAAGCCCAGCCGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTcg          >  2:2149361/1‑151 (MQ=255)
                                                                                                                                          gATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGACGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTcgc         >  2:8494710/1‑151 (MQ=255)
                                                                                                                                             ggCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAATCCCAGCAGCGCGCGGCGGCCGGGCTGGCCAGGCCGCCGTTGTGGATTAGCGGc      >  1:6683608/1‑151 (MQ=255)
                                                                                                                                              gCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCt     <  1:10738773/151‑1 (MQ=255)
                                                                                                                                               cTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTa    <  1:9605687/151‑1 (MQ=255)
                                                                                                                                               cTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTa    <  1:5969835/151‑1 (MQ=255)
                                                                                                                                                tGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTAc   <  2:11202884/151‑1 (MQ=255)
                                                                                                                                                 gCCGATATCGCGCCGACCCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACg  <  2:5972070/151‑1 (MQ=255)
                                                                                                                                                 gCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACg  <  1:6893285/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCAGGGCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTACG  >  NC_008463/2155923‑2156219

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