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.C82.6% 39.9 / ‑2.4 25intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (7/12);  minor base A (4/0);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 3.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4980TG96.0% 54.0 / ‑3.6 25intergenic (+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 (11/13);  minor base C (0/1);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499355‑1499642
                                                                                                                                                   ||                                                                                                                                                   
gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCCGTGCATGGCACCGGGGTTGCCGGGGtttaacgg                                                                                                                                                 >  1:146222/1‑144 (MQ=38)
   aaCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAc                                                                                                                                              <  1:7079/149‑1 (MQ=37)
      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)
                                                                                                                                                   ||                                                                                                                                                   
GACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499355‑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: