Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_tpiA	4,604,698	C→A	A186E (GCA→GAA)	bglJ →	bgl operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_tpiA	4,604,698	0	C	A	100.0%	81.2 / NA	24	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/11);  total (13/11)

TGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTT  >  NC_000913_3_hsa_tpiA/4604563‑4604838                                                                                                                                         |                                                                                                                                             tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCAGTGAACATAAATTTAAt                                                                                                                                 >  1:28144/1‑149 (MQ=255)   aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                               >  1:155613/1‑149 (MQ=255)                   ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAAAAATAGTTCTGTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                               <  1:188889/149‑1 (MQ=255)                          aaGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                        >  2:16862/1‑149 (MQ=255)                             ccGGAGGTTAAGCCCGACTGAGAGCGAAATATTGCGCTTTATGTCGCGTGGCTACGCAATGATACAAATTGCGGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATACATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTg                                                                                                     <  2:78101/149‑1 (MQ=255)                                        gCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttg                                                                                          >  1:240061/1‑149 (MQ=255)                                         cccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgt                                                                                         <  1:16862/149‑1 (MQ=255)                                          cagACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGTTCCGTGAACATCAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                        <  1:142507/147‑1 (MQ=255)                                           cGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTAACGCAGGGttgttg                                                                                       >  1:184897/1‑149 (MQ=255)                                             aCGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAAGGTGATGTCGAAACTGAGCGTCAGTTATGACGCAAGGTTGTTGGa                                                                                     >  1:325650/1‑149 (MQ=255)                                               ggAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                                   <  1:365888/149‑1 (MQ=255)                                                     agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAg                                                                             >  1:149647/1‑149 (MQ=255)                                                                 cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                 <  2:149647/149‑1 (MQ=255)                                                                 cTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtat                                                                 <  2:419980/149‑1 (MQ=255)                                                                                    tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAAc                                              >  1:289693/1‑149 (MQ=255)                                                                                          acacAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                        >  2:74923/1‑149 (MQ=255)                                                                                            acaAATTGCCGAGCAGCTTAAACGCAATATACCACCCGCTTTG‑GAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGTCGGCCGCCGATATGCCGGTATGTATGCGGCATGGAGAAACAAGTAAtg                                      >  1:103486/1‑149 (MQ=255)                                                                                                      gAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAAGGTTGTAGGGGTCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATcc                            >  1:315970/1‑149 (MQ=255)                                                                                                        gcagcTTAAACGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCt                          <  1:88022/149‑1 (MQ=255)                                                                                                         cagcTTAAACGAAATAT‑CAAACCGGTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa                         <  2:184897/149‑1 (MQ=255)                                                                                                                 aCGCAATAT‑CAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg                 >  1:154754/1‑149 (MQ=255)                                                                                                                              aaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCggg     <  2:240061/149‑1 (MQ=255)                                                                                                                                 gATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGttt  <  1:74923/149‑1 (MQ=255)                                                                                                                                 gATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGttt  >  2:287613/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                             TGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAT‑CAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGGGTTT  >  NC_000913_3_hsa_tpiA/4604563‑4604838

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