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.5% 63.8 / ‑6.8 22T18I (ACC→ATC) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (10/11);  total (10/12)
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

GAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_hsa_tpiA/275010‑275207
                                                       |                                                                                                                                              
gAAAATCTCTTTTCTGGTCTGACGGCGTTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                   <  1:209964/149‑1 (MQ=255)
gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGAGGACAAACATGATCTCATATCAGGGACTTGTTCGCCCCTTCCGGAGGCGta                                                   >  1:529315/1‑148 (MQ=255)
        ctTTTCTGGCCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCCCATATCCGGGACATCTTCGCTCCTTCCGGCAGCGTTAtcgtctg                                           >  2:505561/1‑143 (MQ=11)
         ttttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTAAGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTgg                                          <  2:124749/149‑1 (MQ=255)
           ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcg                                        <  2:413564/149‑1 (MQ=255)
                gTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCt                                   >  2:117559/1‑149 (MQ=255)
                  cTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCTCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttt                                 <  2:286413/149‑1 (MQ=255)
                  cTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttt                                 <  1:216375/149‑1 (MQ=255)
                    gccGGCGCTTCCTTCTGGATTCACTGTCGGCGAAGATAACTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAATCATGATCTCATATCCGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCttttt                               <  2:756816/147‑1 (MQ=17)
                    gACGGCGCTTACTGCTGCATTCGCAGTAGGCGAAGATAAGTTGATTACTCATGACGCATCCTGTTCTATGACTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGTCCCttttt                               <  2:927763/149‑1 (MQ=17)
                      cGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGt                             <  1:757876/149‑1 (MQ=255)
                             tACGGCTGAATTCACTGTCGGGGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGat                      <  2:606668/149‑1 (MQ=255)
                                tgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAACCATGATCTCATATCAGGGACTTGTTTGCAACTTCCGGCGGCGTTATGAGCTGGCGGCCCTTTTTGTATGTGatttt                   >  1:776419/1‑147 (MQ=17)
                                tgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGattat                   >  1:795253/1‑149 (MQ=255)
                                  ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTa                                                  <  1:464204/116‑1 (MQ=255)
                                  ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTa                                                  >  2:464204/1‑116 (MQ=255)
                                  ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGAttatta                 >  1:439890/1‑149 (MQ=255)
                                            tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAccc       >  2:420336/1‑149 (MQ=255)
                                            tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACATGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAccc       >  2:948193/1‑149 (MQ=255)
                                              tCGGGGAAGGTAAGTGGGTGTCTCATGATGACCCCTGTTCTATGGCTCCAGATGACCAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGc     <  2:1021158/149‑1 (MQ=21)
                                              tCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGc     <  1:1022015/149‑1 (MQ=255)
                                                 gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGCTTATTCATCCCCACCCGCTAt  >  1:186258/1‑149 (MQ=255)
                                                       |                                                                                                                                              
GAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_hsa_tpiA/275010‑275207

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: