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	95.3%	120.7 / ‑5.0	43	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (21/20);  total (22/21)
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

GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263954‑264241                                                                                                                                                  |                                                                                                                                                ggCTCCCAGGCGCAGTGTAATCGCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGa                                                                                                                                             <  1:72565/149‑1 (MQ=255)     cccAGGGGAAGTGTAATCGCCCTCTGCACGCTCCGGGTAGTTTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGTAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGTTGTATAC‑CGCTAGGACTcc                                                                                                                                         <  2:66720/149‑1 (MQ=255)        agggcctgtgtATTATCCCGCTGTCGCCCCTTGGTAGGGCTGCTCATTACAGGTCGCCACTATTGTCAGTTCAGGGTGATGCTCATCCGGCATGTATAGGTCCCATTGCGGTCCTGGTGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGc                                                                                                                                      <  2:212685/142‑1 (MQ=255)        aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGc                                                                                                                                      <  2:472591/149‑1 (MQ=255)             cAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGatat                                                                                                                                 >  1:418178/1‑149 (MQ=255)                tgtAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCg                                                                                                                              <  1:81571/149‑1 (MQ=255)                    aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACCGTTCTCCCCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGCCTCCGGCGCTATCCGGacg                                                                                                                          >  2:235117/1‑147 (MQ=38)                     tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGAGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGACCg                                                                                                                         >  2:526554/1‑149 (MQ=255)                          cGCTGCAGGCCCCTGGGAGGGTTTCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCt                                                                                                                    <  2:147904/149‑1 (MQ=255)                                   ccccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCTTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgggtg                                                                                                           >  1:418994/1‑148 (MQ=255)                                    cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                          >  2:468578/1‑149 (MQ=255)                                     ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                         >  1:807300/1‑149 (MQ=255)                                      cTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGTGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                        >  2:683809/1‑149 (MQ=255)                                              gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCg                                                                                                                              <  2:290872/119‑1 (MQ=255)                                              gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCg                                                                                                                              >  1:290872/1‑119 (MQ=255)                                              gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGc                                                                                                >  1:286866/1‑149 (MQ=255)                                               ttGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCt                                                                                               >  2:436408/1‑149 (MQ=255)                                                         ccgTTCCCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCAc                                                                                     <  2:471749/147‑1 (MQ=255)                                                         cAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCTGTGTAATGACAGACGCTACCGTGCCAc                                                                                     <  2:471730/149‑1 (MQ=255)                                                              cttcACTATTGTCAGTGCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                <  1:639502/146‑1 (MQ=255)                                                              ctcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTa                                                                                >  2:327000/1‑149 (MQ=255)                                                               tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACAGTGCCACTGGTAt                                                                               >  2:348015/1‑149 (MQ=255)                                                                 cACTATTGTCAGTTCAGGGTGATGCTCATCGGGCAGGTATCCGACATGCGGCGAACATAGCGTCGTGCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACGGTGCCACTGGTAttt                                                                             >  1:408815/1‑149 (MQ=18)                                                                 cACTAGTGTCAGGTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                             <  1:436308/149‑1 (MQ=255)                                                                  aCTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTg                                                                            >  1:177939/1‑149 (MQ=255)                                                                     aTTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                         <  2:723555/149‑1 (MQ=255)                                                                      ttGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCCGACCTCGCACGGTGGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCgg                                                                        >  1:193388/1‑149 (MQ=21)                                                                         tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGCTCCTTACTGCGTTCGGCGTATAC‑CGCGAGGACTCCGGCGAGATCCGGACCGTCCGTGCCGGTGTGACGACAGACGCTACCGTGCCCCAGGCAGTTGCCGGTGc                                                                     >  2:4551/1‑149 (MQ=255)                                                                                             tcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGAc                                                 <  1:590507/149‑1 (MQ=255)                                                                                               aggcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTc                                               <  2:286866/149‑1 (MQ=255)                                                                                                 gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGc                                                                                                                                      <  2:802445/60‑1 (MQ=255)                                                                                                 gcaggTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGc                                                                                                                                      >  1:802445/1‑60 (MQ=255)                                                                                                  caggTATAGGGCCCATTGCGGTCCTGGAGGTTTTCGGCGTACACACGTG‑TGATTCTGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAgg                                            <  1:609750/149‑1 (MQ=17)                                                                                                     gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGCAAGATTCGGGCATCAgcagg                                         >  2:747756/1‑145 (MQ=18)                                                                                                      tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGTACTCAGGATGc                                        <  1:525862/149‑1 (MQ=255)                                                                                                      tataGGGCCAATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATACGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                        >  1:73356/1‑149 (MQ=255)                                                                                                                   gCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCg                           <  2:73356/149‑1 (MQ=255)                                                                                                                     ggTCCTGGCGGCGTTCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCGGGCGGTATATCGTc                         <  1:747690/149‑1 (MQ=255)                                                                                                                             tatcggtCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                 <  2:418178/143‑1 (MQ=255)                                                                                                                                  ttCGGCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCCGGCGGTATATCGTCATTGCCCTGCgct            >  2:439506/1‑147 (MQ=255)                                                                                                                                     ggCGTATAC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTAGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc         <  2:177939/149‑1 (MQ=255)                                                                                                                                          ataC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATa    >  1:703199/1‑149 (MQ=255)                                                                                                                                            aC‑CGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAt  >  1:188376/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                GGCTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATAC‑CGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAT  >  NC_000913_3_hsa_tpiA/263954‑264241

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