Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,276,006 C→A 6.5% intergenic (+58/+48) rlmG → / ← PP_4649 ribosomal RNA large subunit methyltransferase G/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,276,0060CA 6.5% 502.7 / 2.5 159intergenic (+58/+48)rlmG/PP_4649ribosomal RNA large subunit methyltransferase G/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (83/62);  new base A (2/12);  total (85/74)
Fisher's exact test for biased strand distribution p-value = 3.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 162 total aligned reads displayed.

TGGGCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAA  >  NC_002947/5275857‑5276153
                                                                                                                                                      |                                                                                                                                                    
tGGGCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGcc                                                                                                                                                      >  2:3089902/1‑150 (MQ=255)
tGGGCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGcc                                                                                                                                                      <  2:1834748/150‑1 (MQ=255)
tGGGCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGcc                                                                                                                                                      >  1:811465/1‑150 (MQ=255)
 gggCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCa                                                                                                                                                     >  1:2667456/1‑150 (MQ=255)
 gggCTATCACAGCAAGCTGCCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCAATCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTCTGCCa                                                                                                                                                     >  2:5001517/1‑150 (MQ=255)
   gCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATc                                                                                                                                                   <  1:4836977/150‑1 (MQ=255)
     tATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTa                                                                                                                                                 <  1:43060/150‑1 (MQ=255)
       tCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgc                                                                                                                                               >  1:1016218/1‑150 (MQ=255)
       tCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgc                                                                                                                                               >  1:1015867/1‑150 (MQ=255)
         acaGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgcgc                                                                                                                                             <  2:958168/150‑1 (MQ=255)
         acaGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAgcgc                                                                                                                                             <  2:4503294/150‑1 (MQ=255)
            gCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGtt                                                                                                                                          <  2:3163296/150‑1 (MQ=255)
             cAAGCTGGCGCGGTTGTTCAGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTg                                                                                                                                         >  2:3687317/1‑150 (MQ=255)
                gCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCTCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCt                                                                                                                                      <  2:2637137/150‑1 (MQ=255)
                gCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCt                                                                                                                                      >  2:1944038/1‑150 (MQ=255)
                gCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCt                                                                                                                                      >  2:3142959/1‑150 (MQ=255)
                gCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCt                                                                                                                                      >  2:2778761/1‑150 (MQ=255)
                 cTGGCGCGGGTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTg                                                                                                                                     >  2:2278774/1‑150 (MQ=255)
                  tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGccggt                                                                                                                                    <  2:3433154/150‑4 (MQ=255)
                  tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTgt                                                                                                                                    <  1:757130/150‑2 (MQ=255)
                  tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGc                                                                                                                                    <  1:4476585/150‑1 (MQ=255)
                  tGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGc                                                                                                                                    <  1:4439966/150‑1 (MQ=255)
                   ggCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCt                                                                                                                                   >  2:2869540/1‑150 (MQ=255)
                      gcgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCCCGAACACCGGCAACGCCTGTGCCATCTAGCCCGTTGCCTGCTTCg                                                                                                                                >  1:2732396/1‑150 (MQ=255)
                       cgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               <  1:5234293/150‑1 (MQ=255)
                       cgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               <  1:882763/150‑1 (MQ=255)
                       cgGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               <  1:3448785/150‑1 (MQ=255)
                        ggTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTg                                                                                                                              <  2:1768863/150‑1 (MQ=255)
                         tgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             >  2:3909317‑M1/89‑150 (MQ=255)
                         tgagcgtcgcagcgctgggtgtggcggcctgtgtgatgcggctgttcatgggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             >  1:2874642‑M1/89‑150 (MQ=255)
                         gttgttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             <  1:144789/150‑1 (MQ=255)
                         gttgttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             >  1:197112/1‑150 (MQ=255)
                          ttgttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTggg                                                                                                                            >  1:2964139/1‑150 (MQ=255)
                             ttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCaa                                                                                                                         <  1:1623788/150‑1 (MQ=255)
                             ttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCaa                                                                                                                         <  2:250861/150‑1 (MQ=255)
                             ttCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCaa                                                                                                                         <  1:3065252/150‑1 (MQ=255)
                                cGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGcc                                                                                                                     >  2:3576853/1‑150 (MQ=255)
                                 ggggggTTGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGccc                                                                                                                    <  1:3597510/150‑1 (MQ=255)
                                 ggggggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGccc                                                                                                                    >  2:4203608/1‑150 (MQ=255)
                                 ggggggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGccc                                                                                                                    <  1:4899776/150‑1 (MQ=255)
                                   ggggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGc                                                                                                                  >  1:902850/1‑150 (MQ=255)
                                     ggTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCtc                                                                                                                <  1:2782859/150‑2 (MQ=255)
                                           gCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGctc                                                                                                          >  2:502812/1‑150 (MQ=255)
                                            cAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGtt                                                                                                                                          <  2:2536380/118‑1 (MQ=255)
                                               gTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCAGTGCCAACTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGcgcgtag                                                                                                      >  1:4102526/1‑144 (MQ=255)
                                                  gcggcAACGCCCAAGTTCGTCTTTCTGAATCCCCGAATGTTGGGGCGCCACCGCCCCCTTTGGTGCGGGCTTGCCCGCTAAAACCGGCAACGCGTTTGCAATCTCGGGGGTTGCCTGCTTGGTGGGCACGCCCGCTGCCGCTCTTTATCa                                                                                                   <  2:1686884/150‑1 (MQ=255)
                                                  gcggcAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCa                                                                                                   >  1:1740027/1‑150 (MQ=255)
                                                  gcggaAACGCCCACTTTCGTCTTTCTGAACGCCCGCAAGTCGGCGCGCCACCGCCCCCTTCGGAGGGGGCTGGCCGGCTAACACGGGCAACGCCTGTGCAATCTCGCGCGTTGCCTGATTGGTGGGAAAGCCCGCTGCCGCTCTTTATCa                                                                                                   <  2:1735957/150‑1 (MQ=255)
                                                           cccattttCGTAATTATTAAAGCCCGCAAGTCGGCGCGCCCCCCCCCCCGTCGGCGGGGGGTTGCCCGCTAACCCCGGAAATGCCTTTGCCATCTAGCGGTTTGCCTGATTCGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGt                                                                                          <  2:4699457/144‑1 (MQ=255)
                                                           cccAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCAGCGAACACCGGCAACGCCTGTGCCATCTTGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGt                                                                                          <  2:4535414/150‑1 (MQ=255)
                                                              aaGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAg                                                                                       >  2:548365/1‑150 (MQ=255)
                                                              aaGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAg                                                                                       >  2:3142107/1‑150 (MQ=255)
                                                                    gTCATTCTGAAAGCCCGCAAGTTGGCGCGCCACCGCCCCCGTTGGCGCGGGCTGGCCCGCCAACACCGGAAAGGCGTGTGCAATCTCGCGGGTTGCCTGCTTCGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCccgc                                                                                 <  2:2835569/150‑1 (MQ=255)
                                                                         tCTGAACGCCCGAATTTCGGCGCGCCACCGCCCCGGTTGGAGGGGGCTTGCCCGCCAACACCGCAAACGCCTTTGCAATCTCGCGCTTGGCCTGCTTCGTGGGCAACCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccc                                                                            <  2:3833830/150‑1 (MQ=255)
                                                                         tCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccc                                                                            >  2:1483912/1‑150 (MQ=255)
                                                                           tGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccccc                                                                          >  2:3930962/1‑150 (MQ=255)
                                                                           tGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccccc                                                                          >  1:349641/1‑150 (MQ=255)
                                                                           tGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGccccc                                                                          >  2:2135200/1‑150 (MQ=255)
                                                                           ggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             >  2:2329722‑M1/39‑100 (MQ=255)
                                                                           ggggcggcgggcctgaccacacactgataaagagcggcAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTgg                                                                                                                             <  1:2329722‑M1/62‑1 (MQ=255)
                                                                            gAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               <  1:3093055/97‑1 (MQ=255)
                                                                            gAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               <  1:3093435/97‑1 (MQ=255)
                                                                            gAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGt                                                                                                                               >  2:3093435/1‑97 (MQ=255)
                                                                                   cGCATGTAGGCGCGCCACCGCCCCCTTGGGCGGGGGTTGGCCCGCGACCACGGGCAACGCCTGTGCATTCTCGCGGGTTGCCTGCTTCGTGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAg                                                                  <  2:5179821/150‑1 (MQ=255)
                                                                                   cGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGGCTGCTTCGTGGGCTAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAagag                                                                  >  2:3790507/1‑150 (MQ=255)
                                                                                    gCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGc                                                                 >  1:3872596/1‑150 (MQ=255)
                                                                                      aaGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCg                                                               >  1:3892293/1‑150 (MQ=255)
                                                                                            gcgcgcCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcac                                                         >  1:3349403/1‑150 (MQ=255)
                                                                                             cgcgcCACCGCCCCGGTTGGTGGGGGGGTGCCCGCTAACACGGGAAACGCGTTTGCATTCTTGCGGGTTGCCTGCTGCGGGGGCACGCCGGCTGCCGCTCTTTATCAGTGTGTTGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcaca                                                        <  2:3256312/150‑1 (MQ=255)
                                                                                              gcgcCCCCGCCCCGGTTGGTGGGGGTTTGCCCGCTAACACCGGCAACGCCTGTGCCATCTCGGGGTTGGCCTGCTGGGTGGCCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcacac                                                       <  2:1549235/150‑1 (MQ=255)
                                                                                              gcgcCACCGCCCGGTTTGGAGGGTGTTTGCCGGGAACCACCGCAAAGGCCTTTGCAATCTGGCGTTTGGCTTTTTTGTTGGTAATGCCCGCTGCCGCTCTTTATCAGTGTGTGTTCAGGCCCGCCGCCCCCATGAACAGCCGCATcacac                                                       <  2:1549179/150‑1 (MQ=255)
                                                                                               cccccccGCCCCGGTTGGGGCGGGCTTGCCCGCAACCACCGCCCCCGCGTGGGCCATCTCGCGGTTGGCGGGGTTGTTGGCAACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATcacaca                                                      <  1:4764543/145‑1 (MQ=255)
                                                                                                 ccACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAgg                                                    >  1:2200071/1‑150 (MQ=255)
                                                                                                  cACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGgc                                                   >  1:4034853/1‑150 (MQ=255)
                                                                                                     cGCCCCGGTTGGCGGGGGCTGCCCCGCTACACCCGCAACCGCTTGTGCAATCTCGGGGGTTGCGGGTTTGGGGGGAACCCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGgccgc                                                <  1:3224310/150‑1 (MQ=255)
                                                                                                        ccccGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCcac                                             >  1:4186862/1‑150 (MQ=255)
                                                                                                         cccGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCGAGCGCGGTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCGTTATCAGTGTGGGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCcaca                                            >  1:123032/1‑150 (MQ=255)
                                                                                                          ccGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCacac                                           >  2:1427885/1‑150 (MQ=255)
                                                                                                             tcggAGCGGCCTTCCCCGCCAACACCGCAAACGCCTGTGCAAACTAGCGCTTCGCCTGCTTCTTGGCCAACCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCa                                        <  2:2432432/148‑1 (MQ=255)
                                                                                                              tggtGCGGGGTGGCCCGCGAACACCGGCACGGCGTGTGCCTTCTAGGGCGTTGCGTGCTTGGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAg                                       >  1:5364705‑M1/5‑150 (MQ=255)
                                                                                                                 AGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAgcgc                                    <  2:4234275‑M1/150‑1 (MQ=255)
                                                                                                                  GCGGGCTTGCCCGCTAACACCGGAAACGCCTGTGCAATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCt                                   >  2:3713859‑M1/1‑150 (MQ=255)
                                                                                                                    GGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGc                                 <  1:3385206‑M1/150‑1 (MQ=255)
                                                                                                                       gTTGCCCGCCAACACGGGAAACGCGTGTGCAATCTCGCGCGTGGCCGGCTTCGGGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGAc                              >  2:2964139‑M1/2‑150 (MQ=255)
                                                                                                                        TTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACg                             <  1:2130‑M1/150‑1 (MQ=255)
                                                                                                                           CCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTc                          <  1:1189691‑M1/150‑1 (MQ=255)
                                                                                                                               ccaaCACCGGAAACGCCTGTGCCATCTCGCGGGTTGCCTGCTTCGTGGCAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCg                      <  2:197112/148‑1 (MQ=255)
                                                                                                                                gAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTGTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGc                     >  2:1059064/1‑150 (MQ=255)
                                                                                                                                gAACACCGCCAACGCCTGTGCCATCTTGCGGGTTGCCTGCTTCGTGGGAAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGc                     <  2:1538893/150‑1 (MQ=255)
                                                                                                                                 aaCACCGGCACCGGGTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGcc                    <  2:1339498/150‑1 (MQ=255)
                                                                                                                                   ccccGGAACCGCGTGGGCCATCTTGGGCGTGGCCTGGTGCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCa                  <  2:179823/148‑1 (MQ=255)
                                                                                                                                   cacCGGAAACGCCTGTGCCATCTCGCGGGTGGCCTGCTTGGTGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCa                  <  2:4117147/150‑1 (MQ=255)
                                                                                                                                      cGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGAt               >  1:818218/1‑150 (MQ=255)
                                                                                                                                      cGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGAt               >  1:23663/1‑150 (MQ=255)
                                                                                                                                        gCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGccc                                                                                                                    <  2:2175419/47‑1 (MQ=255)
                                                                                                                                           aCGCCTGTGCCTTCTCGCGGTTTGCCTGCTTCGTGGGCACGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCa          <  2:1826937/150‑1 (MQ=255)
                                                                                                                                              gCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAAc        >  2:2938800/1‑150 (MQ=255)
                                                                                                                                              gCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAAc        >  2:2630321/1‑150 (MQ=255)
                                                                                                                                              gCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAAc        >  2:2095304/1‑150 (MQ=255)
                                                                                                                                               ccTGGGCCATCTAGCGCGTTGCCTGCTGCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACa       <  2:3988147/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGGCTATCACAGCAAGCTGGCGCGGTTGTTCCGGGGGGTGGAGCAGGTGGCGGCAACGCCCAAGTTCGTCATTCTGAAAGCCCGCAAGTAGGCGCGCCACCGCCCCCGTAGGAGCGGGCTTGCCCGCGAACACCGGCAACGCCTGTGCCATCTAGCGCGTTGCCTGCTTCGTGGGCAAGCCCGCTGCCGCTCTTTATCAGTGTGTGGTCAGGCCCGCCGCCCCCATGAACAGCCGCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAA  >  NC_002947/5275857‑5276153

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