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.5% 1665.7 / NA 428pseudogene (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/0);  major base G (216/210);  minor base T (0/2);  total (216/212)
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 428 total aligned reads displayed.

CAGGGCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTAC  >  NC_008463/2155924‑2156218
                                                                                                                                                   |                                                                                                                                                   
  gggCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGatat                                                                                                                                                 >  1:5996520/1‑150 (MQ=255)
  gggCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGatat                                                                                                                                                 >  1:3025925/1‑150 (MQ=255)
   ggCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATc                                                                                                                                                >  1:5668292/1‑150 (MQ=255)
   ggCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATc                                                                                                                                                <  1:6696411/150‑1 (MQ=255)
    gcgccgcgccGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATcg                                                                                                                                               <  2:2636428/150‑1 (MQ=255)
    gcgccgcgccGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATcg                                                                                                                                               >  2:5348513/1‑150 (MQ=255)
      gccgcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATcgcg                                                                                                                                             >  2:10775211/1‑150 (MQ=255)
        cgcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGcc                                                                                                                                           <  2:1442429/150‑1 (MQ=255)
        cgcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGcc                                                                                                                                           <  1:5800002/150‑1 (MQ=255)
        cgcgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGcc                                                                                                                                           <  2:7037400/150‑1 (MQ=255)
          cgccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGa                                                                                                                                         <  1:7260754/150‑1 (MQ=255)
           gccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGaa                                                                                                                                        >  2:3389014/1‑150 (MQ=255)
           gccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGaa                                                                                                                                        >  1:2835779/1‑150 (MQ=255)
           gccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGaa                                                                                                                                        <  2:6649347/150‑1 (MQ=255)
            ccgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAc                                                                                                                                       >  1:1663009/1‑150 (MQ=255)
             cgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAcc                                                                                                                                      >  1:5174776/1‑150 (MQ=255)
             cgcTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAAcc                                                                                                                                      >  1:1044528/1‑150 (MQ=255)
               cTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTc                                                                                                                                    <  2:11048153/150‑1 (MQ=255)
                 gAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAg                                                                                                                                  >  2:3029208/1‑150 (MQ=255)
                   aGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGgc                                                                                                                                <  2:1241351/150‑1 (MQ=255)
                    ggAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGgcg                                                                                                                               >  2:8497258/1‑150 (MQ=255)
                      agaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGcgcg                                                                                                                             >  1:5486531/1‑150 (MQ=255)
                       gaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGa                                                                                                                            >  2:1935605/1‑150 (MQ=255)
                       gaACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGa                                                                                                                            >  2:2849630/1‑150 (MQ=255)
                        aaCCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAg                                                                                                                           >  1:4872962/1‑150 (MQ=255)
                        aaCCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAg                                                                                                                           <  1:197676/150‑1 (MQ=255)
                          ccTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAgcg                                                                                                                         >  1:2567013/1‑150 (MQ=255)
                           cTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAgcgc                                                                                                                        <  1:6391674/150‑1 (MQ=255)
                           cTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAgcgc                                                                                                                        <  2:1466213/150‑1 (MQ=255)
                            tGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGCTATCGCGCCGAACCTCAGGCGCGAGCGCt                                                                                                                       <  2:2164926/150‑1 (MQ=255)
                            tGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCt                                                                                                                       <  1:236537/150‑1 (MQ=255)
                            tGGCGGCGGAGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCt                                                                                                                       <  2:894479/150‑1 (MQ=255)
                             ggcggcggcGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTg                                                                                                                      <  1:8477337/150‑1 (MQ=255)
                              gcggcggcgGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCTATATCGCGCCGAACCTCAGGCGCGAGCGCTgg                                                                                                                     <  2:654803/150‑1 (MQ=255)
                              gcggcggcgGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTgg                                                                                                                     <  1:7972249/150‑1 (MQ=255)
                              gcggcggcgGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTgg                                                                                                                     >  2:5593900/1‑150 (MQ=255)
                              gcggcggcgGTGTTGATCCGTGACGGCTGGCCGAAGATCGCCGGCGAGGGCGGGGCGCTGGCCGACCGGATGTGTGGCGTCGGCACCTTCCTGGTCGTGGCGGCGCTGATGGCTGCCGATATCGCGCAGAACCTCAGGCGCGAGCGCTgg                                                                                                                     <  2:4341385/150‑1 (MQ=255)
                               cggcggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGCGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTggg                                                                                                                    >  2:9412117/1‑150 (MQ=255)
                                ggcggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGAGATCGCGCCGAACCTCAGGCGCGAGCGCTgggg                                                                                                                   >  2:7407843/1‑150 (MQ=255)
                                ggcggcggTGTTGATCCGTGACGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTgggg                                                                                                                   >  2:10222553/1‑150 (MQ=255)
                                 gcggcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGc                                                                                                                  >  1:10351966/1‑150 (MQ=255)
                                    gcggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTc                                                                                                               >  2:10502803/1‑150 (MQ=255)
                                      ggTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAg                                                                                                             <  2:6064699/150‑1 (MQ=255)
                                       gtgtTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGc                                                                                                            >  1:10210701/1‑150 (MQ=255)
                                       gtgtTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGc                                                                                                            <  1:9024328/150‑1 (MQ=255)
                                        tgtTGATCCGTGCCGGCTGGCCGAAGCTCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCa                                                                                                           <  2:5848599/150‑1 (MQ=255)
                                        tgtTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCa                                                                                                           >  1:1311528/1‑150 (MQ=255)
                                         gtTGATCCGTGCCGGCTGGCCGCAGCTCTCCGCCGAGGGCGGTGCCCTGGACGACCCGATGTGTGGCGTCGGCACCTTCCGCGTCGAGGCGGCGCTGATGGCGGCCGATAACGCGCCGTACCTCAGGCGCGAGCGCTGGGGCTTCAGCaa                                                                                                          <  2:5042470/150‑1 (MQ=255)
                                         gtTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCaa                                                                                                          <  1:3467747/150‑1 (MQ=255)
                                          ttGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAAt                                                                                                         >  2:8653369/1‑150 (MQ=255)
                                            gATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTg                                                                                                       >  2:832471/1‑150 (MQ=255)
                                              tCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCTCGAGCGCTGGGGCTTCAGCAATTGGc                                                                                                     <  2:4492015/150‑1 (MQ=255)
                                              tCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGc                                                                                                     <  1:1096611/150‑1 (MQ=255)
                                                  tGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTggg                                                                                                 <  2:7785106/150‑1 (MQ=255)
                                                      ggCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCAc                                                                                             <  1:8390552/150‑1 (MQ=255)
                                                      ggCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCAc                                                                                             >  1:7724053/1‑150 (MQ=255)
                                                       gCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACg                                                                                            >  2:10088966/1‑150 (MQ=255)
                                                        cTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGt                                                                                           <  2:9247306/150‑1 (MQ=255)
                                                         tGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTg                                                                                          >  1:2001398/1‑150 (MQ=255)
                                                            ccGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCg                                                                                       >  1:9949261/1‑150 (MQ=255)
                                                             cGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCgg                                                                                      >  1:7409175/1‑150 (MQ=255)
                                                              gAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGgc                                                                                     <  1:3276717/150‑1 (MQ=255)
                                                               aaGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGgcg                                                                                    <  2:7518933/150‑1 (MQ=255)
                                                                          cgAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGgcgcaa                                                                         <  2:7161179/150‑3 (MQ=255)
                                                                             gggCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCt                                                                      >  1:10867555/1‑150 (MQ=255)
                                                                              ggCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTc                                                                     <  2:8030855/150‑1 (MQ=255)
                                                                               gCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTcc                                                                    <  2:2739351/150‑1 (MQ=255)
                                                                                      ccTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGa                                                             >  1:425991/1‑150 (MQ=255)
                                                                                        tGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCACGCTCCACGAGGAAg                                                           >  1:7153566/1‑150 (MQ=255)
                                                                                        tGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAg                                                           <  1:4472707/150‑1 (MQ=255)
                                                                                            cGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCca                                                       <  1:2125511/150‑1 (MQ=255)
                                                                                             gACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCACTTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcag                                                      >  2:2408240/1‑150 (MQ=255)
                                                                                              aCCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagc                                                     <  1:986736/150‑1 (MQ=255)
                                                                                                ccGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagcag                                                   <  2:4384395/150‑1 (MQ=255)
                                                                                                ccGATGTGTGGCGTCGGCACCTGCCTGGTCGTGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCGGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCcagcag                                                   <  2:4149258/150‑1 (MQ=255)
                                                                                                   aTGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgc                                                <  2:7592338/150‑1 (MQ=255)
                                                                                                   aTGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAgcgc                                                <  1:2895518/150‑1 (MQ=255)
                                                                                                      tgtgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGcgcgcg                                             <  1:6953029/150‑1 (MQ=255)
                                                                                                        tgGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcggc                                           <  2:1616252/150‑1 (MQ=255)
                                                                                                         ggCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCgcggcg                                          >  2:9534852/1‑150 (MQ=255)
                                                                                                           cGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCggcggc                                        <  2:824634/150‑1 (MQ=255)
                                                                                                            gTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggcc                                       <  2:4699739/150‑1 (MQ=255)
                                                                                                             tCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggccg                                      >  1:11117311/1‑150 (MQ=255)
                                                                                                             tCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGcggccg                                      >  2:7937634/1‑150 (MQ=255)
                                                                                                                gCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCggccggcc                                   >  2:11081583/1‑150 (MQ=255)
                                                                                                                gCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCggccggcc                                   <  2:4257448/150‑1 (MQ=255)
                                                                                                                      tCCTGGTCGAGGCGGCGCTGAGGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCa                             <  1:1597050/150‑1 (MQ=255)
                                                                                                                            tCGAGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgc                       <  1:9481978/150‑2 (MQ=255)
                                                                                                                               aGGCGGCGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAggccgccgt                    >  1:4043282/1‑146 (MQ=255)
                                                                                                                                ggcggcGCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGtt                  >  1:1849940/1‑150 (MQ=255)
                                                                                                                                  cggcgCTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTtgt                <  1:9239928/150‑1 (MQ=255)
                                                                                                                                    gcgcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTgg              >  2:3947768/1‑150 (MQ=255)
                                                                                                                                    gcgcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAAGTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTgg              >  1:8731311/1‑150 (MQ=255)
                                                                                                                                      gcTGATGGCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGAt            >  2:5194639/1‑150 (MQ=255)
                                                                                                                                            ggCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCgg      >  1:10630325/1‑150 (MQ=255)
                                                                                                                                            ggCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGTGTGGATTCGCgg      >  1:369444/1‑150 (MQ=255)
                                                                                                                                             gCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGc     <  1:3888519/150‑1 (MQ=255)
                                                                                                                                             gCTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGc     <  1:6532691/150‑1 (MQ=255)
                                                                                                                                              cTGCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGTTGTGGATTCGCGGCt    <  1:527273/150‑1 (MQ=255)
                                                                                                                                                gCCGATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCAAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGCCGGGGTGGATTCGCGGCTAc  >  1:8515374/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGGGCGCCGCGCCGCTGAAGGAGAACCTGGCGGCGGCGGTGTTGATCCGTGCCGGCTGGCCGAAGATCGCCGCCGAGGGCGGTGCCCTGGCCGACCCGATGTGTGGCGTCGGCACCTTCCTGGTCGAGGCGGCGCTGATGGCTGCCCATATCGCGCCGAACCTCAGGCGCGAGCGCTGGGGCTTCAGCAATTGGCTGGGGCACGTGCCGGCGCTGTGGCGCCAGCTCCACGAGGAAGCCCAGCAGCGCGCGGCGGCCGGCCTGGCCAGGCCGGCCGTTGTGGATTCGCGGCTAC  >  NC_008463/2155924‑2156218

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: