Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,220	C→T	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	96.9%	73.3 / NA	33	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (19/12);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATG  >  NC_000913/1196081‑1196313                                                                                                                                            |                                                                                              tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATg                                                                                            >  2:378544/1‑143 (MQ=255)      atcgCCCCTGGTGCAAACATCGGGGACGAATGCGCCCTGTTTGCAGCCACCCATGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTATATCCTGGCTCTATTGTTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTg                                                                                       <  2:623047/143‑1 (MQ=255)          ccccTGGTGCAAACATCGGTGACGAATGCGCGCGGTGTGGAGCCACCCACGGTACTGCGCCGAACTATGCCGGTAAGGACAAGGTAAATCCTGGCCCTATTATTCTCTCCGTTGAGGTGGTGCTGCGCCATATGGGTTGGAct                                                                                   >  1:392564/1‑142 (MQ=255)                    aaaCATCGGTGACGAATGCGCCGTGTTTGCAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg                                                                         <  2:345061/143‑1 (MQ=255)                    aaaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg                                                                         >  2:396419/1‑143 (MQ=255)                      aCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGCTGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                                                       <  1:807658/143‑1 (MQ=255)                        aTCGGTGACGAATGCGCCCTGTTTGAGGGCACTCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAAATCGTGGCTGTATTATGCGCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcct                                                                     <  2:827383/143‑3 (MQ=37)                          cGGTGCCGATTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctga                                                                   <  2:668260/143‑5 (MQ=38)                            gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctgatt                                                                 >  2:767153/1‑137 (MQ=255)                               aCGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGCGATGATGCTCCGCCATCTGGGTTGGACTGAAGAGGCTGACCTGCTTGtt                                                              >  2:588963/1‑143 (MQ=255)                                  aaTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                           <  1:163102/143‑1 (MQ=255)                                   aTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                          <  1:468766/143‑1 (MQ=38)                                          cTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGTCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                   <  1:36981/143‑1 (MQ=37)                                             tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCGGATCGCTAAAGGTATGGAAgg                                                >  2:304760/1‑143 (MQ=37)                                             tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATAATACTCCTCGTATTGTGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                                >  1:673394/1‑143 (MQ=25)                                               tGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAACTGAGTGTTAAAGGTGAGGAAGgcg                                              >  2:273423/1‑143 (MQ=37)                                                  ggtcagccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCaa                                           <  1:398750/137‑1 (MQ=25)                                                    ccacccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCCATc                                         >  1:70080/1‑143 (MQ=25)                                                    ccacccacGGTACTGCGCCGAAATATGACGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATc                                         >  1:503309/1‑143 (MQ=25)                                                          acGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgcc                                   >  1:413863/1‑142 (MQ=21)                                                            ggTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCCCCATATGGGTTGGACTGAAGCGTCTGTCCTGATTATTAAGGTTATTGAATGCGCAATaaatgaaaa                                 >  1:302409/1‑137 (MQ=25)                                                             gTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatgccaag                                <  1:378544/143‑8 (MQ=21)                                                             gTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCCGATTGGTAAAGGTATGGAAGGCGCAATCAATGCGaag                                >  2:40530/1‑142 (MQ=25)                                                                   cgcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTa                          >  1:531959/1‑143 (MQ=18)                                                                     ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                        >  1:182844/1‑143 (MQ=17)                                                                      cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                       >  2:726441/1‑142 (MQ=18)                                                                               cgcttCAGGACCAAGTAAATCCTGTCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc              <  2:538996/139‑1 (MQ=14)                                                                               ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTc              >  2:738042/1‑143 (MQ=17)                                                                                   tCAGGACAAAGTAAATCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac          <  1:316572/143‑3 (MQ=14)                                                                                       gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGGTTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGtat      >  1:767488/1‑141 (MQ=14)                                                                                       gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCACAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt      >  1:724039/1‑143 (MQ=14)                                                                                         cAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTCCGAACGTCTga    >  2:818863/1‑143 (MQ=14)                                                                                           aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGCCGCCACATGGGTTGGACTGACGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTCACTTATGACTCTGAACGTctgctg  >  1:392264/1‑143 (MQ=14)                                                                                                                                            |                                                                                              TCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATG  >  NC_000913/1196081‑1196313

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