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% 88.2 / 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 (12/18);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499362‑1499648
                                                                                                                                                          |                                                                                                                                               
cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCAcgccc                                                                                                                                              <  2:128283/149‑2 (MQ=37)
     gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTcc                                                                                                                                         >  5:364679/1‑149 (MQ=25)
      cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTccc                                                                                                                                        <  7:221284/149‑1 (MQ=25)
        cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCcac                                                                                                                                      <  1:114419/149‑1 (MQ=25)
           gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGa                                                                                                                                   >  7:133774/1‑149 (MQ=21)
                            ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACACCCGCCCCCGTGCAAGGCACCCGCTATGCCCGTG‑‑TTAG‑‑CGGGCACCCCCGCTCCCCCAGACAGAATGCCATTAAGCg                                                                                                                  >  6:221192/1‑149 (MQ=255)
                                aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                                               <  8:133774/148‑1 (MQ=21)
                                    aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                                                          >  1:297915/1‑149 (MQ=21)
                                     cgCCTTGGCAGTCGAGCAGGTCTGAGGGTGCCATACGGTGGGAGCGGGGGTGGCCGCGAAACAGGCGCCGGGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                         <  5:221192/148‑1 (MQ=18)
                                     aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                         <  2:297915/149‑1 (MQ=21)
                                     aGCCTTGGCACCACAGCAGGACTGAGGCGGCCAAACGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                                                         <  3:294754/149‑1 (MQ=21)
                                         ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                                                      >  7:2431/1‑148 (MQ=21)
                                           ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCTTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                                                                    <  6:364679/148‑1 (MQ=21)
                                                             aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                                                   <  8:111455/147‑1 (MQ=21)
                                                               gCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGGTTTGGCGGAGGTTTTACCGGGCACCCCCCAACCCACACAATGATTGGCATTAAGCGGGCCTATAGCTACAGTGACGGCCCGttttttt                                                                                  >  1:134977/1‑148 (MQ=255)
                                                                 tGCCAAACGGTGGGAGCGGGCGTGGCCGCGCAACACGCGCCCCGGTGGATTGCACCCGCGTTTCCGGTG‑‑TTAG‑‑CGGGGACGCCCGCGCCCACAGACAGATTGGCAATTAACGGGGCGATAGGGAACATGACGGCCCGGTTTtttttggc                                                                             >  8:14115/1‑147 (MQ=255)
                                                                              gAGGGGGCGGGCCCGCTAAACAGGCCCCGGGGTGCATGGCACGGGCTTTGCGGTTG‑‑TTC‑‑GGGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                                                                <  3:402084/149‑1 (MQ=255)
                                                                               cgCGGGCTTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                               <  8:371006/148‑1 (MQ=25)
                                                                               aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTTAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                                               <  1:11476/149‑1 (MQ=25)
                                                                                     cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                                                         <  7:254113/149‑1 (MQ=25)
                                                                                      gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                                                        <  8:22796/149‑1 (MQ=37)
                                                                                         cccGCTAAACCGGCGCCGCGGTGCAGGGCACGGGCTTGCCGGGTG‑‑TGCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                                                     <  5:40226/149‑1 (MQ=255)
                                                                                          ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                                                      >  8:6996/1‑147 (MQ=37)
                                                                                            gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                                                   >  2:333454/1‑148 (MQ=37)
                                                                                                 aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                                              >  5:90910/1‑148 (MQ=37)
                                                                                                 aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc                                              >  5:131438/1‑148 (MQ=37)
                                                                                                   aGGCGCCGCTGTGCTTGTCAACGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGCGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt                                           <  7:384999/149‑1 (MQ=37)
                                                                                                    ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc                                          <  7:6996/149‑1 (MQ=37)
                                                                                                       gccgcGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                       >  5:389653/1‑149 (MQ=37)
                                                                                                       gccgcGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                       <  8:2431/149‑1 (MQ=37)
                                                                                                           cgGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg                                       >  3:227702/1‑145 (MQ=37)
                                                                                                                                               gcGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGAcg      >  4:149198/2‑148 (MQ=255)
                                                                                                                                                  ggCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACgcggc  <  3:188046/149‑1 (MQ=255)
                                                                                                                                                          |                                                                                                                                               
CTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACGCGGC  >  NC_002947/1499362‑1499648

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: