Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,196,232	C→T	44.9%	T370T (ACC→ACT)	icd →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,232	0	C	T	44.9%	‑4.7 / 10.8	14	T370T (ACC→ACT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (3/4);  new base T (4/2);  total (8/6)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

CTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAT  >  NC_000913/1196093‑1196374                                                                                                                                            |                                                                                                                                               cTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATACGGCGGGCAGGACAAAGTAAATCCTGGCTCTGTTGTTCTCTCCGCTGAGAGGAGGCTGCGCCACAGGGGTTGGGCCGaa                                                                                                                                             >  1:194366/1‑143 (MQ=255)       cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCt                                                                                                                                       <  1:196937/143‑1 (MQ=255)                         gcgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCCGCTGACCTGATGGGTAAAgtg                                                                                                                     >  2:483412/1‑141 (MQ=255)                             ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATGTGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACAGAAGCGGCTGACCTGATTGTTCAAGGTATgg                                                                                                                 >  2:31272/1‑143 (MQ=37)                             ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTCAAGGTATgg                                                                                                                 >  1:134452/1‑143 (MQ=37)                                gTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAg                                                                                                              <  2:520629/143‑1 (MQ=37)                                       gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTAGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGCCCGGAGCGGCTGACTTAATGGTCCAAGGTATCGACGGCGCaac                                                                                                       >  1:156997/1‑142 (MQ=255)                                          acccacGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCAGATTGTTAAAGGTATGGAAGGCGCAATCaa                                                                                                    >  2:501417/1‑143 (MQ=25)                                          acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                                                                                                    <  1:78316/143‑1 (MQ=255)                                                     tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATACGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGAcgg                                                                                         >  2:43404/1‑141 (MQ=18)                                                                                 gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTCACCTATGGCCTCGGGCGTCTgatggc                                                             >  2:497533/1‑142 (MQ=255)                                                                                 gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTgatgga                                                             <  1:329272/143‑1 (MQ=255)                                                                                               tattctTCTCTCCGCTGAGATGAGGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAACAGTATGGAAGGCTCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgc                                               <  1:204706/143‑1 (MQ=255)                                                                                                     tctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                                         <  1:359106/143‑1 (MQ=255)                                                                                                                                            cGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt  >  1:228028/1‑143 (MQ=255)                                                                                                                                            cGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGACCATCGACAACATGTAAt  >  2:40119/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                               CTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAT  >  NC_000913/1196093‑1196374

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