Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	805,853	Δ1 bp	intergenic (+89/+145)	ybhJ → / ← ybhC	putative hydratase YbhJ/outer membrane lipoprotein YhbC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,853	0	T	.	93.8%	28.0 / ‑1.8	16	intergenic (+89/+145)	ybhJ/ybhC	putative hydratase YbhJ/outer membrane lipoprotein YhbC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (10/5);  minor base C (0/1);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATT  >  NC_000913/805723‑805896                                                                                                                                   |                                            gCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGccg                                   >  1:483037/1‑139 (MQ=255)  cGGGTAGTTTTGTTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGGTGCGAAGCTTGCGCGTATAATCAGACCTACAGATCGGATCGAA‑CCGTAGGccgg                                  >  2:215435/1‑138 (MQ=255)    ggTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGAt                                >  2:433845/1‑140 (MQ=21)    ggTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGAt                                >  1:79820/1‑140 (MQ=21)    ggTAGTTTGATTAACCTCAATAAAACCCGACAGATGTCAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGAt                                <  1:211040/140‑1 (MQ=18)       aGTTTGATTAACTTCAAGAAAAACCGTCAGACGTACAAACCGCCATGTCAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGCCCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAg                             <  2:223948/140‑1 (MQ=18)                aaCTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAACC‑GTAGGCCGGATAAGGCGTTTACg                    <  2:646949/140‑1 (MQ=17)                aaCTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGtttacg                    <  1:139220/140‑4 (MQ=18)                  cTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATGCGACCTACTAATAGCAGCGAA‑CCGTAGGCCAGATAAGGCGTTTACgcc                  >  2:206854/1‑139 (MQ=18)                   ttCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACgccg                 >  2:39097/1‑138 (MQ=17)                   ttCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACgccg                 >  2:584388/1‑138 (MQ=17)                           aaaCCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc         <  2:160152/140‑1 (MQ=14)                           aaaCCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc         <  2:57547/140‑1 (MQ=14)                            aaCCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa        >  2:164403/1‑140 (MQ=14)                             aCCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa       >  1:412193/1‑139 (MQ=14)                             aCCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaa       >  2:445114/1‑139 (MQ=14)                                  cAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACAAATCGCATCGAA‑CCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaatagt  >  1:429904/1‑134 (MQ=11)                                                                                                                                   |                                            GCGGGTAGTTTGATTAACTTCAATAAAAACCGTCAGATGTAAAAAGCGCCATGTGAATGTAGGTCGCATTCGGCACTTATTGTCGGATGCGATGCTTGCGCATCTTATCCGACCTACGAATCGCATCGAATCTGTAGGCCAGATAAGGCATTTTCGCAGCATCCGGCACTTATT  >  NC_000913/805723‑805896

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