Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_hsa_tpiA 2,032,046 T→C 9.8% D91G (GAT→GGT)  dcm ← DNA cytosine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA2,032,0460TC 9.8% 121.8 / 4.4 43D91G (GAT→GGT) dcmDNA cytosine methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (24/14);  new base C (2/2);  total (27/16)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

TATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCA  >  NC_000913_3_hsa_tpiA/2031897‑2032195
                                                                                                                                                     |                                                                                                                                                     
tatCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGAt                                                                                                                                                       <  2:94091/150‑1 (MQ=255)
  tCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAGTAGATCg                                                                                                                                                     >  3:348721/1‑150 (MQ=255)
    ttCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCGATTGTCTGTTAAACGCGCGCTGTCGTCCGATCGGTGCAAAACCGCGACGGATGCCGCCAATTCCGGCGGATAGATCGAt                                                                                                                                                   >  3:35724/1‑150 (MQ=255)
      cATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATaa                                                                                                                                                 >  1:246669/1‑150 (MQ=255)
         tAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCCCTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGc                                                                                                                                              <  4:154106/150‑1 (MQ=255)
               gatgCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACg                                                                                                                                        >  1:156528/1‑150 (MQ=255)
               gatgCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACg                                                                                                                                        <  2:42437/150‑1 (MQ=255)
                                      tGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCggt                                                                                                                 <  3:162567/150‑1 (MQ=255)
                                             cTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTgg                                                                                                          <  1:277494/150‑1 (MQ=255)
                                             cTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTgg                                                                                                          <  3:148249/150‑1 (MQ=255)
                                              tttATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGGACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGGTCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTggg                                                                                                         <  2:154984/150‑1 (MQ=255)
                                               ttATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGt                                                                                                        <  1:1945/150‑1 (MQ=255)
                                                 ataAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGCCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGtaa                                                                                                      >  2:317728/1‑150 (MQ=255)
                                                     gTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAAc                                                                                                  <  3:241887/150‑1 (MQ=255)
                                                      tGCGTACCGCATGTTTATTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACg                                                                                                 >  2:125958/1‑150 (MQ=255)
                                                         gTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGttt                                                                                              <  4:348721/150‑1 (MQ=255)
                                                            ccGCATGTTTGTTCCATTCGCTGGTAACCACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGACCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCa                                                                                           >  2:147058/1‑150 (MQ=255)
                                                            ccGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGACCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCa                                                                                           >  2:336752/1‑150 (MQ=255)
                                                                 tgtttgttCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATgg                                                                                      <  4:119073/150‑1 (MQ=255)
                                                                   tttgttCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATgggg                                                                                    >  1:441928/1‑149 (MQ=255)
                                                                        tCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAAct                                                                               >  2:210284/1‑150 (MQ=255)
                                                                             tCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTtc                                                                          >  3:100396/1‑150 (MQ=255)
                                                                             tCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTtc                                                                          >  3:118460/1‑150 (MQ=255)
                                                                                 tGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCAc                                                                      >  4:382859/1‑150 (MQ=255)
                                                                                         acGCACTGTCCGCCAATCGATTCAAAACCGCGAGGGATGCCGCCAATTCCGGCGAATAGAACGATAAAGCGAAACGCATAAAGCGGATGAAGTGCCGGTGGTATTGGTAATAACGTTTGCAGATGGGCGAACACTTTCTCACTTAAAAgg                                                              >  1:80425/1‑150 (MQ=255)
                                                                                            cACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGc                                                           >  1:565504/1‑150 (MQ=255)
                                                                                               tGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCAt                                                        >  2:569217/1‑150 (MQ=255)
                                                                                                    gCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGa                                                   >  4:307730/1‑150 (MQ=255)
                                                                                                         tCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCAt                                              >  1:232032/1‑150 (MQ=255)
                                                                                                              tCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGcc                                         >  1:244133/1‑150 (MQ=255)
                                                                                                              tCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGcc                                         >  1:244306/1‑150 (MQ=255)
                                                                                                                     cgcgACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCAc                                  >  1:224747/1‑150 (MQ=255)
                                                                                                                     cgcgACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCAc                                  >  3:227992/1‑150 (MQ=255)
                                                                                                                                      aTTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCcgcg                 <  1:317728/150‑1 (MQ=255)
                                                                                                                                      aTTCCGGCGAATAGACCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCcgcg                 <  1:147058/150‑1 (MQ=255)
                                                                                                                                      aTTCCGGCGAATAGACCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCcgcg                 <  1:336752/150‑1 (MQ=255)
                                                                                                                                        tCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCt               <  2:156528/150‑1 (MQ=255)
                                                                                                                                        tCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCt               >  1:126541/1‑150 (MQ=255)
                                                                                                                                            gCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAg           >  2:495898/1‑150 (MQ=255)
                                                                                                                                            gCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAg           <  1:125958/150‑1 (MQ=255)
                                                                                                                                              gAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTg         >  1:40457/1‑150 (MQ=255)
                                                                                                                                                 tAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGAtt      <  4:202650/150‑1 (MQ=255)
                                                                                                                                                    aTCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATtctc   >  2:11394/1‑150 (MQ=255)
                                                                                                                                                    aTCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATtctc   >  2:399322/1‑150 (MQ=255)
                                                                                                                                                     tCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCa  <  1:297830/150‑1 (MQ=255)
                                                                                                                                                     tCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCa  >  3:140958/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATCTTCATTAAAATGATGCGTCGCCGGATCGCAATAATGGTTGGCTTTATAAGTGCGTACCGCATGTTTGTTCCATTCGCTGGTAAACACGCACTGTCCGCCAATCGATTCAAAACCGCGACGGATGCCGCCAATTCCGGCGAATAGATCGATAAAGCGAAACGCATAATGCGGATGATGTGCCGGTGGTTTGGGTAATAACGTTTGCAGATGGGCGAACTCTTTCTCACTTAAACGGTGCCATGCCGAGTCATTCGCCAGCGCACGTTTTAAAATTGCCGCGCTCCAGTGATTCTCA  >  NC_000913_3_hsa_tpiA/2031897‑2032195

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