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 36 [0] [0] 36 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)29 (0.890) 24/260 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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_hsa_tpiA/1978356‑1978648
                                                                                                                                                  |                                                                                                                                                  
aaaaTGGGTTCTAACACTTATCCAAGATCTTACCTAAATCAACGCAACTGTACTCGTCACTACACGCACATACAACGTAGCGGTCCTCCGATTATCAACAATGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATgta                                                                                                                                                  <  2:2950/148‑3 (MQ=255)
 aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 >  1:28094‑M1/1‑146 (MQ=255)
 aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 >  1:28356‑M1/1‑146 (MQ=255)
 aaaTGGGTCTTTACACTTATCTAAGATTTTCCCTAAATCGACGCAACTATACTCGTAACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 <  2:235952‑M1/149‑4 (MQ=255)
   aTGGTTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  2:383770‑M1/149‑6 (MQ=255)
             atccttaTCTAACATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                                     <  1:221461‑M1/146‑16 (MQ=255)
              cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGCTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatactttt                                                                                                                                    >  1:248763‑M1/1‑133 (MQ=255)
                  tacccaagatttCTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGATTTCCAACAATGCGTGATGCAGACCACACAAAACACTCAATTACTTAACATAAATGtataaggcacacccttgtct                                                                                                                                <  2:335899‑M1/144‑21 (MQ=255)
                      tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGATTTTCAATAATGCGTGATGCAGATCACACACAACAGTCAGTTACTTAACATAAATGtataagtcatacttttgttttgga                                                                                                                            >  1:618041‑M1/1‑125 (MQ=255)
                          aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                        >  1:239006‑M1/1‑121 (MQ=255)
                           tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                       <  2:452132‑M1/149‑30 (MQ=255)
                                 cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCACTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:565408‑M1/1‑114 (MQ=255)
                                  ttaatcGACGCTACTGTATTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                <  1:509832‑M1/147‑37 (MQ=255)
                                          cGCAACTGTACTCGTCACTACACGCACATACAACGGAGCGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                        <  2:384398‑M1/149‑45 (MQ=255)
                                           gCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCTATTTTCAATAATGCGTGATGCATATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaca                                                                                                       <  2:439556‑M1/149‑46 (MQ=255)
                                              aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                    >  1:467359‑M1/1‑101 (MQ=255)
                                                tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  >  1:12187‑M1/1‑99 (MQ=255)
                                                        tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgattttcttttattacagtctgttaaaaccttttttgctaca                                                                                          >  2:519016‑M1/1‑91 (MQ=255)
                                                           cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                       <  1:576387‑M1/149‑62 (MQ=255)
                                                             acacGCAAATACAACGGAGGGAGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                     <  2:339361‑M1/149‑64 (MQ=255)
                                                                            ggAGGGGCGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                                                                      <  2:643257‑M1/149‑79 (MQ=255)
                                                                                   ggCTGCGATTTTCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                               <  1:352169‑M1/149‑86 (MQ=255)
                                                                                     cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                             <  2:86489‑M1/149‑88 (MQ=255)
                                                                                           ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                       >  1:384221‑M1/1‑56 (MQ=255)
                                                                                                  taatGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatagcaatctgttaaaaagtgtttcgctacgctagcaagcaaaaatgaaacagcaataatcgaaatgggatgt                                                >  1:460309‑M1/1‑49 (MQ=255)
                                                                                                        gTGATGCAGATCACTCAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttacaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                          <  2:137480‑M1/149‑107 (MQ=255)
                                                                                                          gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                        <  1:591365‑M1/149‑109 (MQ=255)
                                                                                                                    acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:24650‑M1/1‑31 (MQ=255)
                                                                                                                    acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattacaatctgtgaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:364532‑M1/1‑31 (MQ=255)
                                                                                                                        aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                          >  1:352903‑M1/1‑27 (MQ=255)
                                                                                                                        aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                          >  1:307761‑M1/1‑27 (MQ=255)
                                                                                                                           acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                       >  1:528605‑M1/1‑24 (MQ=255)
                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttagaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaataggatgttgcgcacagtcaaaataactacccgtaaaca                 >  2:108087‑M1/1‑18 (MQ=255)
                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaacta                 >  1:108001‑M1/1‑18 (MQ=255)
                                                                                                                                          aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagtaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc        <  1:472324‑M1/149‑141 (MQ=255)
                                                                                                                                                aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaccaggattaatggaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:483490‑M1/1‑3 (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 < 13 ≤ ATCG/ATCG < 32 ≤ 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.