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	97.2%	78.4 / NA	36	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/15);  minor base G (0/1);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  NC_000913/1196080‑1196315                                                                                                                                             |                                                                                                atcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCatat                                                                                               >  1:13512/1‑143 (MQ=255)     gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTCTTGGCTGCGCTGAGATGATGCTGCGCCATGTGGGt                                                                                           >  1:608381/1‑143 (MQ=255)     gtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTT‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACCAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGATATGATGCTGCGCCATATGGGtt                                                                                          >  1:659118/1‑143 (MQ=255)         cgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGa                                                                                       <  1:1127976/143‑1 (MQ=255)          gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTATTGCGCCGAAATATGCCTGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAc                                                                                      <  1:1148539/143‑1 (MQ=255)          gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGTCGAAATATGCCGGTCAGGAAAAAGTAAAGCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAc                                                                                      <  1:148032/143‑1 (MQ=255)          gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATGTGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTgggc                                                                                      >  1:1407171/1‑141 (MQ=255)            cccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTAGCTCTATTATTCTCTCCGCCGAGATGATGCTGCGCCATATGGGTTGGActg                                                                                    >  1:1307906/1‑141 (MQ=255)              cTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATCTGCCGGTCAGGACAAAGTAAATCCTGCCTCTATTATTCTCTCCGCTGAGATGATGCTCCGTCGTACGGGTTGGActgca                                                                                  >  2:1307499/1‑139 (MQ=255)                      aaCATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                                          <  1:932284/143‑1 (MQ=255)                                   aaTGTGCCCGGTTGGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                             <  2:563414/143‑1 (MQ=37)                                    aTGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                            <  2:13512/143‑1 (MQ=38)                                       cgcCCTGTTTGAAGCCACCCACGGTACTGCGCCGCAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTCTTACAGGTa                                                         >  2:424902/1‑143 (MQ=37)                                           cTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCGGCTGGCCTGCTTGTTAAGGGTCTggc                                                     >  1:343878/1‑142 (MQ=37)                                              tttGAAGCCACCCACGGTACTGCGCAGAAATATGCCGGTCAGGACAAAGTAAATCCTAGCTCTATTATTCTCTAAGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAGGGAAgg                                                  >  1:175335/1‑143 (MQ=25)                                               ttGAATCCACCCGCGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                                 <  2:1377406/143‑1 (MQ=25)                                               ttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGACTGTTAAAGGTATGGAAGgc                                                 >  2:1156093/1‑143 (MQ=37)                                                tGACGCCACCCACGGTACTGCGCCGAAATGTGCCGGTCAGGACAAAGTAAATCCTGGCTCTATGATTCTCTCCGCTGAGATGATGCTGCGCCAGATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgcg                                                <  2:1247323/143‑1 (MQ=25)                                                    gCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAAt                                            <  2:1388780/143‑1 (MQ=25)                                                          cgctgTAGTGCGCCGAAATATGCAGGTGAAGACAAAGTAAATCCTGTCTTTATTATTCTCTCCGCTGAGATTATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc                                      <  2:823643/139‑1 (MQ=17)                                                          cacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTAACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc                                      <  1:1441143/143‑1 (MQ=21)                                                          cacGGTACTGCGCCGAAAGATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCGGCGCCATAAGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc                                      <  2:1122114/143‑1 (MQ=21)                                                                        taaaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCGCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCGATGCCAAGACCGTAActt                        >  2:1168998/2‑141 (MQ=17)                                                                        gAAATATGCCGGTCAGGTCAAGGTAAATGCTGGCTCGATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAGTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAActt                        <  2:1374002/143‑3 (MQ=17)                                                                            tatGCCGGTCAGGACAAAGTAAATCCTGGCTTTATTATGCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGa                    <  2:1377510/143‑1 (MQ=17)                                                                              tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACt                  >  1:713073/1‑143 (MQ=17)                                                                                 cGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGCCCTGATTGTTAAAGGTATGGACGGCGCAATCAATGCCAAGACCGTAACTTATGACTGCg               >  2:106156/1‑143 (MQ=17)                                                                                   gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGGTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCTCAATCAATGCCAAGAGCGTAACTTATGACTTCGaa             >  1:179124/1‑142 (MQ=17)                                                                                      aGGACAGAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGt          <  2:799097/143‑1 (MQ=14)                                                                                        gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt        >  1:1242282/1‑143 (MQ=11)                                                                                         aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg       >  2:897123/1‑143 (MQ=11)                                                                                         aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTg       >  1:285364/1‑143 (MQ=14)                                                                                             aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgg   >  1:351221/1‑143 (MQ=11)                                                                                              gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGCACGTCTgatgga  >  2:745256/1‑143 (MQ=11)                                                                                              gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACAGTAACTTATGACTTCGAACGTCTgctgga  >  1:186621/1‑143 (MQ=11)                                                                                                                        cGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                                                         <  2:352366/46‑1 (MQ=11)                                                                                                                        cGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                                                         >  1:352366/1‑46 (MQ=11)                                                                                                                                             |                                                                                                ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  NC_000913/1196080‑1196315

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