Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_hsa_tpiA 275,065 G→A T18I (ACC→ATC)  insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_hsa_tpiA275,0650GA95.0% 41.7 / ‑6.4 20T18I (ACC→ATC) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/12);  minor base C (0/1);  total (7/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 = 1.00e+00

AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTAA  >  NC_000913_3_hsa_tpiA/275008‑275210
                                                         |                                                                                                                                                 
aaGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCg                                                        >  1:622373/1‑149 (MQ=25)
  gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                      <  2:606404/149‑1 (MQ=37)
    aaaTCTCTTTTCTGGTCTGAAGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTAt                                                    <  2:657899/149‑1 (MQ=255)
       tctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTCTGAg                                                 >  1:810295/1‑149 (MQ=255)
                ttggcttaCGGCGATTACTGCTGACTTCACTGTCGGTGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGcc                                        <  2:831110/142‑1 (MQ=17)
                   tCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCtt                                     >  2:797610/1‑149 (MQ=255)
                     tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTGATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCAtttt                                   <  2:709987/149‑1 (MQ=255)
                     tGACGGCGCCTGCTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCtttt                                   <  1:606804/149‑1 (MQ=255)
                      gACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTCGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGCGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                                  <  2:1039241/149‑1 (MQ=255)
                           cgcTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATc                             >  1:605556/1‑149 (MQ=255)
                           cgcTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTGCTACGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGAATc                             <  2:1049071/149‑1 (MQ=255)
                           cgcTTACTGCTGAATTCACTGTCGGCCAAGATAAGTTGATGACTCATGATGAACCCTGTTCTAAGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATc                             <  2:1030706/149‑1 (MQ=255)
                                aCTGCTGAATTCACTGTCGGCGAAGCTAAGTTGATGACTCATGATGATCCCTGTGCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatt                        <  1:98635/149‑1 (MQ=255)
                                 gtgctgAATTCACGGTCGGCGAAGATAAGTTGATAGCTCTTGATCAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCTCACCTTCCGGAGGCGTTATTAGCTCGCGGCCCTTTTTGTATCTGatta                       <  2:524024/148‑1 (MQ=17)
                                  tgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTTTTCGCACCTTCCGGAGGCATCATTAACTGCCGGCACTGTTTGTGTCTGattat                      >  2:466178/1‑149 (MQ=14)
                                   gctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACATTCCGGAGGCGTACTGAGCTGGTGGGGCTTTATTTATCTGattatt                     <  2:390181/149‑1 (MQ=14)
                                    ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGAttatta                    >  2:304572/1‑149 (MQ=255)
                                            aCTGTCGGCGAAGATAAGTTGATGACTCCTGAAGATCCCTGTTTTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAc            <  1:209487/149‑1 (MQ=21)
                                                   gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAt     <  1:731091/149‑1 (MQ=255)
                                                      aaGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTaa  >  1:654420/1‑149 (MQ=255)
                                                         |                                                                                                                                                 
AAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTATTAA  >  NC_000913_3_hsa_tpiA/275008‑275210

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