Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,391,109 G→A 5.4% 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,391,1090GA 5.4% 634.6 / 15.9 201R2H (CGT→CAT) rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (79/111);  new base A (5/6);  total (84/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CTGCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTC  >  NC_002947/2390963‑2391256
                                                                                                                                                  |                                                                                                                                                   
    cTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGcccc                                                                                                                                              >  2:75520/1‑150 (MQ=255)
    cTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGACGCATCGcccc                                                                                                                                              <  1:1226268/150‑1 (MQ=255)
     tgatgaTCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCg                                                                                                                                             >  1:2927205/1‑150 (MQ=255)
      gatgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgt                                                                                                                                            >  2:3426021/1‑150 (MQ=255)
       atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtc                                                                                                                                           <  2:943830/150‑1 (MQ=255)
       atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtc                                                                                                                                           <  2:3850450/150‑1 (MQ=255)
       atgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtc                                                                                                                                           >  2:4259402/1‑150 (MQ=255)
        tgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtca                                                                                                                                          <  1:2317659/150‑1 (MQ=255)
        tgatCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtca                                                                                                                                          <  1:2346919/150‑1 (MQ=255)
          atCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagt                                                                                                                                        <  2:785190/150‑1 (MQ=255)
           tCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCgtcagtc                                                                                                                                       >  2:1489673/1‑150 (MQ=255)
            cTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGcccc                                                                                                                                              >  2:2492002/1‑142 (MQ=255)
             tAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCaa                                                                                                                                     >  2:2186123/1‑150 (MQ=255)
                gCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATa                                                                                                                                  <  1:5103069/150‑1 (MQ=255)
                  tGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGc                                                                                                                                <  2:243868/150‑1 (MQ=255)
                    tCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGccc                                                                                                                              <  1:5054540/150‑1 (MQ=255)
                      cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCaa                                                                                                                            >  2:4516678/1‑150 (MQ=255)
                      cTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCaa                                                                                                                            >  2:570476/1‑150 (MQ=255)
                       ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAg                                                                                                                           <  2:693975/150‑1 (MQ=255)
                       ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGACGCATCGcccc                                                                                                                                              >  2:4189388/1‑131 (MQ=255)
                       ttGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGACGCATCGcccc                                                                                                                                              <  1:4189388/131‑1 (MQ=255)
                        tGCCACTGGGCCAGGCACGGATCGGACAGGCGGCCCGGCTTTTCTGCAGCACGTAGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAAAATCACAAgg                                                                                                                          >  1:3766847/1‑150 (MQ=255)
                         gCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgc                                                                                                                         <  1:4609173/150‑1 (MQ=255)
                          ccACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcg                                                                                                                        >  2:1696250/1‑150 (MQ=255)
                           cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcgc                                                                                                                       >  1:4446010/1‑150 (MQ=255)
                           cACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGgcgc                                                                                                                       >  2:190981/1‑150 (MQ=255)
                            aCTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                                      <  1:1014110/150‑1 (MQ=255)
                            aCTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCt                                                                                                                      >  2:4275170/1‑150 (MQ=255)
                             cTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCCGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGGTGGAAAGGTGGACGGGCAAACATCGCGGGGCGATGCGGCGCCCCGTCAGTCCAATAGCCCCGGGCGCTg                                                                                                                     >  2:4593048/1‑150 (MQ=255)
                                 gCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAAt                                                                                                                 >  1:2250263/1‑150 (MQ=255)
                                  ccAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATa                                                                                                                <  1:4636808/150‑1 (MQ=255)
                                   cAGGCACGGCTCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATAc                                                                                                               <  1:3894082/150‑1 (MQ=255)
                                   cAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATAc                                                                                                               >  2:1888430/1‑150 (MQ=255)
                                    aGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGTGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACt                                                                                                              <  1:2479333/150‑1 (MQ=255)
                                     ggCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTc                                                                                                             <  2:2399617/150‑1 (MQ=255)
                                           gATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                                     <  1:4758915/136‑1 (MQ=255)
                                             tCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGa                                                                                                     <  2:513805/150‑1 (MQ=255)
                                              cGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGAc                                                                                                    <  1:76843/150‑1 (MQ=255)
                                                gTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCg                                                                                                  >  1:3418774/1‑150 (MQ=255)
                                                    ggcggcCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTc                                                                                              <  1:1157111/150‑1 (MQ=255)
                                                     gcggcCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCg                                                                                             >  1:1372324/1‑150 (MQ=255)
                                                          ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                        >  1:4694114/1‑150 (MQ=255)
                                                          ccGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                        >  2:4825278/1‑150 (MQ=255)
                                                             gCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCTCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTAc                                                                                     >  1:2054431/1‑150 (MQ=255)
                                                             gCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTAc                                                                                     <  2:787240/150‑1 (MQ=255)
                                                             gCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTAc                                                                                     <  2:3410006/150‑1 (MQ=255)
                                                               ttttCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCt                                                                                   >  2:2045396/1‑150 (MQ=255)
                                                               ttttCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCt                                                                                   >  1:629199/1‑150 (MQ=255)
                                                                   cTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCAcc                                                                               <  1:41216/150‑1 (MQ=255)
                                                                    tGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCt                                                                              <  2:3997016/150‑1 (MQ=255)
                                                                    tGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCt                                                                              >  1:3497748/1‑150 (MQ=255)
                                                                      cagcaCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGt                                                                            >  2:3927363/1‑150 (MQ=255)
                                                                      cagcaCGTTGTGTGTCGCAGAAGGCAACACGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGAGTCGCACCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAAGTCGGACCGTTTCGAAATTTACCTCTCCTgg                                                                            >  2:3644618/1‑149 (MQ=255)
                                                                         caCGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTccc                                                                         >  2:3266055/1‑150 (MQ=255)
                                                                             ttGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAg                                                                     <  1:549979/150‑1 (MQ=255)
                                                                                  tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                <  1:956644/150‑1 (MQ=255)
                                                                                  tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                <  2:2315931/150‑1 (MQ=255)
                                                                                  tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCtt                                                                <  1:3395690/150‑1 (MQ=255)
                                                                                  tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCACGGGGCGACGCATCGCCCCgc                                                                                                                                            <  1:2634045/74‑2 (MQ=38)
                                                                                  tgtCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCACGGGGCGACGCATCGCCCCgc                                                                                                                                            >  2:2634045/1‑73 (MQ=38)
                                                                                      gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                            <  2:2250263/150‑1 (MQ=255)
                                                                                      gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                            <  1:2144274/150‑1 (MQ=255)
                                                                                      gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                            >  1:711662/1‑150 (MQ=255)
                                                                                      gCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGaaa                                                            <  1:578309/150‑1 (MQ=255)
                                                                                        aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                          <  1:2994206/150‑1 (MQ=255)
                                                                                        aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                          <  1:899412/150‑1 (MQ=255)
                                                                                        aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                          <  2:3416161/150‑1 (MQ=255)
                                                                                        aGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGc                                                          <  1:4843987/150‑1 (MQ=255)
                                                                                               cACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCg                                                   >  2:4708673/1‑150 (MQ=255)
                                                                                               cACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCAGGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCg                                                   >  1:4359911/1‑150 (MQ=255)
                                                                                                aCCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGACGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGc                                                  <  2:3785105/150‑1 (MQ=255)
                                                                                                  ccGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCg                                                >  2:1648735/1‑150 (MQ=255)
                                                                                                    gCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGaa                                              <  2:726302/150‑1 (MQ=255)
                                                                                                      aCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                            >  1:4928444/1‑150 (MQ=255)
                                                                                                      aCGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                            <  2:2114983/150‑1 (MQ=255)
                                                                                                       cGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                           <  1:4966073/150‑1 (MQ=255)
                                                                                                       cGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAagag                                           <  1:1746999/150‑1 (MQ=255)
                                                                                                        gCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAgg                                          <  1:4338862/150‑1 (MQ=255)
                                                                                                           tgctgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACtt                                                                                        <  1:3029883/101‑1 (MQ=255)
                                                                                                            gctgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCaa                                      <  1:2789671/150‑1 (MQ=255)
                                                                                                              tgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                    >  2:1072767/1‑150 (MQ=255)
                                                                                                              tgGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAgg                                    >  1:2199043/1‑150 (MQ=255)
                                                                                                                gAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGc                                  >  2:2538466/1‑150 (MQ=255)
                                                                                                                  aaGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCACGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCt                                <  2:2696434/150‑1 (MQ=255)
                                                                                                                       gtACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCCTTGGCCTcttcg                           <  2:1309146/150‑5 (MQ=255)
                                                                                                                           gggCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGa                       <  1:1579916/150‑1 (MQ=255)
                                                                                                                              cAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGa                    <  1:132730/150‑1 (MQ=255)
                                                                                                                              cAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGa                    <  1:3996435/150‑1 (MQ=255)
                                                                                                                               aaaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAg                   >  1:3645850/1‑150 (MQ=255)
                                                                                                                                aaCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAgg                  >  1:1925228/1‑150 (MQ=255)
                                                                                                                                 aCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTg                                                                                                                     <  1:3853128/50‑1 (MQ=255)
                                                                                                                                 aCATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTAACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGt                 <  1:3305449/150‑1 (MQ=255)
                                                                                                                                    tgacgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTgcg              >  1:4772044/4‑150 (MQ=255)
                                                                                                                                      gcgGGGCGATGCGTCGCCGCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGcgcg            <  2:1753908/150‑1 (MQ=255)
                                                                                                                                      acgGGGCGACGCATCGCCCCGTCAGTCAAATAGCCTAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAaga                                            >  2:4767715/2‑118 (MQ=255)
                                                                                                                                       cgGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa           <  2:1514062/150‑1 (MQ=255)
                                                                                                                                        ggggCGACGCATCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAg          >  2:2884324/1‑150 (MQ=255)
                                                                                                                                         gggCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGc         <  2:3062505/150‑1 (MQ=255)
                                                                                                                                             gATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCacac     >  2:4156293/1‑150 (MQ=255)
                                                                                                                                                gCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGTGCTCGCGCACCTc  <  1:437419/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGCCTGATGATCTAAGCTGTCCTTGCCACTGGGCCAGGCACGGATCGGTCAGGCGGCCCGGCTTTTCTGCAGCACGTTGTGTGTCGCAGAAGGCCACCCGCACGCCTGCTGGAAAGGTGTACGGGCAAACATCGCGGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTC  >  NC_002947/2390963‑2391256

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: