Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	264,098	G→T	L50L (CTC→CTA)	ykfH ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	264,098	0	G	T	97.6%	114.1 / ‑3.4	42	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (23/18);  total (24/18)
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.12e-01

GGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_hsa_tpiA/263953‑264242                                                                                                                                                   |                                                                                                                                                 gggCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAgg                                                                                                                                               >  1:469776/1‑149 (MQ=255)      cccAGGCGCAGTGTAATCTCCGTCTCCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTcc                                                                                                                                          <  2:107870/149‑1 (MQ=255)       ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCg                                                                                                                                         >  2:163987/1‑149 (MQ=255)         aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCgcc                                                                                                                                       >  2:149277/1‑147 (MQ=255)         aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                       >  2:149109/1‑149 (MQ=255)                 tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCTACTATTGTCAGTTCAGGGTGGTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCg                                                                                                                               >  2:207725/1‑149 (MQ=255)                   tAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCGACTATTGTCCGTTCCGGGTGTTGTTATTCTGGTAGCTCTAGCGGTGGTTGCGGGCGTGA‑CGGTGTGCGGCGTCTACCGCTAGGACGCCGGCGATAGCCGGa                                                                                                                             >  2:53407/1‑149 (MQ=255)                     aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAAc                                                                                                                           <  1:256079/149‑1 (MQ=255)                        tcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTc                                                                                                                        <  2:231721/149‑1 (MQ=255)                          ccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGGTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTCGGACTCCGGCCATATCCGGAACGTccc                                                                                                                      >  2:377171/1‑149 (MQ=255)                            gCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGGTAGGACTCCGGAGATATCGGGAACGTCCCTg                                                                                                                    >  1:198149/1‑149 (MQ=255)                                   gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                             >  1:41133/1‑149 (MQ=255)                                      ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                          >  1:104029/1‑149 (MQ=255)                                          gTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACa                                                                                                      >  1:257598/1‑149 (MQ=255)                                               gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCGGGGTAAGGACGGACGc                                                                                                 >  1:204925/1‑149 (MQ=255)                                                      aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTg                                                                                          >  1:407285/1‑149 (MQ=255)                                                           aGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCCCt                                                                                     >  1:276127/1‑149 (MQ=255)                                                              tctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGt                                                                                  <  1:53407/149‑1 (MQ=255)                                                               ctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                 <  2:462199/149‑1 (MQ=255)                                                                tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGGATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATAACAGACGCTACCGTGCCACTGGTAt                                                                                >  2:471135/1‑149 (MQ=255)                                                                   aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                 <  2:340258/113‑1 (MQ=255)                                                                   aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCg                                                                                                                 >  1:340258/1‑113 (MQ=255)                                                                      aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                          <  1:446331/149‑1 (MQ=255)                                                                          tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                      <  2:469776/149‑1 (MQ=255)                                                                          tCAGTTCAGGGTGATGCTCAAGAGGGAGGGATAGGGCAGCTTAGTGGCTTGGCCGGCTTATCGCG‑ATTTTGCGTGGTCTCCTTCGATATCCGGAACGTCCCTGCGGGTGTAATAACAGACGCTACCGTGCCACTCGTATTTGCCGGttc                                                                      >  2:581615/1‑147 (MQ=255)                                                                          tCAGTGCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                      >  2:337453/1‑149 (MQ=255)                                                                              ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                  <  1:362577/149‑1 (MQ=255)                                                                              ttCAGGGTGATGCTCATCAGGCAGGGATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACGGCCGTTACCGTGCCACGGGTATTTTCCGGTGCAGtg                                                                  >  1:36450/1‑148 (MQ=255)                                                                               tCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAg                                                                 <  1:207725/149‑1 (MQ=255)                                                                                    gcgATGCTCATCAGGCAGGTATAGGGCCCATTGCTGTCCCGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                            <  1:616/147‑1 (MQ=255)                                                                                    gTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGaaa                                                            <  1:163987/149‑1 (MQ=255)                                                                                          cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTCTACTGACAGCAGATACCGTACCACGGGTCTTTGCCGGTGGAGTAGCGAGAGATTCg                                                      >  2:517712/1‑149 (MQ=255)                                                                                           tcatcaGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATCTCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACCGGTCTTTGCCTGTGCCGTAGCGAAAGATCCgc                                                     >  2:139651/1‑148 (MQ=21)                                                                                                         taGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTg                                       <  1:33686/149‑1 (MQ=255)                                                                                                          cgggCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTgg                                      <  1:218682/148‑1 (MQ=255)                                                                                                           gggCCCATTTCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGc                                     <  2:139176/149‑1 (MQ=255)                                                                                                            ggCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                    <  2:407285/149‑1 (MQ=255)                                                                                                                cATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTAACGATGTAAGGACAGACCCTACACTGCCATCGCTAATTTTCGTTGCATTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGtat                                >  2:280948/1‑149 (MQ=17)                                                                                                                     cGGCCCTGG‑CGGCGTCCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCGGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGt                           <  1:392227/149‑1 (MQ=255)                                                                                                                                       ggCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc          >  1:159610/1‑149 (MQ=255)                                                                                                                                        gCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCGTTGCCCTGCGTTTACt         >  1:154425/1‑149 (MQ=255)                                                                                                                                               ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAtt  <  2:41133/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGG‑CGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATT  >  NC_000913_3_hsa_tpiA/263953‑264242

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