Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,196,220	C→T	44.1%	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.1%	19.6 / 47.4	49	H366H (CAC→CAT)	icd	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base T (11/10);  total (25/24)
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.36e-01

GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCA  >  NC_000913/1196084‑1196359                                                                                                                                         |                                                                                                                                            gtaccgCCCCTGGGGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGt                                                                                                                                       <  2:385820/143‑1 (MQ=255)  tatcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGCTGCTGGGCCACATGGGtt                                                                                                                                      >  1:599219/1‑143 (MQ=255)    tcgCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTgg                                                                                                                                    >  2:623760/1‑143 (MQ=255)     cgCACCTGGTGCAAACATCGGTGACGAATGCCCCCTGTTTGAAGCCCCCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGa                                                                                                                                   <  2:846563/143‑1 (MQ=255)      gCCCCTGGTGCAAACATCGGGGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGAc                                                                                                                                  <  1:118586/143‑1 (MQ=255)        cccTGGTGCGAACATCGGGGACGGATGCGCCGGGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCg                                                                                                                                <  1:516368/143‑1 (MQ=255)          ntGGTGCAAACATCGGTGACGGATGCGCCCTGTTTGAAGCCACCCCCGGTACTGGGCCGAGATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGGTGAGATGGGGCTGAGCCCGACGGGTTGGagacgc                                                                                                                              >  2:629283/2‑138 (MQ=255)            ggTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACACAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGc                                                                                                                            <  2:301323/143‑1 (MQ=255)            ggTGCAAACATCGGTCACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGc                                                                                                                            <  2:774846/143‑1 (MQ=255)                  atcgtCGGTGACGAATGCGCCCTTTTTGAAGCTACCTACTGTATTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGa                                                                                                                      <  1:306596/139‑1 (MQ=38)                  aaCATCGGTGACGAATGCGCCCTGGTTGAAGCCACCCACGGTACTGCGCCGAAATAGGCCGGTCAGGACAAAGTAAGTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGa                                                                                                                      <  2:581643/143‑1 (MQ=255)                    cATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAcc                                                                                                                    <  2:291593/143‑2 (MQ=255)                      tCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAATTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTa                                                                                                                  <  1:495036/143‑1 (MQ=255)                         gTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                               <  1:696507/143‑1 (MQ=255)                         gTGACGAATGCGCCCTGCTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAAtt                                                                                                               <  1:377162/143‑1 (MQ=255)                               aaTGCGCCCTGTTTGAAGCCTGCCATGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTaaa                                                                                                         <  1:1038553/143‑1 (MQ=37)                               aaTGCGCCCTGTGTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                         <  2:458070/143‑1 (MQ=255)                               aaTGCGCCCGGGTTGAAGACACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                                                                                         <  1:22080/143‑1 (MQ=255)                                aTGCGCCCTGTCTGGAGCCACCCACGTTAGTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAg                                                                                                        <  2:658206/143‑1 (MQ=37)                                    gcCCTGTTTGAAGCCACCCACGGGCCTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTAt                                                                                                    <  1:29991/143‑1 (MQ=37)                                     cccTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCCCTATTATTCTCTCCGCTGAGATGATACTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGGTAAAGGTacg                                                                                                   >  2:327072/1‑141 (MQ=37)                                        tGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGGCACATGGGTTGGACCGAAGCGGCTGACTTTATTATTAAAGGTATGGag                                                                                                >  1:517968/1‑142 (MQ=255)                                              aaGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCa                                                                                          <  1:623760/143‑1 (MQ=255)                                                 ccacccacGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGACTGTTCAAGGTATGGAAGGCGCAATc                                                                                       >  2:624557/1‑143 (MQ=25)                                                         ggTACTGCGCCGGAATATGCCGGTCAGGACAAAGTACATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaatgccaa                                                                               <  2:945639/143‑7 (MQ=18)                                                              tGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGCAGCGGCATACCCGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCg                                                                          >  1:770375/1‑143 (MQ=18)                                                                 gcCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTaa                                                                       >  2:423038/1‑143 (MQ=18)                                                                 gcCGAAAGATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCGCCGCTGAGATGATGCGGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa                                                                       <  2:616843/143‑1 (MQ=255)                                                                  ccGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAAc                                                                      >  1:461060/1‑143 (MQ=255)                                                                   cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTCAAGGTATGGAAGGCGCAGTCAATGCCCAGACCGTAAct                                                                     >  1:19350/1‑142 (MQ=18)                                                                   cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCCTAAct                                                                     >  2:957954/1‑142 (MQ=18)                                                                   cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCCGACCCAACTGCTCAAGGTATGGACGGCGCAATCAACGCGAAAACCCTAAcc                                                                     >  2:856588/1‑143 (MQ=255)                                                                   cGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCTGCGACATATAGGTTGGACTGAAACGGCAGACCTGATTGTTAAAGGTATGAAAGGCGCAATCAATGCCAAGACCGTaaac                                                                     >  2:87809/1‑141 (MQ=18)                                                                     aaaTATGCCGGTCAGGACAAAGTAAATCCTGGCGCTATTATTCTCTCCGCTGAGATGATGCGGCGCCATACGGGTTGGACTAAAACGGCTGACCTGATAGTTACAGGTAAGAACGGCGCAATAAATGCCAAGACCGAAACTTa                                                                   >  2:285576/1‑143 (MQ=255)                                                                      aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGCTGCTGCGCCATATGGGTTGGACTGACGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAAcacat                                                                  >  2:54267/1‑139 (MQ=17)                                                                          tGCCGGTCAGGACAAAGTGAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGGCACCTGGGTTGGACCAAAGCGGCTGACTTAATTGTTAAAGGTAGTGAAGCAGCAACCAACGCAACAACCGTAAcccctcact                                                              >  2:1082064/1‑136 (MQ=255)                                                                           gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACtt                                                             >  1:335317/1‑143 (MQ=17)                                                                           gCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACtt                                                             >  2:44594/1‑143 (MQ=255)                                                                            ccGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTc                                                            <  1:44594/143‑1 (MQ=255)                                                                                    gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt                                                    >  1:932292/1‑143 (MQ=255)                                                                                    gACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCt                                                    >  2:981290/1‑143 (MQ=255)                                                                                       aaaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAATGCCAAGACCGTAACTTATGACTTCGAACGTCTgct                                                 >  2:912936/1‑141 (MQ=14)                                                                                                   ggCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAACGAAGCGGCGGACTGAGTTGGTAACGGTAGGGAGGGCGCAATCGACGCGAAAGCCGTAACCAATGACTACGGGAGTCTAGTGGAGGgcgttgc                                     >  1:389101/1‑139 (MQ=255)                                                                                                      tctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGTACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGAATTCGAGCGACTGATGGATGGCGCTAAAct                                  >  2:370592/1‑143 (MQ=255)                                                                                                         attattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAActgct                               <  2:461060/143‑1 (MQ=255)                                                                                                            atgctctcCGCTGAGATGATGCGGCGCCGCGTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGaa                            <  2:694439/140‑1 (MQ=255)                                                                                                              tttctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGGAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAAt                          <  2:131058/141‑1 (MQ=255)                                                                                                                                  tgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGTCGCTAAACTGCTGAAATGTCCAGAAAGTGGTGAcgcg      >  1:1003491/1‑143 (MQ=18)                                                                                                                                      ccACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgccggca  >  2:909999/1‑138 (MQ=255)                                                                                                                                         |                                                                                                                                            GTATCGCCCCTGGTGCAAACATCGGTGACGAATGCGCCCTGTTTGAAGCCACCCACGGTACTGCGCCGAAATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCA  >  NC_000913/1196084‑1196359

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