Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 2,679,525 G→T A1074D (GCT→GAT)  yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,679,5250GT97.0% 78.2 / ‑6.6 33A1074D (GCT→GAT) yphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (14/18);  minor base A (0/1);  total (14/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 = 9.70e-01

TAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGT  >  NC_000913_3_hsa_tpiA/2679379‑2679668
                                                                                                                                                  |                                                                                                                                               
tAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTATTAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTAGGCCAATcc                                                                                                                                               >  2:67446/1‑149 (MQ=255)
  gggTTTTTATTAACAACGGCTTATTCTAATTATTTGGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCgg                                                                                                                                             >  2:530102/1‑149 (MQ=255)
     tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGATGGTGCATAATGCCGCTTTCTGCCCATCCGGAtt                                                                                                                                          >  1:73332/1‑149 (MQ=255)
        ttatttACAAGGTCTTATTATAATTATTGTGTGAGGGGCTGCAGCGCGTGCTTTTAAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAAt                                                                                                                                       <  2:416447/149‑1 (MQ=255)
            taACAACGGCTTAGTGTAATTATTTTGTGATGAGCGGTAGCGCGTGTCGATTAAATAATACAGTTAATAGGGTTTAGTGAACGTGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGcc                                                                                                                                   <  2:544527/149‑1 (MQ=255)
                     cTTATTCTTGTTTGTTTGTGTTGTGCGGCAGCGCGTGCCGCTTAAATATTATGGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCa                                                                                                                          <  2:322585/149‑1 (MQ=255)
                         ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAgcgc                                                                                                                      >  1:75220/1‑149 (MQ=255)
                           cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGa                                                                                                                    <  1:373056/149‑1 (MQ=255)
                                 attTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGAc                                                                                                              >  1:378804/1‑149 (MQ=255)
                                   tttGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTc                                                                                                            <  1:63928/149‑1 (MQ=255)
                                          tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                                                                                                     >  1:399976/1‑149 (MQ=255)
                                             gCGGCAGCGCGTGACGCTCAAATATTACAGTTCCGAGGGCTTAGTGAACGCGGTTAAAAATAAAGGGCATCCCGGTGGTGAATAACGCCGCCTCCGTCCAGTCCGGGTTCATCGACGTCAACATAACGCGAACAGAGTCATCTTgcgaa                                                                                                  >  2:291429/1‑147 (MQ=255)
                                                gCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAcc                                                                                               >  2:1853/1‑149 (MQ=255)
                                                       gTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCAAATAAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGa                                                                                        <  1:271366/149‑1 (MQ=255)
                                                             cgcAAATATGACAGTTAATAGGTTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTCCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCACACTCATATTGCGCAACCTCCCCGAGTCCcg                                                                                  <  1:479245/147‑1 (MQ=255)
                                                                aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCgg                                                                               >  1:113883/1‑149 (MQ=255)
                                                                     ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAcccc                                                                          >  1:261984/1‑149 (MQ=255)
                                                                     ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc                                                                          <  1:53426/149‑1 (MQ=255)
                                                                       aCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTCTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCCTCATACGACACCTCCCCGAGTCCCGCCGCCGCCac                                                                        >  2:232891/1‑148 (MQ=255)
                                                                                              gctgTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAAACCCATTGGGCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTTCCGCGGACGCCAGCATCGCCACCTTCAGGCAttttt                                                 <  1:444682/146‑1 (MQ=255)
                                                                                              gcgGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCAttttt                                                 <  1:530102/149‑1 (MQ=255)
                                                                                                ggTGAAAAATAAGAGGCATCACTGTGGTGAATAATGCCGCTTTCGGCCATTCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCt                                               <  2:331434/149‑1 (MQ=255)
                                                                                                 gTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCtt                                              <  2:75220/149‑1 (MQ=255)
                                                                                                 gTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATACGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCtt                                              <  2:378804/149‑1 (MQ=255)
                                                                                                           aaaGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGctgctg                                    <  1:66034/149‑1 (MQ=255)
                                                                                                           aaaGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGctgctg                                    >  2:66034/1‑149 (MQ=255)
                                                                                                                        tggtgAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCa                       <  1:553832/149‑1 (MQ=255)
                                                                                                                            gAATAATGCCGCTTTCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACa                   >  1:349973/1‑149 (MQ=255)
                                                                                                                                    cccttttCGGCGAATCCGGATTAATCCCCGTAACTTCAGCGCGAGCAGTCTAATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGATGACGTAATTCGCCATACCGCGacaac           <  2:192581/145‑1 (MQ=255)
                                                                                                                                          tCGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCa     <  2:472824/149‑1 (MQ=255)
                                                                                                                                           cGGCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAg    <  2:499371/149‑1 (MQ=255)
                                                                                                                                           cGGCAAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAg    <  2:349973/149‑1 (MQ=255)
                                                                                                                                             gCCAATCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCCGGt  >  2:53214/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAGCATCGCCACCATCAGGCATTTTTCTTTATGCTGCTGTTGTAAATCGCCATACAGCGACAACAAGTCAGGT  >  NC_000913_3_hsa_tpiA/2679379‑2679668

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