Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	503,429	C→T	intergenic (‑191/+47)	ybaL ← / ← fsr	putative transporter YbaL/fosmidomycin efflux pump

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	503,429	0	C	T	90.0%	22.0 / ‑6.1	10	intergenic (‑191/+47)	ybaL/fsr	putative transporter YbaL/fosmidomycin efflux pump
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/4);  minor base A (1/0);  total (6/4)
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

GATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCT  >  NC_000913/503215‑503675                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gATAAGCGGGGTGGCGTGATGCATTCCGTCGCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGTAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTgg                                                                                                                                                                                                                     <  2:362145/250‑1 (MQ=255)        ggggTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATg                                                                                                                                                                                                             <  2:253934/251‑1 (MQ=255)                                                             ggCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAAGTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTTTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAAt                                                                                                                                                         <  1:100244/250‑1 (MQ=255)                                                                              ttACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGAGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAATcac                                                                                                                                         >  2:324176/1‑249 (MQ=255)                                                                                         aGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGa                                                                                                                             <  1:126925/250‑1 (MQ=255)                                                                                                                                            aGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGATATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGAATTTATAGACTAACTCGATGCTGGTGGGATCGGCGATAAACCCCAGAACTGCCGCcccccg                                                                         >  2:8190/1‑249 (MQ=255)                                                                                                                                                             tGCAAAAAAAGCACGGATCGGAAAACAGAAAGCGCTTTAACAGAGGAAAGGCTTATGGTGACTGTCAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGCAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGACTTTATAGACTAACTCGATGCTGGTGTGCTCGGCGATAAGCCCCAGAACTGCCGCTCCCCGACCTACCCTGCCaaa                                                          >  1:401172/1‑249 (MQ=255)                                                                                                                                                                           ggATCGGAAAACAGCAGGCGGTTTAACAGAGGCACGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGACAGCACCGATTTTATAGACTAACTCGATGCTGGTGTGCTCGGCGCTAAGACCCCGCACTGCCGCCCCCCGACCTCCCCTGCCCAAAGCCAAACCGcaaa                                            >  2:172791/1‑249 (MQ=255)                                                                                                                                                                                                         ggAAAGGCTTATGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCCCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTgg             >  2:128080/1‑250 (MQ=255)                                                                                                                                                                                                                   aTGGTGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGGGATCGGCGCTAAGCCCCAGAACTGCCGCTCCCCGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTccc  >  2:12237/1‑250 (MQ=255)                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GATAAGCGGGGTGGCGTGATGCATTCCGTCTCCTTTTCCTGGTGGTTATTGTCCATTTTTGGCCGGGAAAACCAAAATTACAGGTAATAGTTTATGACAATTTCATTGATGATGTTCATGAATAATTATTGAATTTTGCAGAAAAATGGAATTAGCTGCAAAAAAAGCACGGATCGGAAAACAGAAGGCGGTTTAACAGAGGAAAGGCTTATGGCGACTGTAAAACGATGCAGCCAAAGTTTGCCTTTGGCTGCAATGAAATCAGTCTTTATGCCGGTTATCAGGCAGGAATATGGTCAACATCCCCAATAGTGGCAGGAAAGCACAGATTTTATAGACTAACTCGATGCTGGTGTGATCGGCGATAAGCCCCAGAACTGCCGCTCCCAGACCTCCCATGCCAAAAGCAAAACCGAAAAAGAGTCCAGAAACCATACCGATACGTCCTGGAAGCAGCTCCT  >  NC_000913/503215‑503675

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