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_tpiA4,604,6980CA94.7% 58.4 / NA 19A186E (GCA→GAA) bglJbgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (10/8);  minor base G (1/0);  total (11/8)
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

GAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_hsa_tpiA/4604564‑4604834
                                                                                                                                      |                                                                                                                                        
gAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtg                                                                                                                            >  1:184663/1‑149 (MQ=255)
 aaCAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAAtgt                                                                                                                           <  2:384186/149‑1 (MQ=255)
   caagtaGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGa                                                                                                                         <  2:496725/144‑1 (MQ=255)
      tCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATCGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGt                                                                                                                      >  2:434206/1‑149 (MQ=255)
              aCATCAACCAAAGCCTGACGTTAACCCCTACGGAGAGAGAAAAATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATGAATTTAATGTGATGTCGAAACTg                                                                                                              <  1:110895/149‑1 (MQ=255)
                   aaCCAAAGCCGGACGTTAAGCCCGACGGAGGGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCCATATCAAAACGATCAGTGGACATAAATTTAATGTGATGTCGAAACTGGGCGt                                                                                                         >  2:35865/1‑149 (MQ=255)
                               aCGTGAAGCCCGACGGAGAGAGAAATATTGCGCCTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGc                                                                                             <  1:41129/149‑1 (MQ=255)
                                cGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGCTCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGGCTGTTCTGACGCa                                                                                            >  2:55854/1‑149 (MQ=255)
                                                            gcgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGt                                                                <  1:269694/149‑1 (MQ=255)
                                                             cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                                            <  1:195054/120‑1 (MQ=255)
                                                             cgcTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCa                                                                                            >  2:195054/1‑120 (MQ=255)
                                                                            ggCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCCGATATTCTGTTATGTCTGCGGcgtcc                                                >  2:76255/1‑145 (MQ=255)
                                                                                        acacAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTaa                                    >  1:468163/1‑149 (MQ=255)
                                                                                            aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCATCAGTTCTGACGCAGGGTTGTTGGAGACCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAAtgtg                                >  2:168111/1‑149 (MQ=255)
                                                                                            aaaTTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGAACATAAATGTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGGTGTTGGAGGCCGCAGATCTTCTGTTATGTATGCGGCATTTCCAAACAAGTACtgtg                                >  2:164625/1‑149 (MQ=255)
                                                                                                               aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg             <  2:262495/149‑1 (MQ=255)
                                                                                                               aCGCAATATCAAAACGATCCGTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCg             >  1:403549/1‑149 (MQ=255)
                                                                                                                                   gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg  <  2:434001/140‑1 (MQ=255)
                                                                                                                                   gTGAACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCAGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCgg  >  1:434001/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GAACAATCAGTGGTACATTAACCAAAGCCGGACGTTAAGCCCGACGGAGAGAGAAATATTGCGCTTTATGTCGCGTGGCTACTCAATGACACAAATTGCCGAGCAGCTTAAACGCAATATCAAAACGATCCGTGCACATAAATTTAATGTGATGTCGAAACTGGGCGTCAGTTCTGACGCAGGGTTGTTGGAGGCCGCAGATATTCTGTTATGTATGCGGCATTGCGAAACAAGTAATGTGTTGCATCCCTATTAATCCGCATGATGCCGG  >  NC_000913_3_hsa_tpiA/4604564‑4604834

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