Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,673,046	C→T	V211V (GTC→GTT)	rcdB →	putative DNA‑binding transcriptional regulator YhjC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,673,046	0	C	T	96.8%	95.9 / ‑6.0	31	V211V (GTC→GTT)	rcdB	putative DNA‑binding transcriptional regulator YhjC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/9);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGGCGTGTAAACCTGTT  >  NC_000913/3672800‑3673270                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  ccGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCTGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGt                                                                                                                                                                                                                               >  1:269047/1‑251 (MQ=255) ccGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAgg                                                                                                                                                                                                                               >  1:62558/1‑250 (MQ=255) ccGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGt                                                                                                                                                                                                                               <  1:310656/251‑1 (MQ=255) ccGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGt                                                                                                                                                                                                                               >  1:335018/1‑251 (MQ=255)             tatTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTg                                                                                                                                                                                                                  >  2:136524/1‑251 (MQ=255)                    cATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCAt                                                                                                                                                                                                           >  1:314123/1‑251 (MQ=255)                                   tGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCGCTGTCTGACGCGCTTTGGTT‑ATCCGCAAAGGCCCGATGATCTGACTTCGCTCGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCTAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACa                                                                                                                                                                                             <  1:382608/250‑1 (MQ=255)                                    gCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAg                                                                                                                                                                                            <  1:58426/250‑1 (MQ=255)                                           aCGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGaa                                                                                                                                                                                     <  2:120580/250‑1 (MQ=255)                                            cGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGaaaa                                                                                                                                                                                   <  1:249529/251‑1 (MQ=255)                                                                    ggCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCgg                                                                                                                                                           >  1:288940/1‑251 (MQ=255)                                                                      cGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTc                                                                                                                                                         >  1:223525/1‑251 (MQ=255)                                                                                    cccttcggaaaCTGGACCAGGGTCAACGGGTCCAGTCCCCACTATCTGACGCGCTTTGTTT‑ACCCTAAAAGCCCCGAGTATCCGACTTCACACGCAATAGTGCGTTATCCCCCGCCCCTGGGTGTACCTCCGTTAGTTTGTGAGGTTGCCAGCTTTAATGGCGTTCATTGGTTAAA‑TTTTGGCGGGATGTTGACGGTAAAAAGTAGCAAAAACTATCGCACCGCCGGTCTTGCCGGTCTGGGGATTATTCa                                                                                                                                            <  1:116289/243‑1 (MQ=255)                                                                                         ggCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATccc                                                                                                                                      >  1:243848/1‑251 (MQ=255)                                                                                            aaaCTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCcgcg                                                                                                                                   >  2:309017/1‑251 (MQ=255)                                                                                                tG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCAtt                                                                                                                               >  2:321344/1‑251 (MQ=255)                                                                                                     cATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACCCCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCg                                                                                                                           >  2:167132/1‑251 (MQ=255)                                                                                                      ttGGTCAATTTTGCCAGCCCGCACTATTGGACGCG‑TTTGGTTTATCCGCAAGCCCCCGATGATCTGACTTCACACGAAATAGTGCGTTATCCCCCGCCCCTGGGTGTACATCCTTTAGTTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGTTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAg                                                                                                                                           <  2:157969/233‑1 (MQ=255)                                                                                                      aTGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAg                                                                                                                                           >  1:157969/1‑234 (MQ=255)                                                                                                          tCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCAGCAAAGCCACGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCCGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTgcgc                                                                                                                      >  2:29832/1‑251 (MQ=255)                                                                                                                   ccAGTCCGCCCTATCTGCCGCGTTTTGGTT‑ACCCGCAAACCCCCGATTATCTGACTTCACACGCAATAGTGCGTTATCCCCCGCACCTGGGTGTACATCCGTTGGGTTTGGAGGTTGCCAGCGTTAAGGGCGTTCAGTGTTTTAA‑GTCTGGCGGAATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTg                                                                                                             <  2:335018/251‑1 (MQ=255)                                                                                                                           cACTATCTGACGTGTTTTGTTT‑ACCCAAAAAGCCCCTATGATCTGACTTCACCCGAAATAGTGCGTTATACCCCGCCCCTGGGTTTACATCCGTTAGTTTTTGAGTTTGCCAGCGTTAATGGCGTTCAG‑GGTTTAAATTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCCCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCgg                                                                                                     <  2:217763/251‑1 (MQ=255)                                                                                                                              tATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGggcg                                                                                                  >  1:87405/1‑251 (MQ=255)                                                                                                                                tCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGggcggc                                                                                                >  1:360215/1‑251 (MQ=255)                                                                                                                                       cgcTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTg                                                                                          >  1:62973/1‑250 (MQ=255)                                                                                                                                                               cGATGATCAGACTTCCCCCTCAATAGTGCGTTATCCACCGCACCTGGGTTTTCATCCGTTAGTTTTTGAGGTTGCTAGCGTTATTGGCGTCCATTGGTTTAA‑GTTGGGCGGTTTGTGGTCGGTAAACATTGGCGAAAACTATCCCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCg                                                                  <  1:227981/251‑1 (MQ=255)                                                                                                                                                                            tCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCt                                                     >  2:397162/1‑251 (MQ=255)                                                                                                                                                                                  gCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTg                                                                     >  1:126732/1‑229 (MQ=255)                                                                                                                                                                                     aTAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGt                                            >  2:251461/1‑251 (MQ=255)                                                                                                                                                                                                                           aGGTTTTGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGGCGTGTAAAcc      >  1:322459/1‑251 (MQ=255)                                                                                                                                                                                                                               tttGAGGTTGCCAGCGTTAATGGCGTTCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGGCGTGTAAACCTGtt  >  1:297125/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  CCGTCCGGTAGATATTCTTCATGATGGTTTTGATTGCGTGATACGCACTGGCGCGTTACCGGAAGATGGCGTTATCGCCCGTCCCCTCGGCAAACTG‑ACCATGGTCAACTGTGCCAGTCCGCACTATCTGACGCGCTTTGGTT‑ATCCGCAAAGCCCCGATGATCTGACTTCACACGCAATAGTGCGTTATACACCGCACCTGGGTGTACATCCGTTAGGTTTTGAGGTTGCCAGCGTTAATGGCGTCCAGTGGTTTAA‑GTCTGGCGGCATGTTGACGGTAAACAGTAGCGAAAACTATCTCACCGCCGGTCTTGCCGGTCTGGGGATTATTCAGATCCCGCGCATTGCCGTGCGCGAAGCCCTGCGTGCCGGGCGGCTTATTGAAGTATTACCTGGCTACCGTGCCGAGCCGCTCTCCCTTTCGCTGGTTTATCCGCAGCGTCGGGAGCTTTCCCGGCGTGTAAACCTGTT  >  NC_000913/3672800‑3673270

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑