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.C96.2% 69.3 / NA 27intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (11/14);  minor base A (1/0);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4980TG96.3% 72.3 / ‑6.4 27intergenic (+75/‑89)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (12/14);  total (13/14)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCG  >  NC_002947/1499365‑1499646
                                                                                                                                        ||                                                                                                                                                     
aaGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCAcgcccgc                                                                                                                                         <  3:250800/148‑4 (MQ=25)
   cccGTGATGCGGCCATTCAGAGCCTGCACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAATCGGTGTGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTccc                                                                                                                                      <  5:175979/149‑1 (MQ=21)
   cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTccc                                                                                                                                      <  2:300409/149‑1 (MQ=25)
      gTGATGCGGCCATCCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCcaca                                                                                                                                   >  6:275120/1‑149 (MQ=21)
                     agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCCGT‑‑GTTCGCGGGCACGCCCGCGCCCACAGACTGATTGCCaat                                                                                                                     >  6:183341/1‑146 (MQ=21)
                         ccTCAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTccc                                                                                                                                      <  3:350572/127‑1 (MQ=18)
                                  aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                       >  6:285171/1‑149 (MQ=21)
                                        ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                  <  5:275120/148‑1 (MQ=21)
                                         gCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGt                                                                                                <  1:52367/149‑1 (MQ=21)
                                              acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCAAGGCACCGGCATTTCCGGT‑‑GTTATCGGGCACCCCCCCTCCCACAGCCAGATTGCCAATAAGCGCGCCAATACCTACAGAGAccg                                                                                           >  6:113594/1‑147 (MQ=16)
                                                          aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                                                               <  5:113594/149‑1 (MQ=21)
                                                               gCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCGTTGCCGGT‑‑GTTCGCGGGCACCCCCGCTACCCCAGACTGATGTCCATTAAGCGGGCCTAGCGCTACAGTGACGGCCCGTTTGTTGTTGCCt                                                                          >  1:70685/1‑149 (MQ=18)
                                                                  aaaGGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGc                                                                       <  7:326145/149‑1 (MQ=21)
                                                                  aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCCGTGCATTGCACCCGCTTTGCCGGTGAGTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGGTACAGTGACGGCCCGGTTTTTGTTGGCTCg                                                                        >  1:307937/1‑149 (MQ=21)
                                                                              cgggcgggCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATg                                                           <  2:307937/149‑1 (MQ=25)
                                                                                      cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                    >  1:147046/1‑148 (MQ=37)
                                                                                       ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                    >  6:282803/1‑147 (MQ=37)
                                                                                         gcgAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                  <  1:345167/147‑1 (MQ=37)
                                                                                         gcgAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                  >  2:345167/1‑147 (MQ=37)
                                                                                         gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                                 >  7:158436/1‑148 (MQ=37)
                                                                                           gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc                                              <  8:158436/149‑1 (MQ=37)
                                                                                           gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc                                              <  5:119545/149‑1 (MQ=37)
                                                                                                aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                         <  5:282803/149‑1 (MQ=37)
                                                                                                 ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                        <  2:70685/149‑1 (MQ=37)
                                                                                                       gcgGTGCATGGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                     >  5:234819/1‑146 (MQ=37)
                                                                                                                            gCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtct             >  4:179837/1‑149 (MQ=255)
                                                                                                                             ccGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGActtctt            >  4:89159/1‑149 (MQ=255)
                                                                                                                                  gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa         <  5:285171/142‑1 (MQ=255)
                                                                                                                                  gtgttcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCa         <  4:117411/142‑1 (MQ=255)
                                                                                                                                      tcgcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGACTCGCTAACGAATCATGACTCACTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAc     >  4:91102/4‑149 (MQ=255)
                                                                                                                                         cGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgcg  >  6:192033/1‑149 (MQ=255)
                                                                                                                                         cGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcgcg  >  4:2554/1‑149 (MQ=255)
                                                                                                                                        ||                                                                                                                                                     
AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCG  >  NC_002947/1499365‑1499646

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: