Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,122,248	C→A	R41L (CGC→CTC)	wcaF ←	colanic acid biosynthesis acetyltransferase WcaF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,122,248	0	C	A	91.7%	30.8 / ‑2.5	12	R41L (CGC→CTC)	wcaF	colanic acid biosynthesis acetyltransferase WcaF
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (2/9);  total (2/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

TAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGC  >  NZ_CP009273/2122167‑2122326                                                                                  |                                                                               tAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAAt                                                                        >  1:37669/1‑90 (MQ=255)        tttttACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTtgtg                                                                 <  1:382174/90‑1 (MQ=255)              cTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGAcc                                                           <  1:355892/90‑1 (MQ=255)                          tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATAttg                                               <  1:2987/90‑1 (MQ=255)                          tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATAttg                                               <  1:282878/90‑1 (MQ=255)                          tAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATAttg                                               <  2:303511/90‑1 (MQ=255)                                       tCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGc                                  <  2:37669/90‑1 (MQ=255)                                                    ccGAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGc                     <  2:152812/90‑1 (MQ=255)                                                      gAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                   <  2:231332/90‑1 (MQ=255)                                                      gAATAAACGTAATAAAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCAc                   <  2:370204/90‑1 (MQ=255)                                                               taataaAAAAGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGc          <  2:218185/90‑1 (MQ=255)                                                                       aaGCCCGCCAGAGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTgccgc  >  2:275417/1‑90 (MQ=255)                                                                                  |                                                                               TAGGTAATTTTTACTGACGGACGAATAACTACGTTTTTTCCTATTTTTGCTCCGAATAAACGTAATAAAAAAGCCCGCCAGCGATACAATACTTGTGGCGACCAGGCAAATATTGTTGCCTGTACTGCCCACCATAATTGCACTTTAATAGCGTTGCCGC  >  NZ_CP009273/2122167‑2122326

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