Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,508 +C intergenic (+85/‑79) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,5061.C89.5% 47.2 / 1.1 19intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (8/9);  total (9/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499367‑1499648
                                                                                                                                                     |                                                                                                                                               
gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcccgctc                                                                                                                                          <  7:146148/148‑6 (MQ=25)
                                aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                           <  8:35790/148‑1 (MQ=21)
                                              agcagGACTGAGGCGGCCAAAGGGTGGGAGGGGGCGTGCCCGCGAAACAGGCGCGGCGGGGCTTGGCACCGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACGGTGACGGc                                                                                             <  4:17057/148‑1 (MQ=21)
                                                             gCCAAACGGTGGGAGCGGGCGTGCCCGCGCAAAAGGCGCCCCCGGGCATGGCACCCGCGTTTCCCGGGTTAGCCGGGACCCCCGCGCCCACAGACAGAATTGCAATAAACGGGGCGATAGCAAAAGTGCCGGCCCGGTTTTTTTTTGc                                                                              >  2:14940/1‑148 (MQ=255)
                                                              ccAAACCGGGGGGCGCGGGGGGCCCCGCGAAACCGGCGCCGCGGGGCAGGCCCCGGGCTTGGCGGGTGGGC‑‑GCGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                                                             <  1:116648/148‑1 (MQ=255)
                                                                      gTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                                                     <  8:119033/148‑1 (MQ=21)
                                                                                   tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                                        >  2:191712/1‑148 (MQ=25)
                                                                                     cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                                                 >  4:67727/1‑90 (MQ=14)
                                                                                      ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                                <  3:67727/90‑1 (MQ=14)
                                                                                              cgggCGCCGGGGTGCTGGGCCCGGGTTTGCCGGGGTT‑‑GCGCGGCCACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt                                             <  2:90485/146‑1 (MQ=255)
                                                                                               aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGGCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                                            >  4:22774/1‑148 (MQ=37)
                                                                                                         ggggcatgCCACGGGCTTGGCCGGTGTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGa                                   <  1:147121/141‑1 (MQ=255)
                                                                                                             aTGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                           >  1:39435/1‑135 (MQ=37)
                                                                                                              tggcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                          <  2:39435/129‑1 (MQ=37)
                                                                                                                gcaccGGCTTTGCCGGTGTTCGCGGGCACGCCCGTTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                           >  7:47231/5‑132 (MQ=38)
                                                                                                                 caccGGCTTTGCCGGTGTTCGCGGGCACGCCCGTTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                          <  8:47231/129‑1 (MQ=38)
                                                                                                                  accGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGt                         >  4:14158/3‑148 (MQ=38)
                                                                                                                                   gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc            >  3:60536/8‑148 (MQ=255)
                                                                                                                                   gtgttcgcGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTcc            <  1:47526/141‑1 (MQ=255)
                                                                                                                                              gCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGCCGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAACCGTCCCGACTTCTTCCATGACgcggc  >  7:51654/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499367‑1499648

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