Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 59,074 C→A pseudogene (23/228 nt) yabP → pseudogene, pentapeptide repeats‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA59,0740CA94.4% 53.3 / ‑4.6 18pseudogene (23/228 nt)yabPpseudogene, pentapeptide repeats‑containing
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (7/10);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTG  >  NC_000913_3_hsa_tpiA/58928‑59217
                                                                                                                                                  |                                                                                                                                               
gCGTTGATCTTAGAAAGCCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATCCTAACCTGGGCgtgctaccccggacggtgctaacc                                                                                                                                               >  2:957090/1‑149 (MQ=255)
  gTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaacccg                                                                                                                                             >  2:1126299/1‑149 (MQ=255)
    ttatattAGAAACGGTTACCGTGCCGGCGTTGAAATGAGGGATTGTTACCCTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaacccgga                                                                                                                                           <  1:95/144‑1 (MQ=255)
    tGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacgttgctaacccgga                                                                                                                                           <  2:1016619/149‑1 (MQ=255)
      agcTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaacccggacg                                                                                                                                         >  1:536396/3‑149 (MQ=255)
      aTCTTAGAAAGGCTTACCTTCCCGGTGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCgtgctaccccggacggtgctaacccggacg                                                                                                                                         <  1:1152405/149‑1 (MQ=255)
         ttAGAAAGGCTTACCTGCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGctaccccggacggtgctaacccggacggtg                                                                                                                                      <  2:1058248/149‑1 (MQ=255)
          tAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCtaccccggacggtgctaacccggacggtgc                                                                                                                                     >  1:1125186/1‑149 (MQ=255)
             aaaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGTTAGCAATGATATCCTCTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGCCGTAGCTAAcccggacggtgctaa                                                                                                                                  >  2:690664/1‑149 (MQ=255)
             aaaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAAGGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccgggcggtgctaa                                                                                                                                  >  1:788924/1‑149 (MQ=255)
              aaGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCccggacggtgctaacccggacggtgctaac                                                                                                                                 <  1:1046662/149‑1 (MQ=255)
                     aCCCTCCCCGCGTTGACATTAGGGATTGTTACCGTGGTAAAAAAACAATGGAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacg                                                                                                                          <  1:542996/149‑1 (MQ=255)
                      ccTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCccggacggtgctaacccggacggtgctaac                                                                                                                                 <  1:1095424/141‑1 (MQ=255)
                      ccTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCccggacggtgctaacccggacggtgctaac                                                                                                                                 >  2:1095424/1‑141 (MQ=255)
                        ttCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacggtg                                                                                                                       >  2:464726/1‑149 (MQ=255)
                           ccGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacggtgcta                                                                                                                    >  1:1148968/1‑149 (MQ=255)
                               cGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacggtgctaccgt                                                                                                                <  1:123247/149‑5 (MQ=255)
                                 ttGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGGAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacggtgctaccgtga                                                                                                              >  1:1054918/1‑143 (MQ=255)
                                  tGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAAcccggacggtgctaacctggacggtgctaccgtgaa                                                                                                             <  1:952249/149‑8 (MQ=255)
                                        gAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGc                                                                                                       <  2:600091/149‑1 (MQ=255)
                                          gggATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTa                                                                                                     <  1:1022203/149‑1 (MQ=255)
                                          gggATTGTTAACTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTa                                                                                                     <  2:810747/149‑1 (MQ=255)
                                               tGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcc                                                                                                <  1:732002/149‑4 (MQ=255)
                                                 ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt                                                                                              <  1:69001/149‑6 (MQ=255)
                                                 ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt                                                                                              <  1:655847/149‑6 (MQ=255)
                                                 ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt                                                                                              >  2:386339/1‑144 (MQ=255)
                                                 ttACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctt                                                                                              >  2:645900/1‑144 (MQ=255)
                                                  tACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTACTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTAcctcctta                                                                                             <  1:1077697/149‑7 (MQ=255)
                                                                           gTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATtaata                                                     >  2:1139388/1‑149 (MQ=255)
                                                                                                                                             ggctacaccggacggggcgaGACCGGACGGTGATAACCTGGAAGGAGCTACCGTGAACGGCGCATCACCATTAGTCGTTGAGGTAATCTTTATTAATAAAAAAAAACCAAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTg  <  2:720018/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGTTGATCTTAGAAAGGCTTACCTTCCCGGCGTTGACATGAGGGATTGTTACCTTGGTAAAAAAACAATGAAAGGTAGCAATGATATCCTTTATGAGAGACCTGGGTGGAATGCTAACCTGGGCGTGCTACCCCGGACGGTGCTACCCCGGACGGTGCTAACCCGGACGGTGCTAACCTGGACGGTGCTACCGTGAACGGTGCTACCTCCTTATATGATGAGGTAATTATTATTAATAAAATCCCCCCCAAAAAAATTGATACTAAAGGAGTTGCTACTGAAGAAGTTG  >  NC_000913_3_hsa_tpiA/58928‑59217

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: