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.C93.8% 89.3 / 0.6 32intergenic (+83/‑81)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (15/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG  >  NC_002947/1499365‑1499649
                                                                                                                                                  |                                                                                                                                               
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                                                     >  2:345167/1‑147 (MQ=37)
                                                                                         gcgAAACAGGCGCCGCGGTGCATTGCACCGGCTTTGCCGGT‑‑GTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                                                     <  1:345167/147‑1 (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)
                                                                                                                                            gCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggcg  >  3:94048/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑‑GTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGCG  >  NC_002947/1499365‑1499649

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: