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.1%	83.4 / NA	35	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (20/14);  total (21/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

ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  NC_000913/1196080‑1196315                                                                                                                                               |                                                                                                atctgtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCatat                                                                                               <  1:650497/143‑1 (MQ=255)  tcggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTCGCTCTATTATTCTCCCCGCCGAGATGATGCTGCGCCATATg                                                                                              >  1:769089/1‑143 (MQ=255)    ggtatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGCCATCGCGCGACCGT‑CTCCTGCGAAATATGCCTGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATggg                                                                                            <  1:205677/143‑1 (MQ=255)         cgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGa                                                                                       <  1:869603/143‑1 (MQ=255)          gCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCCGAGATGATGCTGCGCCACATGGGTGGGAc                                                                                      >  2:200795/1‑143 (MQ=255)           ccccTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAct                                                                                     <  2:836620/143‑2 (MQ=255)               tGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAg                                                                                 >  1:460804/1‑143 (MQ=255)                      aaCGTCGGTGACGAATGCGCGCTGTTTT‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGAGACATATAGTTTGGACTGAAGCGGCTGa                                                                          <  2:550660/143‑1 (MQ=255)                           cGGGGACGAATGCCCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGGAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcctga                                                                     <  2:14507/143‑5 (MQ=38)                                  gAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaa                                                              <  1:166104/143‑1 (MQ=38)                                   aaTGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                             <  2:666444/143‑1 (MQ=255)                                    aTGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                            <  1:532838/143‑1 (MQ=38)                                     tGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAgg                                                           >  2:511700/1‑143 (MQ=255)                                        gcCCTGTGTG‑AAGCCACCCACGGTACTGCGCCGAAATACGCCGG‑CAGGCAAATGGGAAACCCTGGCTCTAGTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                        <  1:855776/143‑1 (MQ=18)                                            tGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACACAGTACACCCCTTTATTTTTATTCTCTCCGCTGAGGTGATGCTGCGCCATATGGGTTGGACTGACGCGGCTGACCTGATTGTTCACGGTATGGaa                                                    >  1:856222/1‑143 (MQ=16)                                            tGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGCAGCGGCTGACCTGATTGTTAAAGGGATGGaa                                                    >  1:594511/1‑143 (MQ=37)                                             gTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAg                                                   >  1:310542/1‑143 (MQ=37)                                                    aGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCaa                                             <  2:698613/143‑1 (MQ=255)                                                      ccacccacGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAACTCCTGGCTCTATTATTCTCTCCCCTGAGATGATGCTGCGCCATATGAGTTGGACTGAAGCGGCTGACCTAATTGTTAAAGGTATGGAAGGCGCAATc                                           >  1:354571/1‑143 (MQ=25)                                                      ccacccacGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATc                                           >  1:479938/1‑143 (MQ=25)                                                        acccacGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaa                                         <  1:217141/143‑1 (MQ=25)                                                         cccacGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaat                                        >  2:11716/1‑142 (MQ=255)                                                                tACTGCGCCGAGAGATGCCGGGCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgccaaga                                 <  1:227009/143‑6 (MQ=18)                                                                tACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTTAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAGTCAATGCCaaga                                 >  2:687155/1‑141 (MQ=21)                                                                tACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGACTGTTAAAGGTATGGAAGGCGCAATCAAtgccaaga                                 >  2:71997/1‑138 (MQ=21)                                                                   tGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                              >  1:82858/1‑143 (MQ=18)                                                                        cGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAct                         >  1:233744/1‑142 (MQ=18)                                                                           actctGCCGGTCAGG‑GATAAGTAAATCCTGGCTCTATTGTGCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTAt                      <  2:825624/139‑1 (MQ=11)                                                                                 ccGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCTGCTGACCTGATTGTTAAAGGTATGGAGGGCGCAAGCACTGCCAAGACCGTAACTCATGACTTc                >  2:902970/1‑143 (MQ=17)                                                                                  cGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATCTGGGTTGGACTGAAGCGGCTGACCTGATTGTTAGAGGCATGGACGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCg               >  2:483675/1‑143 (MQ=17)                                                                                   ggTCAGG‑ATAAAGTGAATCCTTGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATAGGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGa              <  2:871780/143‑1 (MQ=11)                                                                                     tCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGaac            >  1:570497/1‑141 (MQ=14)                                                                                         g‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTCAAGGTATGGAAGGCGCACTCAATGCCAAGACCGTAACTTATGCCTTCGAACGTCt        >  1:357583/1‑143 (MQ=11)                                                                                         g‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCt        >  1:229535/1‑143 (MQ=14)                                                                                              aaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgatg    >  1:60853/1‑143 (MQ=11)                                                                                                gTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGCTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAGCGTCTgatgga  >  2:134127/1‑143 (MQ=14)                                                                                                                                               |                                                                                                ATCGGTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTG‑AAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGG‑ACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGA  >  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: ‑