Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,173,363 Δ2 bp intergenic (‑1/+1) gatC ← / ← gatC pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC/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.95.2% 85.3 / ‑3.0 21pseudogene (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 (1/0);  new base . (10/10);  total (11/10)
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
*NC_000913_3_hsa_tpiA2,173,3620C.100.0% 88.2 / NA 20pseudogene (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/0);  new base . (10/10);  total (10/10)

AACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913_3_hsa_tpiA/2173216‑2173497
                                                                                                                                                 ||                                                                                                                                       
aaCAGATTTCCACGATGCACGGCGACCGACATCGCCACGAAGAAGCCGATGGTGGCACGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgc                                                                                                                                       >  2:573445/1‑147 (MQ=255)
 aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGATGTATGGCTAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGATGGGGATAAAAATCAGGCTTGCCGATACCACcgcg                                                                                                                                      >  1:270364/1‑149 (MQ=255)
   aGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgt                                                                                                                                    <  1:573915/149‑6 (MQ=255)
     aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACCTCAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccgccgcat                                                                                                                                  >  2:574285/1‑142 (MQ=255)
                             cATCGCCACGAAGAAGGCGATGGTGGCAAGATCGCCAAACGCCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa                                                                                                        <  2:381706/149‑1 (MQ=255)
                              aTCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                       <  1:309414/149‑1 (MQ=255)
                                        agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                             >  1:334008/1‑149 (MQ=255)
                                          aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGaa                                                                                           <  2:347498/149‑1 (MQ=255)
                                                     ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccgc                                                                                <  2:327942/149‑1 (MQ=255)
                                                             ctccAAAACGCAGCACCTGATTACCCCGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                        <  2:270364/147‑1 (MQ=255)
                                                                cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                     >  1:200121/1‑149 (MQ=255)
                                                                cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAACAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc                                                                     <  2:155648/149‑1 (MQ=255)
                                                                 aaaGGGCAGCACAGAATTCCCCGGTAGACAAACAGCAATAAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt                                                                    <  2:518261/149‑1 (MQ=255)
                                                                   aCGGTAGCACCTGATTACCCGGCCCACAAACAGCAATTAAAATGGTGAGTGCGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCTCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGt                                                                  <  2:572021/149‑1 (MQ=255)
                                                                              tGATTACCCGGCACACAAACAGCAATTAAAATCGTCGCACGGAGAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCCCCGAACTTCGCCTGTAAACGACTACg                                                       >  1:145778/1‑149 (MQ=255)
                                                                                      cGGCACAAAAACAGCAATTAAAACGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGGCCCAGCAGCAACGCCAGATGAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCACCTGTAAACGACTACAAGCCTGCt                                               >  2:557393/1‑149 (MQ=255)
                                                                                          acacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGc                                           <  1:573445/149‑1 (MQ=255)
                                                                                            acaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTACGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAGACGACTACGAGCCTGCTTAGCGa                                         >  1:512136/1‑149 (MQ=255)
                                                                                                aCAGCAATTAAAATGGGGAGTGGGATAAAAATCTGGCTTGCCGCTACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCTCAAGGCCAATCAGGAACTCCTGACCGCCGAACCTCCGCTGTCAGCGACGTCGGGCCTGATTAGAGATgac                                     >  1:283523/1‑147 (MQ=255)
                                                                                                 cAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGc                                    >  1:466109/1‑149 (MQ=255)
                                                                                                     aaTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTa                                >  1:52479/1‑149 (MQ=255)
                                                                                                                            aaaTCAGGCCTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGTGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc         <  1:314902/149‑1 (MQ=255)
                                                                                                                              aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAAATTCGCCTGTAAACGCCTACTATCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc       >  1:267602/1‑149 (MQ=255)
                                                                                                                                   gCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGCGTTTACTCACCCGTg  >  2:357990/1‑149 (MQ=255)
                                                                                                                                                 ||                                                                                                                                       
AACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTG  >  NC_000913_3_hsa_tpiA/2173216‑2173497

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: