Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,751,237	A→G	V241V (GTT→GTC)	BW25113_RS13685 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,751,237	0	A	G	92.3%	32.3 / ‑3.0	13	V241V (GTT→GTC)	BW25113_RS13685	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (2/10);  total (2/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTC  >  NZ_CP009273/2751150‑2751312                                                                                        |                                                                            ctcAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGAc                                                                           <  1:9898/90‑1 (MQ=255)         ctGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCat                                                                   <  1:18811/90‑1 (MQ=255)         ctGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCat                                                                   >  1:90685/1‑90 (MQ=255)                           cccATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACt                                                 <  2:290543/90‑1 (MQ=255)                                aGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGt                                            <  2:337581/90‑1 (MQ=255)                                     ttCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAAt                                       <  1:383044/90‑1 (MQ=255)                                       cGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGa                                     <  1:248469/90‑1 (MQ=255)                                       cGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGa                                     <  2:371978/90‑1 (MQ=255)                                        gATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGat                                    <  1:126565/90‑1 (MQ=255)                                                             aaCCTTATTTTTCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCgctttgct               <  2:71001/90‑1 (MQ=255)                                                                       ttCCACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCAccc     >  2:85499/1‑90 (MQ=255)                                                                          cACCCATCATTCGGGACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCcttc  <  1:85499/90‑1 (MQ=255)                                                                          cACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCcttc  <  1:4374/90‑1 (MQ=255)                                                                                        |                                                                            CTCAAACTCTGACATTTGAATAAAATCCCATAGCGGTTCGATAAGTTCATTAACTGCTACAACCTTATTTTTCCACCCATCATTCGGAACAGAGTCATTAATTAGCTCAACCAACTTCAGTTGAATGATATGGTACACTTCCGCTTTGCTTTTGCCACCCTTC  >  NZ_CP009273/2751150‑2751312

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