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	96.9%	78.7 / NA	32	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (16/15);  total (16/16)
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

CTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_hsa_tpiA/263956‑264240                                                                                                                                                |                                                                                                                                               cTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTCAATTCAGGGTGATGCTCATCAGGCAGGAATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGGCt                                                                                                                                          >  1:262170/1‑149 (MQ=255)  tCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACGGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGCGGTTCG‑GCGTAGACCGCTAGGACTc                                                                                                                                         >  2:276823/1‑149 (MQ=38)    ccAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTGCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCg                                                                                                                                       <  2:327757/149‑1 (MQ=255)                      cccGCGGGAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGGGATGCTCATCAGGCAGGTATAGGGCGCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                                     <  1:534731/149‑1 (MQ=255)                                   ccTGTTAGGGTTTGTCCTTACAGCTCCCCACTTTTGTTGGTTGAGGTTGATTCGCAGCAGGCAGGTAGAGGGCCCATTGCTGTCCTGGCGTCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                        <  1:262135/149‑1 (MQ=255)                                    cTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                       <  1:153569/149‑1 (MQ=255)                                          gggTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGCACGTCCCTGCCGGCGTAATGCCAGAc                                                                                                 >  2:320779/1‑149 (MQ=255)                                            gTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGc                                                                                               >  2:509559/1‑149 (MQ=255)                                                  cATTACAGTTCTCCACTATTCTCAGTTCAGGGAGATGCTCATCAGTCCGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGCCTCCGCCGATGTTCGCAACGTCCCTGCCGGCGGCATGACCGCCGCTACCGt                                                                                         >  2:66512/1‑149 (MQ=255)                                                               cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCGCTGGTAttt                                                                            >  1:200049/1‑149 (MQ=255)                                                               cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                            >  1:165195/1‑149 (MQ=255)                                                               cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGAGCCACTGGTAttt                                                                            >  1:555559/1‑149 (MQ=255)                                                                 cTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGc                                                                          >  2:359706/1‑149 (MQ=255)                                                                   aTTTTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                        <  1:377255/149‑1 (MQ=255)                                                                           ttCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGAGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTCATGACAGACGCTACCGTGCCACTGGGATTTGCCGGTGGAGTa                                                                >  1:138501/1‑149 (MQ=25)                                                                           tgcAGGGTGATGCTTATCAGGTAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTa                                                                <  2:326866/147‑1 (MQ=255)                                                                             cAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGc                                                              <  1:235881/149‑1 (MQ=255)                                                                             cAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGGCCTGGCGGCGTTCG‑GCGTATACCGCTAGGCCTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGc                                                              <  1:184728/149‑1 (MQ=255)                                                                                       cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                    <  1:183869/149‑1 (MQ=255)                                                                                        tcctcaGGGTGGTATCAGGACCATTGCGGGCCTGGCGG‑GCTCGCGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCgg                                                   <  1:533375/146‑1 (MQ=17)                                                                                          atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTAA‑CCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGCAAGAATCGGGa                                                 >  1:509981/1‑149 (MQ=21)                                                                                                  ggTAGAGGGCTCTTTGGGGTCCTGGCGGCGGTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGAt                                         <  2:54007/149‑1 (MQ=18)                                                                                                         ggCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                  >  1:466134/1‑149 (MQ=255)                                                                                                          ggccATTGGGGTCCTGGCGGGGTGCG‑GCGGATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCgg                                 <  2:357101/147‑1 (MQ=17)                                                                                                            ccATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGta                               >  1:51594/1‑149 (MQ=255)                                                                                                                tGCGGTCCTGGCGGGGTGCG‑GGGGATACCGCTAGGACTCCGGCGATAGCCGGAACGGCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                           <  1:11821/149‑1 (MQ=17)                                                                                                                tGCGGGCCTGGCGGCGTTCG‑GTGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTACGTGCCGGGGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                           <  2:155608/149‑1 (MQ=17)                                                                                                                          gcggcgTTCG‑GCGTATACCCCTAGGACTCCGGCGATATCCGGAACGGCCCTGCGGGGGTAATGAGGGACGGTACCGTGCCACTGGGATTTGCCGGTGCCGTAGCGAAAGATTCCGGGCTCAGCATGCTCGCGGCCTCACGGCATTGccg                 >  2:549076/1‑148 (MQ=255)                                                                                                                               gTTCG‑GCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGGGCAGTAGCGAAAGATTCGGGACTCAAGATGCTGGCGGTATCTCGTCATGGCCCTGCGt            >  1:47624/1‑149 (MQ=17)                                                                                                                                     gCGAATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGTTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACt       <  1:111437/149‑1 (MQ=255)                                                                                                                                          taCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGGGAGACGCTACCTGGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAGTGCCCTGCGTTTACGGATaa  >  1:80963/1‑149 (MQ=17)                                                                                                                                          taCCGCTAGGACCCCGGCGATATCCGGAACGCCCCTGACGGTGTAATGACAGACGCTACCCTGCCACCGGTATTTGCCCGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa  <  1:276909/149‑1 (MQ=17)                                                                                                                                                |                                                                                                                                               CTCCCAGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCG‑GCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_hsa_tpiA/263956‑264240

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