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.95.8% 97.9 / ‑1.7 24intergenic (‑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 . (12/11);  minor base A (0/1);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCA  >  NC_000913_3_hsa_tpiA/3815661‑3815951
                                                                                                                                                    |                                                                                                                                              
cttGAAGAACTAGGGGCTTGTGCGGCCAGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAg                                                                                                                                                <  1:224759/148‑1 (MQ=255)
   gacgaaATCGGGGCTTCTGCGCCCGGCTTTCAGCGTCACCCCGCCACCCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAgcga                                                                                                                                             <  1:602337/146‑4 (MQ=255)
               gCTTTTGCGCCCGGATTTCAGTGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGc                                                                                                                                <  1:182544/149‑1 (MQ=255)
                       gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACACAAAAAAGCGACCCATCCGTCGCCTGaaaaa                                                                                                                        >  1:248032/1‑149 (MQ=255)
                       gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTaaaaa                                                                                                                        >  2:585867/1‑149 (MQ=255)
                          cGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                     <  1:292702/149‑1 (MQ=255)
                                         aCTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgccgc                                                                                                      >  2:335214/1‑149 (MQ=255)
                                               cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAAATGGCGCTGAAAAGGTTTCATGCAATCGATCCTCATCTTATGTTTCTGCAGACAAAAAAAAGCGTCCCATCATCCGCCATAAATTTCAGTTTTCCAGCGCCGCCTTCTg                                                                                                >  2:63892/1‑149 (MQ=255)
                                                         cacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCAGATATGGTTTCATGCCTTCGCTCCTCATCTCACTTTTCTACAGACAAAAACAAGCGACTCATCAGTCGCCTTAACAATCAGTTTGCAAGCCACCCCTTCTGCGGCCCTACa                                                                                      >  2:286031/1‑149 (MQ=255)
                                                          acCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTCCaa                                                                                     >  1:571772/1‑149 (MQ=255)
                                                           ccTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAt                                                                                    <  2:216145/149‑1 (MQ=255)
                                                              gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGga                                                                                 >  2:468925/1‑149 (MQ=255)
                                                               ctggctACGCGCAAATTCAATAAACGGGTGCTTATATGGTTTCATGCCTTCGTTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgat                                                                                <  2:136612/146‑1 (MQ=255)
                                                                ttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                               <  2:307679/149‑1 (MQ=255)
                                                                      aGCGCAAATTCAATAAACTGGCGCTTATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgattc                                                                         <  2:200345/149‑1 (MQ=255)
                                                                               tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGc                                                                >  2:442797/1‑149 (MQ=255)
                                                                                      aCTGGCGCTGATATGGTGTCATGCCTTCGCTCCTAATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTGCAGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGcc                                                         >  2:569407/1‑149 (MQ=255)
                                                                                       cTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                                        <  1:63892/149‑1 (MQ=255)
                                                                                                atttGGTTTCCTGCCGTCGCTCCTCATCTTACTTTTCTACAGACAAAACAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCTTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGa                                               <  1:79894/146‑1 (MQ=255)
                                                                                                    ggTTTCATGCCTTCGCGCCTCATCTTACTGTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTCGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                           <  2:487470/149‑1 (MQ=255)
                                                                                                    ggTTTCATGCCTTCCCTCCTCATCTTACTTTTCGACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                           <  1:297788/149‑1 (MQ=255)
                                                                                                      tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGct                                         >  2:44331/1‑149 (MQ=255)
                                                                                                        tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACCAACAAAAGCGACTCATCAGCCGCCATAACACCCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGctct                                       >  1:504664/1‑149 (MQ=255)
                                                                                                          aTGCCTTCGCTCCTCATCTTACTTTATTACAGACAAAAAAAAGCGACTCATCGGTCGCCTTAAAAATCAGTTTGACAGCGCCGCCTTCTGCGGCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTc                                     <  2:582493/149‑1 (MQ=255)
                                                                                                                    tcctcATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCCCCGCCTTCGGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGAAGCTCTTCATGGGTGaaa                           >  1:461108/1‑148 (MQ=255)
                                                                                                                                            caaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGTCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg    <  1:585867/148‑1 (MQ=255)
                                                                                                                                            caaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg    >  1:116400/2‑149 (MQ=255)
                                                                                                                                              aaaaaaAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCa  <  1:395717/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTTGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCA  >  NC_000913_3_hsa_tpiA/3815661‑3815951

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: