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.C87.5% 40.3 / ‑1.2 27intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  major base C (5/16);  minor base A (3/0);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 2.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01
*NC_0029471,499,4980TG96.3% 66.7 / ‑4.0 27intergenic (+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 (9/17);  minor base C (0/1);  total (9/18)
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.75e-01

GGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGG‑‑CA‑‑CGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499378‑1499642
                                                                                                                          ||                                                                                                                                                     
ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGAt                                                                                                                        >  2:110908/1‑149 (MQ=25)
    aTTCAGAGCCTCAACGAACCCTTGCCACCACAGCAGGCCTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGcc                                                                                                                    <  5:232901/149‑1 (MQ=18)
      tCAGAGCCTGAACCAAGCCTTGCCACCACAGCCGGACTGCGGCGGCCAACCGGGGGGGGGGGGCGTGCCCGCTAAACAGGCGCCGGGGTGCAGGGCACGGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCa                                                                                                                   <  6:131971/148‑1 (MQ=16)
                    aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                                           <  7:195855/142‑1 (MQ=21)
                     cgCCTTGCCACCACAGCAGGACTGAGGCTCCCAAACGGGGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAg                                                                                                              <  1:334360/137‑1 (MQ=18)
                     aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                   <  4:305657/149‑1 (MQ=21)
                            gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGTTAGCGGG‑‑CA‑‑CGCCCGCGCCCCCAGACTGATTTCCATTAAGCGGGCCGATAGCTACAGt                                                                                            >  3:265072/1‑149 (MQ=21)
                                 acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGGCTTCCAGCGAAAACCGCGCCGCGGTGCATTGCACCCGCGTTGCCCGTGTTGGCGGGG‑‑CC‑‑CCCCGCGACCACACACTGATTGCCATTAACCGGGCCAATCGCGTCAGAGACgg                                                                                       >  6:272630/1‑149 (MQ=255)
                                         acggcggCTGCCAACGGGTGGGAGCGGGCGGGCCCGCTAAACTGGCGCGGGGGTGCAGGGCACCGGCTTGGCGGGGGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                                               <  4:142990/144‑1 (MQ=19)
                                          cTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGTGGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                              <  7:294001/149‑1 (MQ=21)
                                             atgCTGCGAATCGGTGGGTGCGGGCGTGCCCGCGAATCAGGCGCCGCGGTGCAGGGCTCCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑GGCCCTCTCCCACAGTCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                             <  5:272630/145‑1 (MQ=21)
                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCTCA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                                              <  1:110908/148‑1 (MQ=21)
                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                           >  1:225631/1‑149 (MQ=21)
                                              ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                                          <  3:169986/149‑1 (MQ=21)
                                               gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                                                         >  1:228572/1‑149 (MQ=21)
                                                  gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCt                                                                      >  1:42049/1‑149 (MQ=25)
                                                   ccAAACGGTGGGCGCGGGCGTGCCCGCTAAACAGGGGCCGCGGTGCATGGCACGGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCt                                                                      <  4:64777/148‑1 (MQ=25)
                                                   ccAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                                                     <  4:140810/149‑1 (MQ=21)
                                                                   ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                                                     >  2:162143/1‑149 (MQ=25)
                                                                        gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                >  3:253382/1‑149 (MQ=37)
                                                                         cccGCTAAACTGGCGCCGCGGGGCAGGCCCCCGGCTTTGCCGGTGTTCGCGG‑‑CAACCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                               <  5:241240/149‑1 (MQ=255)
                                                                            gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                             >  6:135857/1‑148 (MQ=37)
                                                                             cgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                                           <  2:228572/149‑1 (MQ=37)
                                                                               aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg                                         <  2:225631/149‑1 (MQ=37)
                                                                                  cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa                                      <  3:326800/149‑1 (MQ=37)
                                                                                       gccgcGTTGCATGGCACCGGCTTTGCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                 <  3:60397/149‑1 (MQ=37)
                                                                                                               gCCGGTGTTCGCGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtc          <  5:135857/148‑1 (MQ=255)
                                                                                                                       ttcgcGGG‑‑CA‑‑CGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGa  >  7:365014/5‑149 (MQ=255)
                                                                                                                          ||                                                                                                                                                     
GGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGG‑‑CA‑‑CGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGA  >  NC_002947/1499378‑1499642

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