Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,679 C→T A216T (GCC→ACC)  yjhU ← putative DNA‑binding transcriptional regulator YjhU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,6790CT92.9% 27.1 / ‑3.6 14A216T (GCC→ACC) yjhUputative DNA‑binding transcriptional regulator YjhU
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (7/6);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TTCCCCACCAGCAAGGAAGATCTTCCGTTCAATCTGCCTTAGCTGAGCCAGACTGATTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGCTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCACGTGCAGCGAGGCGACACGGCCTGAAAGCTTGTAGGCGATTTGCGATGAGGATTCAATGATCGTCAGGTCTGGCTGAGATGGCATGGAACCGAGCATCTGTATAACGGTGACATTTTTACATGACTTCGGCGTCATGGTTTTGCTCATCTG  >  NC_000913/4519446‑4519925
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
ttCCCCACCAGCAAGGAAGATCTTCCGTTCAATCTGCCTTAGCTGAGCCAGACTGATTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGGCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTGGTACAGCGAg                                                                                                                                                                                                                                          <  1:74332/248‑1 (MQ=255)
ttCCCCACCAGCAAGGAAGAGCGTCCGTTCAATCTGCCTTAGCTGTGCCAGACTGATTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCATCACGCCCACCCTGATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGCTTTCGTACAGCGAGtn                                                                                                                                                                                                                                        <  1:38994/250‑2 (MQ=255)
                                         gCTGAGCCAGACTGCTTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAATtcgat                                                                                                                                                                                                                      >  1:110992/1‑227 (MQ=255)
                                                                               tCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGgc                                                                                                                                                                             >  1:13204/1‑230 (MQ=255)
                                                                                                               gAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGCGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTAGTACAGc                                                                                                                                                                                                                                             <  1:221957/134‑1 (MQ=255)
                                                                                                               gAAGCGAANANNNANNGCACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTNTTCAGCACGCCCACCNGCNTCAGNGGTNTTTCATCAAGCGCGTTACCAACGACAAAAAAGGTTTTCGTACAGc                                                                                                                                                                                                                                             >  2:221957/1‑134 (MQ=255)
                                                                                                                         caAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCtt                                                                                                                                                                         >  1:216861/1‑192 (MQ=255)
                                                                                                                                cacCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATGTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCGCTTCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCTCGTGCAGCGAGGCGACACGGCCTGAAAGCTTGTAGGCGATTTGCGATGAg                                                                                                       <  1:141565/251‑1 (MQ=255)
                                                                                                                                                                                  aCGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGc                                                                                                                                                                                                          <  1:326468/102‑1 (MQ=255)
                                                                                                                                                                                   cGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTACAGCGCGTCAGACCGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCACGTGCAGCGAGGCGACACGGCCTGCAAGCTTGTAGGCGATTTGCGATGAGGATTCAATGATCGTCAGGTCTGGCTGAGATggggatcggaagggcacacgt                                                    >  1:47329/1‑232 (MQ=255)
                                                                                                                                                                                        aCCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGTTTTCGTCCAGCGAGTCAGAACGTCAAAAt                                                                                                                                                                                                                           >  1:159411/1‑79 (MQ=255)
                                                                                                                                                                                                                    cgTTACCCNCGACNANNNAGGTTNTCGTACAGCNAGTCAGAACGTCAAAATTCNATCGaa                                                                                                                                                                                                                  >  2:31783/1‑60 (MQ=255)
                                                                                                                                                                                                                    cgTTACCCACGACAAAAAAGGTTTTCGTACAGCGAGTCAGAACGTCAAAAGTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCACGTGCAGCGAGGcn                                                                                                                                           <  1:31783/131‑2 (MQ=255)
                                                                                                                                                                                                                                     aaGGTTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCACGAGCAGCGAGGCGACACGGCCTGAAAGCTTGTAGGCGATTTGCGATGAGTATTCAATGATCGTCAGGTCTGGCTGAGATGGCATGGAACCGAGCATCTGTATAACGGTGACATTTTTACATGACTTCGGCGTCATGGTTTTGCTCATCTg  <  1:323629/251‑1 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
TTCCCCACCAGCAAGGAAGATCTTCCGTTCAATCTGCCTTAGCTGAGCCAGACTGATTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCATTTCCTTCTTATTCAGCACGCCCACCCGAATCAGGGGATTTTCATCAAGCGCGTTACCCACGACAAAAAAGGCTTTCGTACAGCGAGTCAGAACGTCAAAATTCGATCGAATAATAGGCTCTGCCTGTAACTCCATGGCCAGTCTGGCGCTTGATACCACGGCGGGCACGTGCAGCGAGGCGACACGGCCTGAAAGCTTGTAGGCGATTTGCGATGAGGATTCAATGATCGTCAGGTCTGGCTGAGATGGCATGGAACCGAGCATCTGTATAACGGTGACATTTTTACATGACTTCGGCGTCATGGTTTTGCTCATCTG  >  NC_000913/4519446‑4519925

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