Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,196,220	C→T	44.5%	H366H (CAC→CAT)	icd →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,196,220	0	C	T	44.5%	‑0.3 / 10.5	16	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (4/5);  new base T (5/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAT  >  NC_000913/1196082‑1196335                                                                                                                                           |                                                                                                                    cggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                                                                 <  1:157305/143‑1 (MQ=255)        gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATACGCCGGTCAGGACAAAGTAAATCCTGGCTCGATTATTCCATCCGCTGAGATGATGCGGCGCCACATggactgagc                                                                                                          >  2:345913/1‑136 (MQ=255)            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)                                                                                                                                           |                                                                                                                    CGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAT  >  NC_000913/1196082‑1196335

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: ‑