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.C89.7% 73.7 / 3.9 31intergenic (+74/‑90)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (10/16);  total (13/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 = 1.00e+00
*NC_0029471,499,4980TG96.7% 78.3 / ‑3.9 31intergenic (+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/17);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG  >  NC_002947/1499362‑1499644
                                                                                                                                             ||||                                                                                                                                                     
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)
                                                                                                                                             ||||                                                                                                                                                     
CTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCTGAAACGTCCCGACTTCTTCCATGACG  >  NC_002947/1499362‑1499644

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: