Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,498 1 bp→CG intergenic (+75/‑89) 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,4971.C81.2% 31.1 / 4.6 17intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (6/7);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4980TG94.1% 43.5 / ‑6.8 17intergenic (+75/‑89)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (6/10);  minor base C (1/0);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC  >  NC_002947/1499367‑1499638
                                                                                                                                         ||                                                                                                                                                
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)
                                                                                                                                         ||                                                                                                                                                
GCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCC  >  NC_002947/1499367‑1499638

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: