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

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913_3_hsa_tpiA 1978503–1979269 1979278 10–776 42 [0] [0] 42 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)41 (0.820) 36/268 0.5 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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_hsa_tpiA/1978356‑1978648
                                                                                                                                                  |                                                                                                                                                  
gAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGta                                                                                                                                                  >  2:277736‑M1/1‑147 (MQ=255)
   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  1:606966‑M1/149‑6 (MQ=255)
            tACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactt                                                                                                                                      >  2:712435‑M1/1‑135 (MQ=255)
              cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                    <  1:302878‑M1/149‑17 (MQ=255)
                  tATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttt                                                                                                                                >  2:468190‑M1/1‑129 (MQ=255)
                    tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                              >  1:295549‑M1/1‑127 (MQ=255)
                    tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                              >  1:506557‑M1/1‑127 (MQ=255)
                     ctgagatTTTTCCTAAGGCGGCGCAACTGTGCTCGTCACGACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  1:540020‑M1/146‑24 (MQ=255)
                      tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  2:767719‑M1/1‑125 (MQ=255)
                      tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                            >  1:578173‑M1/1‑125 (MQ=255)
                          aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                        >  1:492694‑M1/1‑121 (MQ=255)
                                 cTAAATCGACGCAATTGTACCCGTCACTAAACGCACATAATACTGAACGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 <  2:416828‑M1/149‑36 (MQ=255)
                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 <  1:840275‑M1/149‑36 (MQ=255)
                                   aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                               <  1:812838‑M1/149‑38 (MQ=255)
                                         aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                         <  1:554051‑M1/149‑44 (MQ=255)
                                         aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaca                                                                                                         >  1:393077‑M1/1‑106 (MQ=255)
                                                  tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                <  1:69243‑M1/149‑53 (MQ=255)
                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                             >  1:133679‑M1/1‑94 (MQ=255)
                                                         cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                         <  2:470099‑M1/149‑60 (MQ=255)
                                                             acacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                     <  1:467714‑M1/149‑64 (MQ=255)
                                                               acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                   <  1:651321‑M1/149‑66 (MQ=255)
                                                               acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                   <  1:468190‑M1/149‑66 (MQ=255)
                                                                cGCACAAACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                                  <  2:104978‑M1/149‑67 (MQ=255)
                                                                        cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                          <  2:13045‑M1/149‑75 (MQ=255)
                                                                                gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                                  <  1:211598‑M1/149‑83 (MQ=255)
                                                                                      tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                                 <  2:425476‑M1/144‑84 (MQ=255)
                                                                                      tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                                 >  1:425476‑M1/1‑61 (MQ=255)
                                                                                                          gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                        <  2:492694‑M1/149‑109 (MQ=255)
                                                                                                                gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                  <  1:388047‑M1/149‑115 (MQ=255)
                                                                                                                 aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  1:268977‑M1/1‑34 (MQ=255)
                                                                                                                  tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaaggggatgttgcgaccagtcaaaat                                >  2:492801‑M1/1‑33 (MQ=255)
                                                                                                                    acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  2:506557‑M1/149‑119 (MQ=255)
                                                                                                                    acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  2:578173‑M1/149‑119 (MQ=255)
                                                                                                                      aaaaaaaactAATTTACCTAACAGAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgttacgctagcaagtaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                            <  2:175622‑M1/147‑121 (MQ=255)
                                                                                                                        aaaaCACTCTATTACTTAACATTAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgatatgggatgttgcgcacagtcaaaataactca                          <  2:508399‑M1/149‑123 (MQ=255)
                                                                                                                              cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                    >  2:27840‑M1/1‑21 (MQ=255)
                                                                                                                               tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                   <  1:27840‑M1/149‑130 (MQ=255)
                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                 >  1:727159‑M1/1‑18 (MQ=255)
                                                                                                                                       ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc           >  2:312868‑M1/1‑12 (MQ=255)
                                                                                                                                       ttAACATAAATGgataaggcatacttttgttttgggtgtatttcaatctgttacaaagtttttccctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc           <  2:98826‑M1/149‑138 (MQ=255)
                                                                                                                                         aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg         <  2:727159‑M1/149‑140 (MQ=255)
                                                                                                                                                aTGtataagtcatacttttgttttgggtgtatttcaatccgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  <  1:712435‑M1/149‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_hsa_tpiA/1978356‑1978648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.