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.C91.7% 54.8 / 0.5 25intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (9/13);  total (10/15)
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

GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499361‑1499642
                                                                                                                                                        |                                                                                                                                         
gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcg‑‑cc                                                                                                                                        >  2:97398/1‑149 (MQ=37)
  tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcg‑‑cccg                                                                                                                                      >  7:283673/1‑147 (MQ=25)
                  gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCCCGGTGCATGGCACCCGCGTTGCCCGTGTTAGCGGGCAAGCCCGCTCCCACAGACAGAtt                                                                                                                       >  8:245191/1‑149 (MQ=21)
                             ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCg                                                                                                            <  8:283673/149‑1 (MQ=21)
                             ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCCCGAAACACGCGCCGCGGTGCAAGGCACCCGCGTTGACCGTGTTAGCGGGCACCC‑‑CCCCACCCACAGACAGAATGACATTAAGCg                                                                                                            >  8:380041/1‑149 (MQ=255)
                                   cGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCCCAGCCTGCTTGCCATTAAGCGGGCCTa                                                                                                      >  5:317267/1‑149 (MQ=21)
                                     aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                    <  7:380041/149‑1 (MQ=21)
                                                 cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCAtt                                                                                                                 <  1:356361/124‑1 (MQ=18)
                                                 cacaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCCGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCCCTCCCACAGCCGGAATGCCATTAAGCGGGCCTATCGCGACAGTGACg                                                                                        >  7:127304/1‑149 (MQ=18)
                                                   cagctgGACTGAGGCTGCCAAAGGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGc                                                                                      <  6:317267/149‑1 (MQ=21)
                                                     gcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGGGTTGCCGGTGTTAGCGGGCCCGCCCGGGCCCACAGACTGCTTGCCATTAAGCGGGGCTATCGCTACAGTGACGGccc                                                                                    >  6:204379/1‑149 (MQ=21)
                                                      cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCg                                                                                   >  5:133090/1‑149 (MQ=21)
                                                       agGACTGAGGCTGCCAAACGGTGGGAGCGGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATTGCACCGGCGTTGCCGGTGTTAGCGGGCCCGCCCGCTCCCACAGACTGATTTGCATTAATCGGGCCTATCGCTACACTGAAGGcccc                                                                                   >  1:240375/1‑148 (MQ=255)
                                                        ggACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtt                                                                                 <  3:310397/149‑1 (MQ=21)
                                                        ggACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtt                                                                                 <  1:97398/149‑1 (MQ=21)
                                                                    ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCt                                                                      <  5:437357/148‑1 (MQ=21)
                                                                                    gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                         >  5:353754/1‑130 (MQ=21)
                                                                                    gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                         <  6:353754/130‑1 (MQ=21)
                                                                                         tGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGTCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                                                 <  2:280110/149‑1 (MQ=25)
                                                                                          ccccGCTAAACCGCCCCCGCGGTGAAGGCCACGGGCTTTGCGGTTGTCCGCGGCAAGGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                <  3:221257/148‑1 (MQ=255)
                                                                                           cccGCGAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGCGCGGGCACGG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGAGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                <  2:146222/148‑1 (MQ=37)
                                                                                                 aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACAGATTGCCATTAAGCGGGCCTATCGCTa                                                                                                <  4:22372/94‑1 (MQ=17)
                                                                                                 aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACAGATTGCCATTAAGCGGGCCTATCGCTa                                                                                                >  3:22372/1‑94 (MQ=17)
                                                                                                           cgcgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa                               <  2:240375/149‑1 (MQ=37)
                                                                                                                                          ttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGa  <  8:127304/145‑1 (MQ=255)
                                                                                                                                             gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt                          >  7:188077/2‑122 (MQ=255)
                                                                                                                                             gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt                          <  8:188077/121‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                         
GCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499361‑1499642

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: