Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 2,173,361 Δ2 bp 100% pseudogene (1‑2/442 nt) gatC ← pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,173,3610C.94.6% 153.4 / 1.7 37pseudogene (2/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (18/17);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.
*NC_000913_3_hsa_tpiA2,173,3620C.94.6% 150.2 / ‑0.2 37pseudogene (1/442 nt)gatCpseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base . (18/17);  minor base A (0/1);  total (18/19)
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.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGG  >  NC_000913_3_hsa_tpiA/2173215‑2173498
                                                                                                                                                  ||                                                                                                                                        
gAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATCAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccg                                                                                                                                         >  2:592446/1‑148 (MQ=255)
      aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgtat                                                                                                                                   >  2:584234/1‑142 (MQ=255)
                tGCACGCCGACCGCCATCGCCCCGAAGAAGCCGATGGTGGCAAGTTCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCcagcag                                                                                                                       <  1:52317/149‑1 (MQ=255)
                 gCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTCGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCagcagc                                                                                                                      >  2:774827/1‑149 (MQ=255)
                   aCGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCG‑‑ACCCCACCGCCGTATGTCCCAGCAGCaa                                                                                                                    <  2:424807/149‑1 (MQ=255)
                        gACCTCCAACGCCAGAAAGAACCCCAGGGCTGGAAGAGCGTCAAACGACAGCACCTGATTACCCGGCACACAAAGAGCAAGTAAAATGGTGAGTGGGATAAAAATGAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                               <  2:342326/149‑1 (MQ=255)
                        gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAACATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTCGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                               <  1:518836/149‑1 (MQ=255)
                          ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                             <  1:207301/149‑1 (MQ=255)
                           cggcagcGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCGCAGCAGCAACGCCGGAt                                                                                                            <  2:199656/143‑1 (MQ=255)
                                tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAgg                                                                                                       >  1:418655/1‑149 (MQ=255)
                                    cACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCaa                                                                                                   >  1:306250/1‑149 (MQ=255)
                                               cGATGGTGGCAAGATCGCTAAACGGCAGCAACTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTcc                                                                                        <  2:103315/149‑1 (MQ=255)
                                                       gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgcc                                                                                >  1:82543/1‑149 (MQ=255)
                                                           gATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGTCCAATCAGGAACTCCTGACCGCCGAAc                                                                            >  1:201477/1‑149 (MQ=255)
                                                              ccccAAACGGCAGCACCTTATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                         <  2:318223/147‑1 (MQ=255)
                                                               gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCg                                                                        >  1:776556/1‑149 (MQ=255)
                                                                ccAGACGGCAGGCCCTGAGTACCCGGCGCACAAACAGCGATGAAAATGGTGAGTGGGATAAAAATCGCGCTTGGCGATAGCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAGCAGGAACTCCTGACCGCCGAACTTCGc                                                                       <  2:482446/149‑1 (MQ=255)
                                                                 taaaCGGCCGCACCTTATTACTCGGCCCACAAACACTAATTAAAAGGGTGAGTGGGATAAAAATCAGGCTTGCCGACACCA‑‑CCGCCGTATGTCCCAGGAGCAACGCCGGATCAATGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                      <  2:83588/148‑1 (MQ=255)
                                                                       gcagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGAGGCTGAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                                <  1:24434/149‑1 (MQ=255)
                                                                        cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAAc                                                               <  1:58900/149‑1 (MQ=255)
                                                                          gcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAACTAGCAAACCCCGACCGCGGAACTTCGCCTGTAAACGa                                                             >  1:530027/1‑149 (MQ=255)
                                                                                        ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTACCCGACTCCGATCcttctt                                               >  1:93894/1‑149 (MQ=255)
                                                                                        ggCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACCCCGAGCCTGCtt                                               >  2:251314/1‑149 (MQ=255)
                                                                                          cacaaaAACTGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTTACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg                                             <  2:515777/149‑1 (MQ=255)
                                                                                           acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                            >  1:30803/1‑149 (MQ=255)
                                                                                             acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGa                                          >  1:575910/1‑149 (MQ=255)
                                                                                                     cAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGtt                                  >  2:759887/1‑149 (MQ=255)
                                                                                                      aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                 <  1:669180/149‑1 (MQ=255)
                                                                                                           aaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACAA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAACACGAGCCTGCTTGGCGATGGGCGTTAAAcca                            >  1:296444/1‑149 (MQ=255)
                                                                                                             aTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc                          <  1:674843/149‑1 (MQ=255)
                                                                                                                gTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatccat                       <  2:418655/149‑1 (MQ=255)
                                                                                                                  gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                     <  1:543116/149‑1 (MQ=255)
                                                                                                                     tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGGCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                  <  1:581597/149‑1 (MQ=255)
                                                                                                                     tGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTACTGCCCGCCGTAGTTCGCCTGTAAACGACAGTGGGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg                  <  2:809785/149‑1 (MQ=255)
                                                                                                                        gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAAGCCAATCGGACACTCCTGACCGCCGACCTTCGCCCGCAACTGTCAACTACCCGGCGTACCGATGGGCGTTAAACAATCCGTGATGGGtt               >  2:174717/1‑149 (MQ=255)
                                                                                                                         aTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGttt              >  2:474425/1‑149 (MQ=255)
                                                                                                                                    gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGAAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCAGGCTGGGTTTAATCACCCCTg   >  2:184046/1‑149 (MQ=255)
                                                                                                                                    gCTTGCCGATACCA‑‑CCGCCATATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg   >  2:224365/1‑149 (MQ=255)
                                                                                                                                     cTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAAGCCAGAGCAAGGCCAACCAGTTAGCCCGGACCGCCGAACTTCGACTGTAAACGACTACGAGCCTGCTTAGCGATGGCCGTTAAACCATCCATGATGGGTTTAAACACCCGTgc  >  2:726729/1‑148 (MQ=255)
                                                                                                                                                  ||                                                                                                                                        
GAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGG  >  NC_000913_3_hsa_tpiA/2173215‑2173498

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