Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,196,232	C→T	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	96.9%	76.6 / NA	32	T370T (ACC→ACT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/14);  total (18/14)
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

TGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  NC_000913/1196094‑1196315                                                                                                                                           |                                                                                    tGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAg                                                                                 <  1:238931/143‑1 (MQ=255)    tGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                              <  1:91539/143‑1 (MQ=255)    gtcataCAGCGGGGACGAATGCGCCCTGTTTCAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCgg                                                                              <  1:703159/138‑1 (MQ=255)     gCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCCGCGATGCTGCTGCGCCATATGGGTTGGACTGAagcagc                                                                             >  1:274267/1‑143 (MQ=255)      cAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGCAATCTGCCGGCCAGGACATAGTAAATCCTGGCTCTATTATTCTCTCCGCCGAGATGATGCTCCGCCATATGCGTTGCACCGACACGGCt                                                                            >  2:509291/1‑143 (MQ=255)         tcttctgTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc                                                                         <  1:488571/137‑1 (MQ=255)              ggTGGCGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctgat                                                                    <  2:468501/143‑6 (MQ=38)                    gAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGCCAGGACAAAGTAGATCGTGGCTCTATCACTCTCTCCGCTGCGATGATGCTGCGCCATATGGGTTGGACTGAAACGGCTGACCTGATTGTTaa                                                              >  1:657632/1‑143 (MQ=37)                     aaTGCGCCCCGTTTGACGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGTTGCGCCATATGTGTTGGACTGAAGTGGCTGACCTGATTGTTaaa                                                             <  2:245811/143‑1 (MQ=38)                          tcccTGTTTGAATCCTCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                        <  2:328170/142‑1 (MQ=37)                          gcCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                        <  1:583481/143‑1 (MQ=37)                           cccTGTGTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATg                                                       <  2:208957/143‑1 (MQ=37)                            ccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATgg                                                      >  2:439132/1‑143 (MQ=255)                             cTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                     >  1:700218/1‑143 (MQ=37)                             cTGTGTGAGGCCACCCGCGGTACTGCGCCGAGATAGGCCGGTCCGGACGAAGTAAATCCTGGCTCTATGATGCGCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGa                                                     <  2:739672/143‑1 (MQ=21)                              tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGaa                                                    <  1:915714/143‑1 (MQ=37)                                tttGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAgg                                                  >  1:101786/1‑143 (MQ=37)                                 ttGAAGCCACCCACGGTACTGCGCCGAAATATGCCTGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGgc                                                 <  2:450800/143‑1 (MQ=37)                                   gAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCGATAAGGGATGGACTGAAGCGGCTGACCTGCATGTTAAAGATATGGAAGgcgc                                               >  2:1114890/1‑143 (MQ=37)                                     aGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCaa                                             >  1:808353/1‑143 (MQ=25)                                       ccacccacGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATGGTTACAGGTATGGAAGGCGCACTc                                           >  2:361318/1‑143 (MQ=25)                                       ccacccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTGAAGGTATGGAAGGCGCGATc                                           >  1:321713/1‑143 (MQ=25)                                         acccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGACGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGCGAAAGGGATGGAGGGCGCAATCaa                                         >  2:894686/1‑143 (MQ=25)                                                        cgggaaTGTGCCGGTCAGGGCAGAGTAAATCCTGGCTCTAGTATTGTCTCCGCTGAGATGATGCTGCGGCATATGGGTTGGACTGAAGCGGCTGACCTGAGTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAc                          <  2:186260/139‑1 (MQ=14)                                                         cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                         >  1:24818/1‑142 (MQ=18)                                                            aaTATGCCGTTCAGGACAAAGTAATTCCTCGCTCTATTATTCTCTCCGCTTAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTTATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                      <  2:290136/143‑1 (MQ=17)                                                                          gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt        >  2:582388/1‑143 (MQ=11)                                                                          gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt        >  1:955942/1‑143 (MQ=11)                                                                          gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt        >  1:227867/1‑143 (MQ=11)                                                                           aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGACCgggtg       >  2:763587/1‑139 (MQ=14)                                                                                gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatgga  >  1:1045309/1‑143 (MQ=11)                                                                                gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGCCCTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCCgatgga  >  1:904325/1‑143 (MQ=11)                                                                                                                                           |                                                                                    TGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  NC_000913/1196094‑1196315

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