Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,889,668 C→T G313G (GGC→GGT tnaA → tryptophanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,889,6680CT93.8% 35.0 / ‑3.6 16G313G (GGC→GGTtnaAtryptophanase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (6/9);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTC  >  NC_000913/3889423‑3889893
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
gCTGAAAACGCCTATTTCATCAAGTAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGCAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGtgcgct                                                                                                                                                                                                                              >  1:106500/1‑249 (MQ=255)
     aaaCGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCAGCGAGCAGAACACCCGCGAATCCTACAAATATGCCGATAGGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCTGCCTGCTGTGCATGAAAGACGATAGCTGCTCTGATGTGTACACCGAGTGCAGAACCCGTTGCGTGGTGCAGGAAGGCGGCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGc                                                                                                                                                                                                                         <  1:380687/251‑1 (MQ=255)
                 catcaATCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCGTGCTGTGCATGAAAGACGACAGCTTCTTTGAGGTGTACACCTATTGCAGAACCCTTTGCTTTGTGCGGGAAGGCTTCCCGACATATTGCGGCCTGGAAGGCGGTGCGATGGAGCGTCGGGCGGTAg                                                                                                                                                                                                             <  1:310864/251‑1 (MQ=255)
                                                   aCCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACtggctggc                                                                                                                                                                           <  1:365932/251‑1 (MQ=255)
                                                            cAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGGAAGATGCGCTGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCAGACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTATGGCGGTAGGTCTGTATGACGGCATGaa                                                                                                                                                                                          >  1:66055/1‑227 (MQ=255)
                                                                                                        gTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCTGTGCGATTGAGCGTCTGGCGGTAGGTCTGTATGAc                                                                                                                                                                                                  <  1:375187/175‑1 (MQ=255)
                                                                                                                        gCGATGGTGCCGATGGGCGGCCTGCGGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAAGCCTTTGCGTGGTGCAGGAAGGCTTCCCATCATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGGCTGTATGACGGCATGAATCTCGACGGGCGGGGTTAGCGTATCGCGCAGTTTCAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAgg                                                                                                       <  1:448564/250‑1 (MQ=255)
                                                                                                                                                                                                     ttGCGTGTTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGGCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTc                                                                                                                                            <  1:366170/136‑1 (MQ=255)
                                                                                                                                                                                                                  gAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcgggcg                                                                                                 <  1:308714/166‑1 (MQ=255)
                                                                                                                                                                                                                       cTTCCCGACATATGGCGGCCTGGAAGGCGGTGCGGTGGAGCGTCTGGCGGTAGGTCTGTACGACGCCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTCCAGTATCTGGTAGATGGTCTGGAGGAGATTGGCGTTGTCTGCCCGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCAACTCCCGGCAGACCAGTTCCCGGCCCAGGCGCTGGCCTGCGCGCTGtat        >  1:201343/1‑250 (MQ=255)
                                                                                                                                                                                                                            cGACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTATGGCGGTAGGTCTGTATGACGGCATGa                                                                                                                                                                                           >  1:173770/1‑66 (MQ=255)
                                                                                                                                                                                                                             gACATATGGCGGCCTGGAAGGCGGTGNGATGGANNGTNTGGCGGTAGGTCTGTATGANGGNNNGNNNCNAGACTNGCNNGCTTATCGTNTCGCGCAGGTACAGTATCTGGTTGATGGTCTGGCAGAGATTGGCGTTGTCTGCNAGCcggcgg                                                                                                    >  2:98322/1‑152 (MQ=255)
                                                                                                                                                                                                                             gACATATGGCGGCCTGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTTGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAggcgg                                                                                                    <  1:98322/152‑1 (MQ=255)
                                                                                                                                                                                                                              aCNTATGGNGGCCNGGAAGGCGGTGCGATNGAGCGTCTGGCGGTAGGTCNNTNTTACGNCATGa                                                                                                                                                                                           <  2:173770/64‑1 (MQ=255)
                                                                                                                                                                                                                                          ntGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGgc                           >  1:273164/2‑212 (MQ=255)
                                                                                                                                                                                                                                          cTGGAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACGGGCTGGCTTATCGTATCGCGCAGGTGCAGGATCTGGTCTATGGTCTTGAAGAGAGTGGCGTTGTCTGCCAGCAGGGGGGCGGTCACGCGGCATTCGTTGATGCCGGGAAACTGCTGCCGGATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTc  <  1:230088/237‑1 (MQ=255)
                                                                                                                                                                                                                                            ggAAGGCGGTGCGATGGAGCGTCTGGCGGTAGGTCTGTCTGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGCTGGTCTGGAAGCGATTGGCGTTGTCTGCcag                                                                                                           >  1:64728/1‑130 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
GCTGAAAACGCCTATTTCATCAAGCAGCGTGAAGCAGAATACAAAGACTGGACCATCGAGCAGATCACCCGCGAAACCTACAAATATGCCGATATGCTGGCGATGTCCGCCAAGAAAGATGCGATGGTGCCGATGGGCGGCCTGCTGTGCATGAAAGACGACAGCTTCTTTGATGTGTACACCGAGTGCAGAACCCTTTGCGTGGTGCAGGAAGGCTTCCCGACATATGGCGGCCTGGAAGGCGGCGCGATGGAGCGTCTGGCGGTAGGTCTGTATGACGGCATGAATCTCGACTGGCTGGCTTATCGTATCGCGCAGGTACAGTATCTGGTCGATGGTCTGGAAGAGATTGGCGTTGTCTGCCAGCAGGCGGGCGGTCACGCGGCATTCGTTGATGCCGGTAAACTGTTGCCGCATATCCCGGCAGACCAGTTCCCGGCACAGGCGCTGGCCTGCGAGCTGTATAAAGTC  >  NC_000913/3889423‑3889893

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