Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 3,815,810 Δ1 bp intergenic (‑42/+24) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA3,815,8090G.92.9% 104.3 / 3.0 28intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (9/17);  minor base A (1/1);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

TTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC  >  NC_000913_3_hsa_tpiA/3815662‑3815947
                                                                                                                                                   |                                                                                                                                           
ttGAAGAAATAGGGGCTTTTACGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCATTCCCTCCTCCTCTTACTTTTCTACAGACACAAAAAAga                                                                                                                                            >  1:301714/1‑148 (MQ=255)
gtGAAGAACTAGGGGCTTTTGCGCTCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAgc                                                                                                                                            <  1:841581/148‑2 (MQ=255)
   aagaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAgcgac                                                                                                                                         >  2:622298/1‑145 (MQ=255)
      aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACtca                                                                                                                                     >  1:608604/1‑149 (MQ=255)
      aaaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGGTTCATGCCTCCGCTCCTCATATTACTCTTCTACAGACAGAAAAAAGCGACtca                                                                                                                                     >  1:609442/1‑149 (MQ=255)
                        ccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCGATCGCCTTAAAATTc                                                                                                                   >  2:705924/1‑149 (MQ=255)
                                            gCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATCAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTc                                                                                               <  2:785160/149‑1 (MQ=255)
                                                       acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                    <  1:622298/149‑1 (MQ=255)
                                                         acCTGCCTGTTAAGCGCAAATTCTATAAACTGGTGCTGATATGGTTCCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCTGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACaa                                                                                  <  2:301714/149‑1 (MQ=255)
                                                             gcttggtAAGGGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                              <  2:101146/149‑1 (MQ=255)
                                                             gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGt                                                                                          <  1:500289/137‑1 (MQ=255)
                                                             gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGt                                                                                          >  2:500289/1‑137 (MQ=255)
                                                                  caaaGCGCAATTTAAATTATCTTGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACCCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                         <  2:134487/147‑1 (MQ=255)
                                                                            agtCAATAAACGTGCGCTAATATGAAATCAGGCCTTCGCTCCTCATCTGACTGTAAAACAAAAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                               <  2:82375/147‑1 (MQ=255)
                                                                            aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                               <  1:335799/149‑1 (MQ=255)
                                                                             ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                           <  1:250594/104‑1 (MQ=255)
                                                                             ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                           >  2:250594/1‑104 (MQ=255)
                                                                                             tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGAt                                              <  1:79977/149‑1 (MQ=255)
                                                                                               atatGGTTTCATGCCTTCGCTCCCCATCTTACTTTTCTACAGACAAAAAAAAGCGACTTAGCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTAGGGCCTTTGCCAACAAGATGa                                            <  2:734826/149‑1 (MQ=255)
                                                                                                    gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCACCAAGCTGGGCAGc                                      >  2:221390/1‑149 (MQ=255)
                                                                                                            ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATg                              >  2:137426/1‑149 (MQ=255)
                                                                                                                     ctcATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                               >  2:515955/1‑139 (MQ=255)
                                                                                                                     ctcATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCAt                               <  1:515955/139‑1 (MQ=255)
                                                                                                                        agtttACTTTTCTACAAACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                  <  2:210794/146‑1 (MQ=255)
                                                                                                                                 ttCTACAGACAAAAAAAAGCGACTCATCAGTCGTCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGACCACGCTGGTGTCATGGGTGAACGGCTCGCCTTCTGc         <  2:608604/149‑1 (MQ=255)
                                                                                                                                 aacgaacgACAAAGAAAAGCGACTATTCAGTGGTCTTAAAAATCAGTTTGCCAGCGCCGTCTTGTGTGTCGCTACAAGGGATTCGATTCCCCTCGGGCCAGAGCCAAGAAGATGTGCACTTGGTCGTGGGTGAACGGCTCGCCTTCTGc         <  2:609442/142‑1 (MQ=255)
                                                                                                                                  tCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc        <  1:137426/149‑1 (MQ=255)
                                                                                                                                  tCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATAGGTGAACGGCTCGCCTTCTGcc        <  1:783165/149‑1 (MQ=255)
                                                                                                                                     aCAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTc     >  2:82860/1‑149 (MQ=255)
                                                                                                                                       cgACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTccc   <  1:440657/148‑1 (MQ=255)
                                                                                                                                         aCAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCATTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc  >  1:560134/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCC  >  NC_000913_3_hsa_tpiA/3815662‑3815947

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: