Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,948 G→T 6.3% A256S (GCG→TCG)  rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,389,9480GT 6.3% 201.8 / 3.0 64A256S (GCG→TCG) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (34/26);  new base T (2/2);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

TGCTGCCGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGCCGGCCTGGGC  >  NC_002947/2389806‑2390091
                                                                                                                                              |                                                                                                                                               
tgctgcCGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCtgtg                                                                                                                                          <  2:505516/150‑1 (MQ=255)
 gctgcCGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTgg                                                                                                                                         >  1:345218/1‑150 (MQ=255)
  ctgcCGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGgc                                                                                                                                        >  2:1765840/1‑150 (MQ=255)
     ccGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGATTGGCTCGGCCAGGTGCCGGCGCTGTGGCGCa                                                                                                                                     >  2:452318/1‑150 (MQ=255)
       gAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAg                                                                                                                                   <  2:345218/150‑1 (MQ=255)
          ggtggtGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCtgt                                                                                                                                           <  1:89180/139‑1 (MQ=255)
          ggtggtGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCtgt                                                                                                                                           >  2:89180/1‑139 (MQ=255)
             ggtgCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCAt                                                                                                                             <  2:1457366/150‑1 (MQ=255)
             ggtgCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCAt                                                                                                                             <  2:799982/150‑1 (MQ=255)
              gtgCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATg                                                                                                                            >  1:1747756/1‑150 (MQ=255)
                    tGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCGGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTAGCTCTGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGGCGAAg                                                                                                                      >  2:108342/1‑150 (MQ=255)
                    tGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGaa                                                                                                                       <  2:949596/149‑1 (MQ=255)
                    tGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGaa                                                                                                                       >  1:949596/1‑149 (MQ=255)
                     ggCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAgc                                                                                                                     <  2:1886115/150‑1 (MQ=255)
                         gACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAg                                                                                                                 <  2:736135/150‑1 (MQ=255)
                            ccGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCtgtg                                                                                                                                          >  1:1448877/1‑122 (MQ=255)
                            ccGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCtgtg                                                                                                                                          <  2:1448877/122‑1 (MQ=255)
                            ccGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGgc                                                                                                                                        >  1:430714/1‑124 (MQ=255)
                            ccGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGgc                                                                                                                                        <  2:430714/124‑1 (MQ=255)
                             cGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgc                                                                                                             >  2:1893433/1‑150 (MQ=255)
                               aTGTGCGGCGTGGGTACCTTCTTGGTCGAGGCCGCGATGATCGCCGCCGATATCGGGCTCAACCTAAAGCGCGAACGCTGGGGCTTCAGTGCTTGTCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCCGgcgcgc                                                                                                           >  2:163250/1‑150 (MQ=255)
                               aTGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACCAGGCTCAGgcgcgc                                                                                                           >  1:400464/1‑150 (MQ=255)
                                tgtgCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgcgcg                                                                                                          >  2:1550077/1‑150 (MQ=255)
                                         tGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCCAAGCGAAGCGCAGGCGCGCGCGCAggccg                                                                                                 <  1:1507157/150‑1 (MQ=255)
                                          gggTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAggccgg                                                                                                >  1:1696156/1‑150 (MQ=255)
                                           ggTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAggccggc                                                                                               <  1:1628864/150‑1 (MQ=255)
                                             tACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCt                                                                                             <  2:1696156/150‑1 (MQ=255)
                                                       gTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCa                                                                                   >  2:1254967/1‑150 (MQ=255)
                                                          gAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGgc                                                                                                                                        >  1:1291785/1‑94 (MQ=255)
                                                          gAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGgc                                                                                                                                        <  2:1291785/94‑1 (MQ=255)
                                                             gCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGTCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgc                                                                                                             >  1:364489/1‑118 (MQ=255)
                                                             gCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGTCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgc                                                                                                             >  1:364406/1‑118 (MQ=255)
                                                             gCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGTCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgc                                                                                                             <  2:364489/118‑1 (MQ=255)
                                                             gCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGTCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGgcgc                                                                                                             <  2:364406/118‑1 (MQ=255)
                                                             gCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCtg                                                                             <  2:1503719/150‑1 (MQ=255)
                                                              ccGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCtgt                                                                            <  2:589063/150‑1 (MQ=255)
                                                               cgcgATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCtgtg                                                                           <  1:1550077/150‑1 (MQ=255)
                                                                     gatCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGCCTTCAGTGCTTGGCTCGGCCATGTGCCCGCGCTGTGGCGCAAGGTGCATGACGAAGCGCCCGCGCCCGCGCAGCCCGGCCTCGCCAACCCACCGCTGTGGATCCg                                                                     >  2:402549/1‑150 (MQ=255)
                                                                      ntCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTNAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGt                                                                    >  1:306747/2‑150 (MQ=255)
                                                                        cgccgcCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTgg                                                                  <  1:1254967/150‑1 (MQ=255)
                                                                                  atCGCCCCCAACCTCAAGCGTGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGC‑‑TGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGa                                                      <  1:563381/150‑1 (MQ=255)
                                                                                    cGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGGGCAGGCCTGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGa                                                      <  1:824178/150‑1 (MQ=255)
                                                                                    cGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGa                                                      <  1:1738107/150‑1 (MQ=255)
                                                                                     gCCCCCAACCTCCAGCTCTAACGCTGGGGCTTCAGTGCTTGGCTCGCCCAGGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAACGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGAc                                                     <  2:306747/150‑1 (MQ=255)
                                                                                         ccAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCAc                                                 >  1:633216/1‑150 (MQ=255)
                                                                                          cAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGACCGGCCTGGCGAAGCCACCGCCGTGGAGGCGTGGCTATGAGGCTggcacccg                                                >  2:907693/1‑143 (MQ=255)
                                                                                              cTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTg                                            >  1:1502836/1‑150 (MQ=255)
                                                                                              cTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTAGCCAAGCAACCGCTGTGGATCCTTGGCTATGAGGCTGACCCACGGCTg                                            >  1:1883875/1‑150 (MQ=255)
                                                                                               tCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGa                                           >  1:1683462/1‑150 (MQ=255)
                                                                                                            cTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGc                                     <  1:1666111/143‑1 (MQ=255)
                                                                                                            cTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGc                                     >  2:1666111/1‑143 (MQ=255)
                                                                                                             tGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGc                             >  2:91562/1‑150 (MQ=255)
                                                                                                               gggCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCTGAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGGCGcaa                           >  1:1590468/1‑150 (MQ=255)
                                                                                                               gggCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGcaa                           <  1:1765840/150‑1 (MQ=255)
                                                                                                                   ttCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCaacaac                       >  2:717824/1‑150 (MQ=255)
                                                                                                                   ttCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCaacaac                       >  1:1489256/1‑150 (MQ=255)
                                                                                                                        tGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCACGACGAAGCGCAGGCGCGCGCGCAg                                                                                                     >  1:33184/1‑67 (MQ=255)
                                                                                                                        tGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCACGACGAAGCGCAGGCGCGCGCGCAg                                                                                                     <  2:33184/67‑1 (MQ=255)
                                                                                                                              gCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGc            >  2:877511/1‑150 (MQ=255)
                                                                                                                                   gCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACAACAATGGATACGAGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGCCGGcc       >  1:1789986/1‑150 (MQ=255)
                                                                                                                                        gtgCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGCCGGCCTGGGc  >  2:246291/1‑150 (MQ=255)
                                                                                                                                        gtgCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGCCGGCCTGGGc  >  2:1010682/1‑150 (MQ=255)
                                                                                                                                           ccGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCt                                                                              >  1:1376876/1‑71 (MQ=255)
                                                                                                                                           ccGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCt                                                                              <  2:1376876/71‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
TGCTGCCGAAGGTGGTGCGCTGGCCGACCCGATGTGCGGCGTGGGTACCTTCCTGGTCGAGGCCGCGATGATCGCCGCCGATATCGCCCCCAACCTCAAGCGCGAACGCTGGGGCTTCAGTGCTTGGCTCGGCCATGTGCCGGCGCTGTGGCGCAAGGTGCATGACGAAGCGCAGGCGCGCGCGCAGGCCGGCCTGGCCAAGCCACCGCTGTGGATCCGTGGCTATGAGGCTGACCCACGGCTGATCCAGCCGGGCCGCAACAACGTCGAGCGTGCCGGCCTGGGC  >  NC_002947/2389806‑2390091

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