Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_pgi 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_pgi2,173,3610C.96.2% 102.6 / ‑3.0 26pseudogene (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/1);  new base . (13/12);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
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_pgi2,173,3620C.96.2% 102.6 / ‑3.0 26pseudogene (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);  new base . (13/12);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
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.

ACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913_3_hsa_pgi/2173234‑2173493
                                                                                                                               ||                                                                                                                                   
aCGGCGACCGCCATCGCCACGAAGAAGCCGCTGGTGACAAGCTCGCCAAACGGCAGCACCGGATTACCCAGCACACAACCAGCAATTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCA‑‑CCGCCGGACGTCCcagcgcacc                                                                                                               >  2:300301/1‑144 (MQ=255)
     gACGGCCATCGGCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCg                                                                                                          <  2:20326/149‑1 (MQ=255)
       ccGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa                                                                                                        <  1:60385/149‑1 (MQ=255)
            agcGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg                                                                                                   <  2:616861/147‑1 (MQ=255)
                      agaagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGTTAACAGTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAgg                                                                                         <  1:427208/149‑1 (MQ=255)
                            cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTCAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCGGGAACTcc                                                                                   >  1:583993/1‑149 (MQ=255)
                            cGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCAAAAGGACAAGCAGGAACTcc                                                                                   >  1:252841/1‑149 (MQ=255)
                                  tggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGAccg                                                                             >  2:159071/1‑149 (MQ=255)
                                        gATCGCCAAACGGCCGCACCTGATTGCCCGGTACACAATCGGCAATTAAAATGGTTAGTGGGATAAAAGTCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCGATCAGGAACTCCTGACCGCCGAAc                                                                       <  2:179867/149‑1 (MQ=255)
                                           cGCCAAACTGCAGCACGTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc                                                                    <  2:413488/149‑1 (MQ=255)
                                             gtcaaCGCCAGCACGTGAGTGCGCGGCAGACAAACAGCAATTAAAATGTTTATTGGTATAAATCTGAGGCTTGCCGATACCA‑‑CCGCCGTATCGCCCACCAGCAACGCCGCATCAAGGCCAATCAGGAGCTCCTGACCGCCGAACTTCGc                                                                  <  2:2256/146‑1 (MQ=255)
                                             ccAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGc                                                                  >  2:130846/1‑149 (MQ=255)
                                                    gcagcaCCTGATTACCCGGCACACAAACAAAACAATGAATTCCTAGTGGGCTAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTaaa                                                           >  1:493260/1‑149 (MQ=255)
                                                       ggacgtGAGTACCGGGCACACAAACAGCAATGAAAATGGTGAGGGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGGATGGCCCGGCAGCAACGCCGGATCAAGGCCAATCAGGTGCTCCTGACCGCCGAACTTCGCCTGTAAACGa                                                        <  2:335828/144‑1 (MQ=255)
                                                          ccTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCAACATTAAACGACTa                                                     >  2:435860/1‑149 (MQ=255)
                                                           cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGGCTAc                                                    >  1:192173/1‑149 (MQ=255)
                                                                      gCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTACTTa                                         >  1:648929/1‑149 (MQ=255)
                                                                         cacaAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGCACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCg                                      >  1:412815/1‑149 (MQ=255)
                                                                                        tatgggTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAcca                       <  1:482868/145‑1 (MQ=255)
                                                                                               tggTGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTTCACGACGCGCCTCCGGATCAAGGCCAATCAGGAACTCCTTACCTCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                <  2:473284/147‑1 (MQ=255)
                                                                                               gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCCGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga                <  1:55177/149‑1 (MQ=255)
                                                                                                  tGGGATAAAAGTCAGGCTTGCCGATCCCA‑‑CCGCCGGATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCACGATgg             <  2:340926/149‑1 (MQ=255)
                                                                                                  tGGGATAAAAAGCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATgg             <  2:261386/149‑1 (MQ=255)
                                                                                                     gATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAGCGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt          >  2:347933/1‑149 (MQ=255)
                                                                                                       tAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtttg        >  1:140613/1‑148 (MQ=255)
                                                                                                             tCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCAAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTCGCGATGGGCGTTAAACCCACCATGATGGGTTTAATCAcc  >  2:369410/1‑149 (MQ=255)
                                                                                                                               ||                                                                                                                                   
ACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACC  >  NC_000913_3_hsa_pgi/2173234‑2173493

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: