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

GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATG  >  NC_000913_3_hsa_tpiA/4604560‑4604827                                                                                                                                           |                                                                                                                                  gacTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTCAAACGCAATCCACAAACGATCTGTGAACATAAAttt                                                                                                                         >  1:530663/1‑149 (MQ=255) gacTGAACAATCAGTGGTAAATTAACCAAAGCCGGACGCAAAGCCCGACGGAGAGAGAAATAATGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAACTTAAACGCAAAATCACAACGATCCGTGAAAATCAAttt                                                                                                                         >  2:147147/1‑149 (MQ=255)  acTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGGGCTTTATTTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAAGGCAATATCAAAACGATGCGTGAACATGAAtttg                                                                                                                        >  2:143138/1‑148 (MQ=255)    tGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAGATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAt                                                                                                                      >  1:290581/1‑149 (MQ=255)        cAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTCACTGCGATAGCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                  <  1:564685/149‑1 (MQ=255)         aaTCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGAt                                                                                                                 >  2:559310/1‑149 (MQ=255)           tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                               >  1:296232/1‑149 (MQ=255)                ggaaCATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGCGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGaaa                                                                                                          <  1:538482/146‑1 (MQ=255)                      ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                    >  1:321944/1‑149 (MQ=255)                      ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                    >  1:56681/1‑149 (MQ=255)                      ttAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGc                                                                                                    >  1:464372/1‑149 (MQ=255)                          ccAAAGCCGGCCGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATTAATTTATTGTGATGTCGAAACTGGGCGACa                                                                                                >  1:342738/1‑149 (MQ=255)                           cAAAGCCGGACGTTAAGGCCGACGGAGAGAGAAATATTGCGCTTTATGGCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATAGCAAAACGATCCGTGAAAATAAATTTAATGGGATGTCAGAACTGGGCGTTAg                                                                                               >  2:162828/1‑149 (MQ=255)                                 cGGATGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                         <  2:239363/149‑1 (MQ=255)                                 cGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGa                                                                                         <  1:34381/149‑1 (MQ=255)                                      gTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAg                                                                                    >  1:518426/1‑149 (MQ=255)                                        tAAGCCCGACGGAGGGAGAAATATTGCGATTTAGGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATCAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                  <  1:67124/149‑1 (MQ=255)                                        tAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAggg                                                                                  <  2:276590/149‑1 (MQ=255)                                           gCCCGACGGAGAGAGAAATATTGCGCTATATGTCGCGTGGCTACGCAATGACACAAATTGCCGAGCAACTTAGACGCAATATCAAAACGATCCGTGACCATAAAGTTAATGTAATGTCGAAACTCGGCGTCAGTTTTGAAGCGGGcttg                                                                               >  2:42957/1‑149 (MQ=255)                                              cGACGGAGAGAGAAAAGTGTCGCTTTATGTCGCGTGGCTAATCAATGATACAAATTGCCGAGCCGCCTAAACGCAATATCAAAACGATCCGTGAACATACATTCAATGTGATGTCGAAAGTGGGCGTCAGTTGTGACGCAGGGttgttg                                                                            <  2:516592/149‑1 (MQ=255)                                                    agagagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGATCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGcc                                                                      <  1:200173/149‑1 (MQ=255)                                                      agagaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGAGTCAGTTATCACGCAAGGCTGTTGGAAGCCga                                                                    >  2:278474/1‑148 (MQ=255)                                                        agaAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGAGGTCAGTTCTGACTCAGGGTTGTTGGAGGCCGctg                                                                  >  2:211913/1‑147 (MQ=255)                                                                                  gctatgaaaTGACACAAATTGCAGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGTGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGc                                        <  2:539731/142‑1 (MQ=255)                                                                                     aCTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAGGTAATGTCGAAACTGGGCGTCAGTTCTAACGCCGGATGATTGAAGGCAGCGGATATTCTTTTATGTCAGCGGCATTGccga                                     >  2:199216/1‑146 (MQ=255)                                                                                         aaTGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACaa                                 <  1:231282/149‑1 (MQ=255)                                                                                             acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                             >  1:395783/1‑149 (MQ=255)                                                                                               acaAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtg                           >  2:184239/1‑149 (MQ=255)                                                                                                  aaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATgtgt                        >  1:292088/1‑149 (MQ=255)                                                                                                   aTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGtt                       <  1:354779/149‑1 (MQ=255)                                                                                                    gtGCCAAGCAGCCGAAACACAAAACCAAAACGATACGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCACATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTg                      <  1:460147/148‑1 (MQ=255)                                                                                                            cagaTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGACAGCAAATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTa              <  2:377288/149‑1 (MQ=255)                                                                                                                  aaaCGAAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTTGCAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATc        <  1:415759/149‑1 (MQ=255)                                                                                                                      gCAATATCAAAACGAGCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCa    <  1:52264/149‑1 (MQ=255)                                                                                                                       cAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATGCCTATTAATACGCat   >  1:142553/1‑149 (MQ=255)                                                                                                                        aaTATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACCAGTAATGTGTGGCATACCTATTAATCCGCatg  >  2:411716/1‑149 (MQ=255)                                                                                                                                           |                                                                                                                                  GACTGAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATG  >  NC_000913_3_hsa_tpiA/4604560‑4604827

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