Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,187 G→A 6.0% R2H (CGT→CAT)  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,1870GA 6.0% 236.9 / 4.6 85R2H (CGT→CAT) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (37/42);  new base A (2/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

GTACTGCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGG  >  NC_002947/2389038‑2389336
                                                                                                                                                     |                                                                                                                                                     
gTACTGCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCg                                                                                                                                                       <  2:1887790/150‑1 (MQ=255)
     gCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGcc                                                                                                                                                  >  1:2418176/1‑150 (MQ=255)
       cTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACCGGCAAACATCGCGGGGCGACGCATCGcccc                                                                                                                                                <  2:1409663/150‑1 (MQ=255)
        tgatgaTCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCg                                                                                                                                               >  2:2413290/1‑150 (MQ=255)
         gatgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                              >  2:2265322/1‑150 (MQ=255)
           tgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTCCGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtca                                                                                                                                            <  1:2016167/150‑1 (MQ=255)
            gatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcag                                                                                                                                           <  2:1698661/150‑1 (MQ=255)
                tAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGTAGCACGTTGTGTGCCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaa                                                                                                                                       >  2:1538887/1‑150 (MQ=255)
                     tGTCCTTGCCACTGGGAAAAGCACCATTCGGTTCGGTGGTCCGTCTTTTCGGCAGCACGTTTTGTGTCGCAGAAGGCAACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGAGGCCTCGCCCCGTCAGTCAAATAGc                                                                                                                                  >  2:1207372/1‑150 (MQ=255)
                      gTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGcc                                                                                                                                 <  1:887711/150‑1 (MQ=255)
                       tCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                                >  2:8437/1‑150 (MQ=255)
                       tCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGGTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCTGGGCGATGCGTTGCCCCGTCAGTCAAATAGccc                                                                                                                                >  2:1219526/1‑150 (MQ=255)
                         cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCaa                                                                                                                              <  2:16553/150‑1 (MQ=255)
                         cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCATCGCCCCGTCAGTCAAATAGCCCaa                                                                                                                              >  2:2028876/1‑150 (MQ=255)
                          ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAg                                                                                                                             >  1:713808/1‑150 (MQ=255)
                              cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGAGTTGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcgc                                                                                                                         >  1:21862/1‑150 (MQ=255)
                               aCTGGGCCAGGCACGGATCGGGCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                                        >  1:798587/1‑150 (MQ=255)
                                cTGGGCCAGGCGCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGGGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                                       <  2:149542/150‑1 (MQ=255)
                                 tGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAATGTGTACGGGCAAGCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTgg                                                                                                                      <  2:1326000/150‑1 (MQ=255)
                                     ccAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACAACGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATa                                                                                                                  >  2:1455613/1‑150 (MQ=255)
                                           aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCCCCCGCACGCCTGCTGGACAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATCCTCATg                                                                                                            >  2:803858/1‑150 (MQ=255)
                                           aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATg                                                                                                            >  1:1050953/1‑150 (MQ=255)
                                           aCGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATg                                                                                                            >  1:1131909/1‑150 (MQ=255)
                                            cGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCACGt                                                                                                           <  1:346124/150‑1 (MQ=255)
                                             ggATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGCCAAATAGCCCAAGGCGCTGGAATACTCATGTc                                                                                                          <  2:1345550/150‑1 (MQ=255)
                                                tctgTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCg                                                                                                                                                    >  2:1270163/4‑105 (MQ=255)
                                                tctgTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCg                                                                                                                                                    <  1:1270163/102‑1 (MQ=255)
                                                 cGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAc                                                                                                      <  2:1634531/150‑1 (MQ=255)
                                                      aGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGttt                                                                                                 <  2:2418176/150‑1 (MQ=255)
                                                      aGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGttt                                                                                                 >  2:88530/1‑150 (MQ=255)
                                                           acccGACTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAAc                                                                                            >  2:1588637/2‑150 (MQ=255)
                                                             ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                          >  1:215349/1‑150 (MQ=255)
                                                             ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                          >  1:2355286/1‑150 (MQ=255)
                                                             ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                          >  2:920383/1‑150 (MQ=255)
                                                              cGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACttt                                                                                         >  1:1129607/1‑150 (MQ=255)
                                                                gCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTAc                                                                                       <  1:2028876/150‑1 (MQ=255)
                                                                 cTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCACGTCGGACCGTTTCGAACTTTAcc                                                                                      <  1:2265322/150‑1 (MQ=255)
                                                                  ttttCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCt                                                                                     >  2:1626710/1‑150 (MQ=255)
                                                                    ttCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCa                                                                                   <  1:2413290/150‑1 (MQ=255)
                                                                           gcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                            <  2:713808/150‑1 (MQ=255)
                                                                             aCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGaa                                                                                             >  1:622052/1‑131 (MQ=255)
                                                                             aCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGaa                                                                                             <  2:622052/131‑1 (MQ=255)
                                                                                  gtgtgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGc                                                                     >  1:2014589/1‑150 (MQ=255)
                                                                                  gtgtgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGc                                                                     <  1:1295139/150‑1 (MQ=255)
                                                                                    gggtCGCAGAAGGCCAACCGCACGCATGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCt                                                                   <  2:21862/148‑1 (MQ=255)
                                                                                     tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                  <  1:1588637/150‑1 (MQ=255)
                                                                                     tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                  >  2:2052536/1‑150 (MQ=255)
                                                                                      gtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTg                                                                 >  1:217987/1‑150 (MQ=255)
                                                                                         gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGACTACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                              <  1:1974287/150‑1 (MQ=255)
                                                                                         gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                              <  1:1538887/150‑1 (MQ=255)
                                                                                         gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                              <  2:1050953/150‑1 (MQ=255)
                                                                                         gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                              <  2:2124592/150‑1 (MQ=255)
                                                                                         gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGGTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                              >  1:1706451/1‑150 (MQ=255)
                                                                                             aaGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCt                                                          <  1:2052536/150‑1 (MQ=255)
                                                                                                 ccACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTc                                                      >  1:982260/1‑150 (MQ=255)
                                                                                                  cACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTc                                                                             >  2:2086660/1‑126 (MQ=255)
                                                                                                  cACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTc                                                                             <  1:2086660/126‑1 (MQ=255)
                                                                                                   aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGt                                                                                                           <  1:815632/95‑1 (MQ=255)
                                                                                                   aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGt                                                                                                           >  2:2044872/1‑95 (MQ=255)
                                                                                                   aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGt                                                                                                           <  1:2044872/95‑1 (MQ=255)
                                                                                                   aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGt                                                                                                           >  2:815632/1‑95 (MQ=255)
                                                                                                   aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGc                                                    >  2:846093/1‑150 (MQ=255)
                                                                                                    cccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGcc                                                   <  1:870767/150‑1 (MQ=255)
                                                                                                     ccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGa                                                                                              >  1:174183/1‑106 (MQ=255)
                                                                                                     ccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGa                                                                                              <  2:174183/106‑1 (MQ=255)
                                                                                                       gCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGaa                                                <  2:2014589/150‑1 (MQ=255)
                                                                                                         aCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                              <  1:740437/150‑1 (MQ=255)
                                                                                                           gCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAgg                                            <  2:1837577/150‑1 (MQ=255)
                                                                                                                 tgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                            >  1:2087884/1‑112 (MQ=255)
                                                                                                                 tgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTcc                                                                            <  2:2087884/112‑1 (MQ=255)
                                                                                                                    aaaGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGcc                                   >  2:1999135/1‑150 (MQ=255)
                                                                                                                         tgtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGc                              <  1:846093/150‑1 (MQ=255)
                                                                                                                               ggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGAc                        <  2:108507/150‑1 (MQ=255)
                                                                                                                                gCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACg                       <  2:493790/150‑1 (MQ=255)
                                                                                                                                 cAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGa                      <  1:8437/150‑1 (MQ=255)
                                                                                                                                  aaaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAg                     <  2:2339032/150‑1 (MQ=255)
                                                                                                                                  aaaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAg                     >  1:2396571/1‑150 (MQ=255)
                                                                                                                                  aaaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAg                     <  2:798587/150‑1 (MQ=255)
                                                                                                                                         acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                              >  1:251310/2‑118 (MQ=255)
                                                                                                                                         acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                              <  2:251310/117‑1 (MQ=255)
                                                                                                                                          cgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa             <  2:1129607/150‑1 (MQ=255)
                                                                                                                                           ggggCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAg            <  1:1455613/150‑1 (MQ=255)
                                                                                                                                             ggCGTTGCGACGGCACGCCAGTCAGATAGCCCAAGGCGCTGGAATACTCATGTCGCACCGTTTGGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCAAGca          <  2:1131909/150‑1 (MQ=255)
                                                                                                                                                 aTGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGaa                                                >  2:1334997/1‑108 (MQ=255)
                                                                                                                                                 aTGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGaa                                                <  1:1334997/108‑1 (MQ=255)
                                                                                                                                                 aTGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGCGCGAGCACAcc      >  2:1617791/1‑150 (MQ=255)
                                                                                                                                                     gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCgg  <  1:920383/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTACTGCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGG  >  NC_002947/2389038‑2389336

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: ‑