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	100.0%	132.8 / NA	45	L50L (CTC→CTA)	ykfH	uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/23);  total (22/23)

AGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_hsa_tpiA/263961‑264240                                                                                                                                          |                                                                                                                                               aGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                     >  1:136354/1‑149 (MQ=255)   gcgcAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGa                                                                                                                                   >  1:634068/1‑149 (MQ=255)            aaTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGaa                                                                                                                          <  1:448305/149‑1 (MQ=255)                tcCCGCTGCAGGCCCCTGGTAGGGTTCCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTc                                                                                                                      <  1:566449/149‑1 (MQ=255)                 ctcGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATCTCCGGAACGTcc                                                                                                                     >  2:463841/3‑149 (MQ=255)                        cAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTCCCgg                                                                                                              >  2:120126/1‑149 (MQ=255)                           gCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCCATATCCGGGACGTCCCTGCCGgtgt                                                                                                           >  1:388377/1‑149 (MQ=255)                            ccccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTa                                                                                                          <  2:634068/149‑1 (MQ=255)                             cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                         >  1:186501/1‑149 (MQ=255)                              ccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAAt                                                                                                        <  2:136354/149‑1 (MQ=255)                               cTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCATGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                       >  1:94132/1‑149 (MQ=255)                               cTGGTAGGGTTGCTCATTACAGTTCTCCACTATCGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATg                                                                                                       >  1:824327/1‑149 (MQ=255)                                   tAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                     <  2:354657/115‑1 (MQ=255)                                   tAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGc                                                                                                                                     >  1:354657/1‑115 (MQ=255)                                    aGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATGAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGa                                                                                                  >  1:650639/1‑149 (MQ=255)                                            tCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCg                                                                                          >  1:493106/1‑149 (MQ=255)                                             cATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGt                                                                                         <  2:493106/149‑1 (MQ=255)                                              aTTACAGTTCTCCACTATTGCCAGTTCAGCGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGCCCTGACGCCCACCTGCGTATACCGCTAGTACACCGGCGATATCCGGCACGTCCCGGCCGGGGCAATTACAGACGCACCCGtc                                                                                        >  2:504598/1‑148 (MQ=255)                                               ttACAGTTCTCCACTATTGTCAGTTCAGGATCATCCGCATCAAGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGGTATCCGGAACGTCCCTGCCGGTGTGATGACAGACGCTACAGTGc                                                                                       >  1:199476/1‑149 (MQ=255)                                                    gTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTACAGGGCCCATTGCGGTCCTGGCGGCGTTCTGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTg                                                                                  <  1:721954/149‑1 (MQ=255)                                                             tATTGTCAGTTCAGGGTGATCCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGGTATCCGGAACGTCCCTGCCGGTGTAATGACAGCCGCTCCCGTGCCACTGGTATTTGcc                                                                         >  2:715074/1‑149 (MQ=21)                                                                 gTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTg                                                                     <  2:526386/149‑1 (MQ=255)                                                                  tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                         >  1:372084/1‑112 (MQ=255)                                                                  tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                         <  2:372084/112‑1 (MQ=255)                                                                    aGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                                                                                                         >  2:238208/1‑78 (MQ=255)                                                                    aGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTc                                                                                                                                         <  1:238208/78‑1 (MQ=255)                                                                     gTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                           >  2:702660/1‑107 (MQ=255)                                                                     gTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGgtgt                                                                                                           <  1:702660/107‑1 (MQ=255)                                                                                  cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                      <  2:205610/67‑1 (MQ=255)                                                                                  cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCgg                                                                                                                                      >  1:205610/1‑67 (MQ=255)                                                                                  cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                    <  1:506749/149‑1 (MQ=255)                                                                                     atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTCCACCGCTAGGACTCCGGCGATATCCGGAACGTCCCCGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                 <  2:278713/149‑1 (MQ=255)                                                                                     atcaGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGa                                                 <  1:438537/149‑1 (MQ=255)                                                                                          gcaggTATAGGGCCCATTGCGGTCTTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGATCGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCTGTAGGGAAAGATTCGGGACTCAg                                            <  1:761471/149‑1 (MQ=21)                                                                                              gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACGCAGGATg                                        >  1:144568/1‑149 (MQ=255)                                                                                              gTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGAGACTCAGGGTg                                        >  1:144300/1‑149 (MQ=255)                                                                                               tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                       <  2:650639/149‑1 (MQ=255)                                                                                                         aTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTGAGGATGCTGGCGGtata                             >  2:810260/1‑149 (MQ=255)                                                                                                           tGCGGCCCTGGCGGCGTTCGGCGTGTACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                           <  1:463841/149‑1 (MQ=255)                                                                                                           tGCGGCCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                           <  2:807346/149‑1 (MQ=255)                                                                                                                      cggcgTTCTGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                <  2:262566/149‑1 (MQ=255)                                                                                                                               gcgcacaCCGCTAGGTCTCCGGGGAGATCCGGACCGTCCCTGCCGGGGTAATGACAGACGCCACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACt       <  1:258432/143‑1 (MQ=14)                                                                                                                                    tcccGCTAGGACTCCGGAGATATCCGGAACGTCCCTGGCGGTGTAAGGACAGACGCTACCGTGCCGCTGGTATTTGCCGGTGCAGTAGCGTAAGATTCGGGGCTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa  <  1:502255/147‑1 (MQ=14)                                                                                                                                      ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                    <  2:501291/97‑1 (MQ=255)                                                                                                                                      ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                    >  1:501291/1‑97 (MQ=255)                                                                                                                                          |                                                                                                                                               AGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAA  >  NC_000913_3_hsa_tpiA/263961‑264240

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