Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	885,533	C→A	39.9%	intergenic (‑258/‑12)	ybjL ← / → ybjM	putative transporter/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	885,533	0	C	A	39.9%	8.9 / 2.8	10	intergenic (‑258/‑12)	ybjL/ybjM	putative transporter/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (2/2);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGA  >  CP009273/885288‑885731                                                                                                                                                                                                                                                      |                                                                                                                                                                                                       cttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTAAGCGCACCTAAACGATAGTAAAGTAACGATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGt                                                                                                                                                                           <  1:43717/275‑1 (MQ=255)                                                                                                                                    aTTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTAAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCtt                                                                                               >  1:54769/1‑220 (MQ=255)                                                                                                                                        ggTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTAAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCg                                                                                           <  2:54769/220‑1 (MQ=255)                                                                                                                                              aaCGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCAAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTNTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGAATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGCAGCAGTCGCCAGc                              >  1:509949/1‑275 (MQ=255)                                                                                                                                                                         tatTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTAAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTc      >  2:183577/1‑272 (MQ=255)                                                                                                                                                                             aGCATTGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTAAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGa  <  1:183581/272‑1 (MQ=255)                                                                                                                                                                                              ggCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCagg                                                                                                                                                                                                     >  2:591931/1‑60 (MQ=255)                                                                                                                                                                                                  ttATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAggagga                                                                                                                                                                                                 <  1:591942/60‑1 (MQ=255)                                                                                                                                                                                                               cGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTg                                                                                                                     >  1:114235/1‑123 (MQ=255)                                                                                                                                                                                                                   acacTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGa                                                                                                                 <  2:114234/123‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                       CTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAAAATGTGTTCGATAGACGCAGTGTCAGGAGGAACGAGTGAAACATAAACAACGTTGGGCGGGGGCAATCTGCTGTTTTGTCCTCTTCATTGTGGTGTGCCTTTTTCTGGCGACGCACATGAAAGGCGCTTTTCGGGCTGCCGGGCATCCTGAAATCGGCTTGCTATTTTTCATTCTTCCTGGAGCAGTCGCCAGCTTCTTTTCACAGCGTAGAGAAGTCCTGA  >  CP009273/885288‑885731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑