Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	1,757,058	G→C	16.7%	G454A (GGT→GCT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	1,757,058	0	G	C	16.7%	54.0 / 7.0	29	G454A (GGT→GCT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (15/9);  new base C (3/2);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

CTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑TTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGA  >  NC_000913_3_hsa_tpiA/1756914‑1757205                                                                                                                                                 |                                                                                                                                                     cTCATCAGTTGGGACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCa                                                                                                                                                  <  1:368295/149‑1 (MQ=255)  tcatcaGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCAc                                                                                                                                                 >  2:339896/1‑149 (MQ=255)           ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGa                                                                                                                                        <  1:540653/149‑1 (MQ=255)                       aGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACATTGTAGTTATGGGTTCTGGTGCACTGGTACCCAGCGGCActgcta                                                                                                                            >  2:161488/1‑149 (MQ=255)                                      gCTTGTTAAAGAGATCACTTCAACTGATGATTGCTACAGTCTGAGTAAAGAACTGGATCTGCAGAGCTATCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                                                                                             <  2:608611/149‑1 (MQ=255)                                      gCTTGTTAAACAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTc                                                                                                             <  1:276989/149‑1 (MQ=255)                                                   aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCTGAAGGTTTGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAata                                                                                                <  1:454609/149‑1 (MQ=255)                                                        ttCTACTGATGATTTCTACCGTCTGGGTAAAGCACTGGCTCTGAAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAAAACCGCATCTGTTAGCGTCCTGTAATATTGCt                                                                                           >  1:74782/1‑149 (MQ=255)                                                        ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGTCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCt                                                                                           >  1:111108/1‑149 (MQ=255)                                                        ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCACCTGTTCACGTCCTGTAATATTGCt                                                                                           >  1:273678/1‑149 (MQ=255)                                                              tgatgaTTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTtgtg                                                                                     >  2:111636/1‑149 (MQ=255)                                                              tgatgaTTTCTACCAACTGGGTAAAGAACTAGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCGTCTGCTCCCGTCCTGTAATATTGCatttata                                                                                     >  2:111267/1‑146 (MQ=255)                                                                atgatTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACGCCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaa                                                                                   >  1:422896/1‑149 (MQ=255)                                                                     ttCTACCGTCTGGGTAAAGAACTGGCTCTGCAGATCGTGCTTGCCCCAACTGCTGGCGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGCTCACGTCCTGTAATATTGCTTTTGTGaattaat                                                                              <  1:544929/149‑1 (MQ=255)                                                                       cTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAttt                                                                            >  1:447391/1‑149 (MQ=255)                                                                          ccGTCCGGGTAAAGAACTGGCTCTGCAGAGTTTTCTGGCACACAAAGGTGACGTAGTAGTTATGGTTTCTGGTGCACGGGTACCGAGCGGCACTACTAACACCGCATCTGTACACGTCCTGTAATATTGCTGTTGTGTATTAATTTGtg                                                                         >  2:468939/1‑148 (MQ=255)                                                                            gTCTGGGTAAAGAACTGGCTCTGCAGAGAGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCAATGGTACGGAGGGGAACTACTAACACGGCATCTGTTCACATCCTTGAATATTGCTTTAGTGAATTAATTAGtcta                                                                       >  2:67157/1‑149 (MQ=255)                                                                                  gTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGAAATATTGCTTTTGTGAATTAATTTGTATATCGAAg                                                                 >  1:379003/1‑149 (MQ=255)                                                                                   tAAAGAACTGGGTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGGATTAATTTGTATATCGAAgc                                                                >  1:377555/1‑149 (MQ=255)                                                                                                 tGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAACTTGCATGTCGAATCGCCCTCATTGgcgc                                                  >  2:435081/1‑149 (MQ=255)                                                                                                   cAGGGCGGTCTGGCACACAAAGGGGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑tt                                               <  2:550246/149‑1 (MQ=255)                                                                                                    agagCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCATGTAATAATACGTTTGTGAATTAATTTTTATATCGAAGCGCCCTGTTGTGCGC‑ttt                                              >  1:596057/1‑149 (MQ=255)                                                                                                           gttGGCACGAAAGGGCGACGTTGTACTTATGGTTCGTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑TTTTtttatt                                       <  2:17078/147‑1 (MQ=255)                                                                                                                              gTTGTAGTTATGGTTTCTGCTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGGATATCTAAGCGCACTGATGGGCGC‑TTTTTTTATTTAATCGATAACCAGAAGCa                    >  1:144916/1‑149 (MQ=255)                                                                                                                               ttGTAGTTATGGATTCTGGTGCACTGGTACCGAGCGGCACGACTAACAACGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGCTTTTTTTTATTTGATCTATAACCACAACCa                    >  2:279820/1‑149 (MQ=255)                                                                                                                                tGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑TTTTTTTATTTAATCGATAACCAGAAGCAAt                  <  1:358645/149‑1 (MQ=255)                                                                                                                                     ttATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTCATTTGTATATCGAAACTCCCTGAGGGGCGC‑TTTTTTTATTTAATCGATAACCAGACGCaccaacaa             >  2:118008/1‑142 (MQ=255)                                                                                                                                          gTTTCTGGTTCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑TTTTTTTATTTAATCGATAACCAGAAGCAATAAAaaatcaa        <  1:187566/149‑1 (MQ=255)                                                                                                                                             tCTGCTGGATTGGTACTGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑GTTTTTTATCTAAGCACTTGCTTGAAGCAATAAAaaatcaaatc     <  2:442606/149‑1 (MQ=255)                                                                                                                                                ggTGCAGTGGGACCGAGCGTCATTTCTAACACCGCTTCTTTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATTTCGATGCGCCCTGATTGGCGC‑TTTTTTTATTTAATCGTTAACCAGAGGCAATAAAAAATCATATCGGa  <  2:447767/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAATTTGTATATCGAAGCGCCCTGATGGGCGC‑TTTTTTTATTTAATCGATAACCAGAAGCAATAAAAAATCAAATCGGA  >  NC_000913_3_hsa_tpiA/1756914‑1757205

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