Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,156,146 C→A pseudogene (733/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,1460CA99.5% 1804.9 / NA 427pseudogene (733/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/2);  new base A (215/209);  total (216/211)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 427 total aligned reads displayed.

AGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTGGGTGAAGATCT  >  NC_008463/2156000‑2156294
                                                                                                                                                  |                                                                                                                                                    
    cGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCa                                                                                                                                               <  2:4880415/150‑1 (MQ=255)
     ggTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCAc                                                                                                                                              <  1:1234507/150‑1 (MQ=255)
      gTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACg                                                                                                                                             >  1:7387607/1‑150 (MQ=255)
           cTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGaa                                                                                                                                        <  1:2719290/150‑1 (MQ=255)
           cTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGaa                                                                                                                                        <  2:514754/150‑1 (MQ=255)
            tGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAg                                                                                                                                       >  1:6838859/1‑150 (MQ=255)
             ggCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGc                                                                                                                                      >  2:7476396/1‑150 (MQ=255)
                  aCCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagc                                                                                                                                 <  2:7397048/150‑1 (MQ=255)
                   cccGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGACGCCcagca                                                                                                                                >  2:816509/1‑150 (MQ=255)
                    ccGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagcag                                                                                                                               >  2:3762398/1‑150 (MQ=255)
                      gATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcg                                                                                                                             <  2:4024295/150‑1 (MQ=255)
                       aTGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgc                                                                                                                            <  1:746676/150‑1 (MQ=255)
                         gtgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgcgc                                                                                                                          >  1:3045070/1‑150 (MQ=255)
                          tgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGcgcgcg                                                                                                                         <  1:7624741/150‑1 (MQ=255)
                           gtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcgg                                                                                                                        <  1:962553/150‑1 (MQ=255)
                            tgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcggc                                                                                                                       <  2:2065082/150‑1 (MQ=255)
                               cGTCGGCACCTTCCTGGTCGAGGCGGCGCTGCTGGCTGCCGATCTCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCggcggc                                                                                                                    <  2:9151239/150‑1 (MQ=255)
                               cGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCggcggc                                                                                                                    <  2:4784716/150‑1 (MQ=255)
                                   ggCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggccggc                                                                                                                >  1:2379588/1‑150 (MQ=255)
                                     cACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGACGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCt                                                                                                              >  2:4417286/1‑150 (MQ=255)
                                       ccTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTgg                                                                                                            >  1:4493773/1‑150 (MQ=255)
                                        cTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGc                                                                                                           <  2:3886436/150‑1 (MQ=255)
                                        cTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGc                                                                                                           <  2:3725520/150‑1 (MQ=255)
                                          tCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCa                                                                                                         >  1:6516866/1‑150 (MQ=255)
                                          tCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCa                                                                                                         <  2:2917008/150‑1 (MQ=255)
                                           ccTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAg                                                                                                        <  1:409399/150‑1 (MQ=255)
                                           ccTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAg                                                                                                        >  1:1777255/1‑150 (MQ=255)
                                            cTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAgg                                                                                                       >  2:10667969/1‑150 (MQ=255)
                                            cTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAgg                                                                                                       >  2:10068847/1‑150 (MQ=255)
                                              ggTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggcc                                                                                                     <  1:7282876/150‑1 (MQ=255)
                                                 cGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgcc                                                                                                  <  2:2631387/150‑3 (MQ=255)
                                                   aGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgccgt                                                                                                <  1:4830950/150‑5 (MQ=255)
                                                   aGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgccgt                                                                                                >  2:9380128/1‑146 (MQ=255)
                                                   aGGCGGAGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgccgt                                                                                                >  1:1364759/1‑146 (MQ=255)
                                                      cggcgCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTtgt                                                                                            >  1:9503790/1‑150 (MQ=255)
                                                        gcgcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTgg                                                                                          <  2:9811264/150‑1 (MQ=255)
                                                        gcgcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTgg                                                                                          >  1:2505376/1‑150 (MQ=255)
                                                         cgcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGa                                                                                         >  1:2253727/1‑150 (MQ=255)
                                                          gcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGAt                                                                                        <  1:8922468/150‑1 (MQ=255)
                                                             gATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTcg                                                                                     >  1:4608860/1‑150 (MQ=255)
                                                               tGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTcgcg                                                                                   >  1:8565199/1‑150 (MQ=255)
                                                                ggCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCgg                                                                                  <  1:5753137/150‑1 (MQ=255)
                                                                  cTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCt                                                                                >  1:6674796/1‑150 (MQ=255)
                                                                   tGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTa                                                                               <  2:671244/150‑1 (MQ=255)
                                                                    gCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTAc                                                                              <  1:7986087/150‑1 (MQ=255)
                                                                     ccgatatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGGGTGGTTTCGCGGCTACg                                                                             >  1:3347156/4‑150 (MQ=255)
                                                                     ccgatatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCGGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACg                                                                             <  1:10114477/147‑1 (MQ=255)
                                                                      cgatatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCAGGCCAGGCCGCCGGTGTGGAGTCGCGGCAACGa                                                                            >  1:5357620/3‑150 (MQ=255)
                                                                       gatatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAg                                                                           >  1:4305102/2‑150 (MQ=255)
                                                                        ctatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCAGGCGCTGTGGCGCAAGCTCCTCGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAgg                                                                          >  1:4154085/2‑150 (MQ=255)
                                                                        atatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAgg                                                                          >  2:9724258/1‑150 (MQ=255)
                                                                         tatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGc                                                                         <  2:7666614/150‑1 (MQ=255)
                                                                         tatCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGc                                                                         <  2:1787113/150‑1 (MQ=255)
                                                                             gcgcCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGAc                                                                     <  1:10372523/150‑1 (MQ=255)
                                                                               gcCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGAccc                                                                   <  1:7759053/150‑1 (MQ=255)
                                                                               gcCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGAccc                                                                   >  2:8773493/1‑150 (MQ=255)
                                                                                ccGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGATGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGGTTCGCGGCTTCGAGGCCGGCCcg                                                                  >  1:1847340/1‑150 (MQ=255)
                                                                                  gAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCcgcg                                                                >  1:5054756/1‑150 (MQ=255)
                                                                                   aaCCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCgcgc                                                               >  1:9426140/1‑150 (MQ=255)
                                                                                   aaCCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCgcgc                                                               >  2:5726567/1‑150 (MQ=255)
                                                                                     ccTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCt                                                             <  2:7074749/150‑1 (MQ=255)
                                                                                      cTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTg                                                            <  1:8786557/150‑1 (MQ=255)
                                                                                           gcgcGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCa                                                       <  1:8528121/150‑1 (MQ=255)
                                                                                             gcgAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGc                                                     <  1:4012503/150‑1 (MQ=255)
                                                                                             gcgAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGc                                                     <  1:872370/150‑1 (MQ=255)
                                                                                              cgAGCGCTGGGGCTACAGCACTTGGCTGGGGCACGTGCCGGCGCGGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGAGTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGcc                                                    <  2:10566408/150‑1 (MQ=255)
                                                                                                  cgcTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGgcg                                                <  2:8082255/150‑1 (MQ=255)
                                                                                                     tGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGgcgcgc                                             <  1:974369/150‑1 (MQ=255)
                                                                                                     tGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGAGCCAGCCGgcgcgc                                             >  1:3843976/1‑150 (MQ=255)
                                                                                                       gggCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaa                                           <  1:2117661/150‑1 (MQ=255)
                                                                                                       gggCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaa                                           <  1:150049/150‑1 (MQ=255)
                                                                                                         gCTTCAGCACTTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGATGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcacca                                         >  2:1403261/1‑150 (MQ=255)
                                                                                                         gCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCCGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGccaca                                         >  1:3007297/1‑150 (MQ=255)
                                                                                                         gCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaaca                                         >  1:5183133/1‑150 (MQ=255)
                                                                                                         gCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaaca                                         >  1:6065959/1‑150 (MQ=255)
                                                                                                          cTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaacaa                                        >  2:5024000/1‑150 (MQ=255)
                                                                                                          cTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGcaacaa                                        <  2:5588912/150‑1 (MQ=255)
                                                                                                            tCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCACCAata                                      >  1:6144574/1‑150 (MQ=255)
                                                                                                            tCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAata                                      >  2:1685939/1‑150 (MQ=255)
                                                                                                               gCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCg                                   >  2:7842043/1‑150 (MQ=255)
                                                                                                                cAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGa                                  <  1:3592492/150‑1 (MQ=255)
                                                                                                                 aaTTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAg                                 <  2:8638166/150‑1 (MQ=255)
                                                                                                                         ggggCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACCAGCTCGAGCGTGCCgg                         <  2:528501/150‑1 (MQ=255)
                                                                                                                         ggggCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCgg                         >  2:1707417/1‑150 (MQ=255)
                                                                                                                         ggggCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCgg                         >  2:6436445/1‑150 (MQ=255)
                                                                                                                          gggCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGa                        <  2:10392181/150‑1 (MQ=255)
                                                                                                                          gggCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGa                        <  1:10482897/150‑1 (MQ=255)
                                                                                                                               cGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGc                   <  1:7036510/150‑1 (MQ=255)
                                                                                                                               cGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGc                   <  2:7158409/150‑1 (MQ=255)
                                                                                                                                 tGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCg                 >  1:5903926/1‑150 (MQ=255)
                                                                                                                                  gCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGa                <  1:6466661/150‑1 (MQ=255)
                                                                                                                                   ccGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGAc               <  1:9282592/150‑1 (MQ=255)
                                                                                                                                      gcgcTGTGGCGCAAGCTCCACGAGGAAGCCCAGCCGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGGGGATGCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCACCAAGATCGAGCGTGCCGGAGTGGCCGACAgg            >  1:10260601/1‑150 (MQ=255)
                                                                                                                                      gcgcTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTgg            <  2:7803084/150‑1 (MQ=255)
                                                                                                                                        gcTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGCGTGGCCGACTGGGt          >  2:4617142/1‑150 (MQ=255)
                                                                                                                                         cTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTGGGTg         >  2:597161/1‑150 (MQ=255)
                                                                                                                                           gtgGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTGGGTGaa       >  1:2759588/1‑150 (MQ=255)
                                                                                                                                           gtgGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCACGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCGGGCGCGCAACACTATCGAGCGTGCCGGAGTGGCCGACTGGGTGac       >  2:8450244/1‑149 (MQ=255)
                                                                                                                                            tgGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCACCAATATCGAGCGTGCCGGAGTGGCCGACTGGGTGAAg      >  2:1031324/1‑150 (MQ=255)
                                                                                                                                                gcCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTGGGTGAAGATCt  <  2:6823361/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTACGAGGCCGACCCGCGCCTGATCCAGCCGGCGCGCAACAATATCGAGCGTGCCGGAGTGGCCGACTGGGTGAAGATCT  >  NC_008463/2156000‑2156294

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