Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	534,609	G→A	V232I (GTT→ATT)	gcl →	glyoxylate carboligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	534,609	0	G	A	90.9%	22.5 / NA	11	V232I (GTT→ATT)	gcl	glyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (5/5);  minor base C (1/0);  total (6/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 = 1.00e+00

GAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGT  >  NC_000913/534423‑534854                                                                                                                                                                                           |                                                                                                                                                                                                                                                      gAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCCGATCGACAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCCCTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTCTCCCGGACGATCATGAACTGATGGCCgg                                                                                                                                                                                       >  1:114231/1‑251 (MQ=255)         gTTTGATCCTGACATGTGCGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAAGGTTAATCCAGGCCGAACGTCCGGTGATTTTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATTGGCTGGGGCTGTATCCCGGACGATCATGAACTGAGGGCCGGGAtggtgg                                                                                                                                                                               <  1:124354/251‑1 (MQ=255)           ttGATCCTGACATGTTCGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTggg                                                                                                                                                                              <  1:7067/250‑1 (MQ=255)              aTCCTGACATGTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTAAGCAGCTAGAAAAAGCTGTCGAAATGTTAATCCAGGCCGACCGTCCGGTGATTGTTGCCGGGGGCGGGGTAAATAATGCTGACGCAGCTGCACAGTTACAACAGTTTGCTGAACTGAACAGCCTTCCGGAGAACCCAACGATAATGGGCTTGGGCTAAATACCGGACGATCATGAACTGATGGACGGGATGGTGGGtatg                                                                                                                                                                          >  1:120330/1‑248 (MQ=255)                     cATGTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCgg                                                                                                                                                                                                             <  1:69874/209‑1 (MQ=255)                        gTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCTGGGTAATTAATGCTGACGCAGCTTCACTGTTACAACATTTTGCTGAACTGACCAGCATTCCGGTGATCTCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACcgcn                                                                                                                                                                 <  1:32989/250‑2 (MQ=255)                                                                                     gAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGCAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGttt                                                                                                   >  1:54960/1‑251 (MQ=255)                                                                                                                                   gggTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCATTCCGGTGATCCCAACGCTCATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGgcgc                                                     >  1:54061/1‑251 (MQ=255)                                                                                                                                                        cTGCANNGTNNNAANNGTTTGCTGAACTGACCAGCATTCNGGTGATCCCAACGCTAATGGGCTGGGGCTNTNTNCAGGNNGATNCTGAACTGATNGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGNTNCTGGCGTCTGACATGGTGTTTGGTATCGGTANCCGTTTTGCTAACCGTCATACCGGNTCGGNnggn                                                                       >  2:127573/1‑210 (MQ=255)                                                                                                                                                                            gCTGAACTGACCAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTGCATATTGATATTGAGCCGACGCAAATTGGTc             <  1:126592/250‑1 (MQ=255)                                                                                                                                                                                      ccAGCATTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTgtgt  >  1:128892/1‑251 (MQ=255)                                                                                                                                                                                           |                                                                                                                                                                                                                                                      GAAATCGAGTTTGATCCTGACATGTACGAACCGCTGCCGGTCTACAAACCTGCTGCCAGCCGTATGCAGATCGAAAAAGCTGTAGAAATGTTAATCCAGGCCGAACGTCCGGTGATTGTTGCCGGGGGCGGGGTAATTAATGCTGACGCAGCTGCACTGTTACAACAGTTTGCTGAACTGACCAGCGTTCCGGTGATCCCAACGCTAATGGGCTGGGGCTGTATCCCGGACGATCATGAACTGATGGCCGGGATGGTGGGTCTGCAAACCGCGCATCGTTACGGTAACGCAACGCTGCTGGCGTCTGACATGGTGTTTGGTATCGGTAACCGTTTTGCTAACCGTCATACCGGCTCGGTAGAGAAATACACCGAAGGGCGCAAAATCGTTCATATTGATATTGAGCCGACGCAAATTGGTCGCGTGCTGTGT  >  NC_000913/534423‑534854

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