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,0650GA100.0% 69.5 / NA 27T18I (ACC→ATC) insH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/13);  total (14/13)

GGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_hsa_tpiA/275003‑275207
                                                              |                                                                                                                                              
gggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACGGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCgg                                                          <  2:530332/149‑1 (MQ=14)
 ggACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGa                                                         <  2:62943/149‑1 (MQ=17)
       gAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGtt                                                   >  1:616849/1‑149 (MQ=37)
        aaaaTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTa                                                  >  1:945727/1‑149 (MQ=38)
           aTCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGCTGAACCCTGTTCTATGGCTCCAGCTGACAAACATGATCTCATATCAGAGACATGTTCGCACCTTCCGGAGGCGTTATGa                                               >  2:147794/1‑149 (MQ=14)
            tctcTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAg                                              <  1:359924/149‑1 (MQ=255)
             ctctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCCGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAgg                                             >  1:275767/1‑148 (MQ=255)
                  ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcg                                        >  2:1000539/1‑149 (MQ=255)
                    cTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATTGGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTggcggc                                      >  2:701883/1‑149 (MQ=255)
                      ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGccc                                    <  1:911827/149‑1 (MQ=255)
                            aCGGCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGGGGCGTTATGAGCTGGCGGCCCTTTTTg                              >  1:422110/1‑149 (MQ=255)
                             tggCGCTTACTGGTGAATTCACTGTCGGGGAAGATAAGTTGATTGCTTATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGTACCTTCCGGAGGCGTTTTTAGCTGGCGGCCCTTTTTGt                             <  2:215086/148‑1 (MQ=17)
                             cggtgtttACTGCTGAATTCACTGTCGGCGAACATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCATATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGt                             <  2:94981/143‑1 (MQ=17)
                              ggCGCTTACTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTCTGAGCTGGCGGCCCTTTTTGTa                            >  2:420657/1‑149 (MQ=255)
                                     aCTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatt                     <  1:741448/149‑1 (MQ=255)
                                     aCTGCTGAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGatt                     <  2:921201/149‑1 (MQ=255)
                                       tgctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGattat                   <  2:79201/149‑1 (MQ=255)
                                         ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGAttatta                 >  2:585304/1‑149 (MQ=255)
                                         ctgAATTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCATTTTTGGATCTGAttatta                 >  2:1028090/1‑149 (MQ=255)
                                            aaTTGACTGTCGGCGAAGATAAGATGATGACTCATGATGAACCCTATTCTATGGCTCCAGATCACAAACATGATCTCATATCAGGGACTTGCTCGCACATTCCGTAGGCGTTATGAGCTCGCGGCCCTTTTTGTATCTGATTATTAATc              <  1:886769/149‑1 (MQ=18)
                                            aaTTCACTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATc              <  2:616849/149‑1 (MQ=255)
                                               tCATTGTCGGCGAAGATAAGGTTGTGTCTCAGGATGAATCCTGTTCTATGGCTCCATATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATcccc           <  1:484913/149‑1 (MQ=21)
                                                 aCTGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGCCTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTTTCTGACAATTCATCCCCAc         >  2:350910/1‑149 (MQ=18)
                                                   tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAccc       >  1:691156/1‑149 (MQ=255)
                                                   tGTCGGCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAccc       >  1:685214/1‑149 (MQ=255)
                                                   ggTCGGAGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACAGGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCAccc       <  1:387851/148‑1 (MQ=255)
                                                        gCGAAGATAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATATCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCAGCTAt  >  2:8118/1‑149 (MQ=255)
                                                              |                                                                                                                                              
GGGACAAGAAAATCTCTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCGGAGGCGTTATGAGCTGGCGGCCCTTTTTGTATCTGATTATTAATCCCCACCCGCTAT  >  NC_000913_3_hsa_tpiA/275003‑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: