Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,050,838	A→T	59.8%	V186E (GTA→GAA)	glnA ←	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,050,838	0	A	T	59.8%	2.5 / 13.0	15	V186E (GTA→GAA)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base T (4/5);  total (6/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

AGCTTTTTTGGTCATGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGAACCGGAGATAGATGATCCGAAACGGA  >  CP009273/4050696‑4050979                                                                                                                                               |                                                                                                                                              aGCGTTATTGGTCATGGTATTGACGCGGGTAGCCACTTCGTTCTGACGAGCAGTCGCTACTTCGTGGCGATGGGCTTCAACCACCAGACCCATCTGTTCCATCTCCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTgg                                                                                                                                        <  1:150059/150‑1 (MQ=255)                             gTAGCCACTTCGTTCTTACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGcc                                                                                                            <  1:253367/150‑1 (MQ=255)                                            tGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGCAGGTGAca                                                                                                 >  2:101594/1‑145 (MQ=255)                                                cagcagTCGCTACTGCGTGGTGAGGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGCCTACCGGTGGAACCGGTAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTac                                                                                         <  1:105961/150‑1 (MQ=255)                                                    agTCGCTACTTCGTGGTTATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTaccacc                                                                                     <  1:218977/150‑1 (MQ=255)                                                               cGTGGTGATGGGTTTCAACCACCGGACCCTTCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCCTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTggg                                                                          <  1:241162/150‑1 (MQ=255)                                                                     gATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTggagg                                                                    <  1:54993/150‑1 (MQ=255)                                                                                         cccATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTTCCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCttg                                                    >  2:7030/1‑145 (MQ=255)                                                                                          ccATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCg                                                   <  2:55453/146‑1 (MQ=255)                                                                                                   ccATCACCAGACACATTTCAGAACGAATATCCTGAGCCGATTCTTCCGGTGGAACCGGTAATTAACCGCCTTTCACTGCCTGACGTTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAg                                      <  1:222673/150‑1 (MQ=255)                                                                                                             acacATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGa                            <  1:8196/150‑1 (MQ=255)                                                                                                                 aTTTCAGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGAGCCTTCGCTATCGTCGAGCGCAACGTggggagcg                        >  1:118007/1‑145 (MQ=255)                                                                                                                      aGAACGAATATCCTGAGCCGAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGAACCgg                       >  2:142523/1‑146 (MQ=255)                                                                                                                                         gAGTCTTCCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGAACCGGAGCTAGATGATCCGAACCg    >  2:13801/1‑146 (MQ=255)                                                                                                                                           gTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGAACCGGAGATAGATGACCCGAAAAGGa  >  2:31466/1‑146 (MQ=255)                                                                                                                                               |                                                                                                                                              AGCTTTTTTGGTCATGGTATTGAAGCGGGTAGCCACTTCGTTCTGACCAGCAGTCGCTACTTCGTGGTGATGGGCTTCAACCACCAGACCCATCTGTTCCATCACCAGACACATTTCAGAACGAATATCCTGAGCCGAGTCTACCGGTGGAACCGGGAAGTAACCGCCTTTCACTGCCGGACGGTGACCTTTGTTACCACCTTCGTATTGGGTGGAGGAGTTCCATGCGCCTTCGATATCGTCGATAGCAACGTGGGAACCGGAGATAGATGATCCGAAACGGA  >  CP009273/4050696‑4050979

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