Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC NC_000913_3_hsa_tpiA 1,978,503 Δ776 bp 100% insB1insA insB1, insA

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_hsa_tpiA 1978503 1979278 776 46 [1] [1] 46 insB1–insA insB1,insA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 19785020 (0.000)39 (0.930) 34/258 NT 100% intergenic (‑305/+16) flhD/insB1 flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
?NC_000913_3_hsa_tpiA 1979279 = 0 (0.000)intergenic (‑64/‑474) insA/uspC IS1 repressor TnpA/universal stress protein

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_hsa_tpiA/1978363‑1978646
                                                                                                                                           |                                                                                                                                                
gTCTCTACACTTATCTAACATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                         <  1:358386‑M1/149‑10 (MQ=255)
    ttACACTTATCTAAGATTTTTCCTAAATCTACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATCCGTGATTCAAATCACACAAAACACCTAATTACTTAACATAAATGtataagtcatact                                                                                                                                     >  2:215286‑M1/1‑136 (MQ=255)
            aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                             <  1:427456‑M1/149‑22 (MQ=255)
              cTAAGATTTTTCCAAAATCGACGCATCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                           <  2:651805‑M1/149‑24 (MQ=255)
               tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttgtggg                                                                                                                          >  1:490383‑M1/1‑125 (MQ=255)
                   aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                      >  2:815342‑M1/1‑121 (MQ=255)
                          cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                               >  1:345411‑M1/1‑114 (MQ=255)
                          cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAAAGGAGGGGGGCTGAGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATAAATTCCTTAACATAAATGtataagtcataattttgttttgggtgtctttaaat                                                                                                               >  2:184955‑M1/1‑114 (MQ=255)
                             aaTTGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatacgtcatacttttgttttgggtgtatttcaatctg                                                                                                            <  2:366800‑M1/149‑39 (MQ=255)
                               tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCGGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtatttcaatctgtt                                                                                                          >  1:201586‑M1/1‑109 (MQ=255)
                                    gCAACTGTTCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                     <  1:746166‑M1/149‑46 (MQ=255)
                                       aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtc                                                                                                  >  2:276963‑M1/1‑101 (MQ=255)
                                         tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                >  1:788733‑M1/1‑99 (MQ=255)
                                          gTACTCGTTAATACACGGACATACAACGGTGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               <  2:235925‑M1/149‑52 (MQ=255)
                                          gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               <  1:495797‑M1/149‑52 (MQ=255)
                                          gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               >  1:679437‑M1/1‑98 (MQ=255)
                                           tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                              <  1:309466‑M1/149‑53 (MQ=255)
                                              tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCACTAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                           >  2:33777‑M1/1‑94 (MQ=255)
                                              tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACCAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaatgttttcacg                                                                                           >  1:59598‑M1/1‑94 (MQ=255)
                                              tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                           >  2:534806‑M1/1‑94 (MQ=255)
                                               cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                          <  1:201915‑M1/149‑57 (MQ=255)
                                                     tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctactctac                                                                                    >  2:380264‑M1/1‑87 (MQ=255)
                                                               tACAACGGAGGGGGGCTGCGATGTTCAATAATGCGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttccgctacgctagcaagcaaaaa                                                                          <  2:262345‑M1/149‑73 (MQ=255)
                                                                aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtatacgtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                         >  1:327737‑M1/1‑76 (MQ=255)
                                                                  aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcttttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       <  2:561070‑M1/149‑76 (MQ=255)
                                                                  aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       >  1:657776‑M1/1‑74 (MQ=255)
                                                                   acggctgggggctgcgaTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                      >  2:312045‑M1/7‑73 (MQ=255)
                                                                       aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                  >  1:321791‑M1/1‑69 (MQ=255)
                                                                         gggggCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacagga                                                                >  2:814861‑M1/1‑67 (MQ=255)
                                                                                            aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttcaatctgttaaaaagttcttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatggcataca                                             >  2:838400‑M1/1‑48 (MQ=255)
                                                                                              tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtgtaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                           >  2:788785‑M1/1‑46 (MQ=255)
                                                                                              tGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                           <  2:528605‑M1/149‑104 (MQ=255)
                                                                                               gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                          >  2:605975‑M1/1‑45 (MQ=255)
                                                                                                       cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                                  >  1:471751‑M1/1‑37 (MQ=255)
                                                                                                         gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacggtcaaa                                >  2:190527‑M1/1‑35 (MQ=255)
                                                                                                          aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                               >  1:634904‑M1/1‑34 (MQ=255)
                                                                                                                  aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                       <  2:56262‑M1/149‑124 (MQ=255)
                                                                                                                   aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                      >  2:658118‑M1/1‑25 (MQ=255)
                                                                                                                          aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata               >  1:631591‑M1/1‑18 (MQ=255)
                                                                                                                                ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc         <  1:215286‑M1/149‑138 (MQ=255)
                                                                                                                                  aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaataatcatctg       >  1:262117‑M1/1‑10 (MQ=255)
                                                                                                                                    cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     >  2:215725‑M1/1‑8 (MQ=255)
                                                                                                                                    cATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     <  1:17441‑M1/149‑142 (MQ=255)
                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta    >  1:519171‑M1/1‑7 (MQ=255)
                                                                                                                                     aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgctaaaaagtttttcgctacgctagcaagcacaaatgaaacagcaataatcgaaatgggaagttgcacacagtcaaaatacctcaacgtaaataatcatcctcta    >  2:682799‑M1/1‑7 (MQ=255)
                                                                                                                                       attggggTAATGACTCCAACTTATTGACAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTTTGCCGCTc  >  1:177056/5‑149 (MQ=2)
                                                                                                                                           |                                                                                                                                                
GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_hsa_tpiA/1978363‑1978646

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: