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%	143.8 / NA	41	A186E (GCA→GAA)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (16/25);  total (16/25)

ACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_hsa_tpiA/4604561‑4604834                                                                                                                                          |                                                                                                                                         acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTaa                                                                                                                              >  1:36012/1‑150 (MQ=255)       cAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATTACAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                        <  2:144887/150‑1 (MQ=255)          tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTc                                                                                                                     >  1:344189/1‑150 (MQ=255)              tGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                                 <  1:377908/150‑1 (MQ=255)                gTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACt                                                                                                               <  2:364579/150‑1 (MQ=255)                  ataTTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTgg                                                                                                             <  1:285277/148‑1 (MQ=255)                     ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTTAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCACAACGATCCGTGAACATAGATTTAATGTGATGTACAAACTGGGCg                                                                                                          >  1:541785/1‑150 (MQ=255)                           aaggCCTGATGTTAAGCCCGACGGAGCGAGAAATATTGCGGTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGtt                                                                                                    <  2:582108/147‑1 (MQ=255)                             aGCCGGACGTTAAGCCCGACGGAGATAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                                                                                  <  2:528544/150‑1 (MQ=255)                             aGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCt                                                                                                  <  1:412669/150‑1 (MQ=255)                                   aCGTTAAGCCCGACGTAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                                            <  1:85606/150‑1 (MQ=255)                                   aCGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                                            <  2:209349/150‑1 (MQ=255)                                        aaGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCTAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgtt                                                                                       <  1:575734/150‑1 (MQ=255)                                           cccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGgttgtt                                                                                    >  2:535084/1‑150 (MQ=255)                                            ccGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                   <  2:36012/150‑1 (MQ=255)                                            ccGACAGAGAGAAAAAAATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCAAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGttgttg                                                                                   <  2:451176/150‑1 (MQ=255)                                              gACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                 <  1:535084/150‑1 (MQ=255)                                                cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                               <  2:356218/150‑1 (MQ=255)                                                cGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATCGCCGAGCAGCTTAAACGCAATCTCAAAACGACCCGTGAACATAAATTTAATGTGCTGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAgg                                                                               <  2:627554/150‑1 (MQ=255)                                                  gagagagaAATATTGCGCTTTAGGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGcc                                                                             <  1:239460/150‑1 (MQ=255)                                                     agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCa                                                                          <  1:294233/150‑1 (MQ=255)                                                       agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTAAACATAAAATTCATGTGAAGCCGAAACTGGGCGTCAGCCCTAACGCAGGGTTGTTGGAGGCCGCAGa                                                                        >  2:440063/1‑150 (MQ=255)                                                                    tttATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTtatgt                                                           <  1:69826/150‑1 (MQ=255)                                                                          tCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATGCTGTTATGTATGCgg                                                     >  1:604905/1‑150 (MQ=255)                                                                            gcgTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAAGGATGTTGTGCGCCGCAGATATTCTGTTATGTATGCggga                                                   >  1:617619/1‑148 (MQ=255)                                                                              gTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCACTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGATGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCAtt                                                 >  1:410719/1‑150 (MQ=255)                                                                                ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTCATGTGATGACGAGAGTGGGCGGCAGTTCTGACGCAGGGTTGTTGGAGGCCGCGGATATTCTGTTATGTATGCGGCATTGc                                               >  2:72863/1‑150 (MQ=255)                                                                                   tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                          <  1:46494/104‑1 (MQ=255)                                                                                   tACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAgg                                                                                          >  2:46494/1‑104 (MQ=255)                                                                                      tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                       <  2:199280/88‑1 (MQ=255)                                                                                      tCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCa                                                                                                       >  1:199280/1‑88 (MQ=255)                                                                                       aaaTCACTCAATTTGCCGAGTAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGGGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTGGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                        <  1:615097/149‑1 (MQ=255)                                                                                             cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                 >  2:460583/1‑103 (MQ=255)                                                                                             cacaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGa                                                                                 <  1:460583/103‑1 (MQ=255)                                                                                                 aaTTGACGCGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                              <  2:482383/150‑1 (MQ=255)                                                                                                           cagcTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTAt                    >  2:538624/1‑150 (MQ=255)                                                                                                              cTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTaa                 >  1:90481/1‑150 (MQ=255)                                                                                                                          atCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGATGGAGGCCGCAGATATTCTGTTATGTATGCGGCACTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGc     >  1:503381/1‑150 (MQ=255)                                                                                                                             aaaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg  <  1:538624/150‑1 (MQ=255)                                                                                                                              aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                    >  2:31202/1‑99 (MQ=255)                                                                                                                              aaaCGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGc                                                    <  1:31202/99‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         ACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_hsa_tpiA/4604561‑4604834

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