Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,800,111	C→A	26.8%	E75* (GAA→TAA)	ansB ←	glutaminase‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,800,111	0	C	A	26.8%	17.1 / 5.8	15	E75* (GAA→TAA)	ansB	glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base A (2/2);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

GATGCCGTCCACGTCCTTGCTTTCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGACCACGTTGGCCAGCTTCTGTTGGGTTTCGGC  >  NC_002947/2799977‑2800252                                                                                                                                       |                                                                                                                                              gATGCCGTCCACGTCCTTGCTTTCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGAt                                                                                                                                <  2:111338/150‑1 (MQ=255)                     tttCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCt                                                                                                            >  1:176341/1‑150 (MQ=255)                     tttCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTACGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCt                                                                                                            >  2:163322/1‑150 (MQ=255)                     tttCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTACGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCt                                                                                                            >  2:407288/1‑150 (MQ=255)                                   gCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTaa                                                                                              >  1:192553/1‑150 (MQ=255)                                    cGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTACGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAg                                                                                             <  1:163322/150‑1 (MQ=255)                                    cGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTACGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAg                                                                                             <  1:407288/150‑1 (MQ=255)                                                  aGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGc                                                                                                                  >  1:253793/1‑115 (MQ=255)                                                  aGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGc                                                                                                                  <  2:253793/115‑1 (MQ=255)                                                                   aTTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGTCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGcc                                                              >  1:268948/1‑150 (MQ=255)                                                                                              cATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGAt                                   >  2:465002/1‑150 (MQ=255)                                                                                                 cACTTGCTCGACACGTGCATTGGAGATGTCGGACACTTCCGTCACGCCTTCTATAAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGAc                                <  2:53916/150‑1 (MQ=255)                                                                                                                  cATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGACCACGTTGGCCAGCTTCt               <  2:491841/150‑1 (MQ=255)                                                                                                                            gTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGACCACGTTGGCCAGCTTCTGTTGGGTTTc     <  2:268948/150‑1 (MQ=255)                                                                                                                               ggCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGACCACGTTGGCCAGCTTCTGTTGGGTTTCGGc  >  1:325958/1‑150 (MQ=255)                                                                                                                                       |                                                                                                                                              GATGCCGTCCACGTCCTTGCTTTCGGCCAGCTCGGCGACGCGCTTGCCCAGCTTCAGCAGGTCGTCATTGCTGATGCTCTCGGAGGCGATCTGCATCACTTGCTCGCCACGTACATTGGCGATGTCGGCCAGTTCCGGCACGCCGGCGATGAGCTTGTCGACGCCCAGCTTGGCAGCCTGGTAAGTGGCGCTGTTGGCAGCGCTGGCACCGGCGCCGGCAATGGTGCCGCCGGTGGCGAGGATGACCACGTTGGCCAGCTTCTGTTGGGTTTCGGC  >  NC_002947/2799977‑2800252

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