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–1979268 1979278 11–776 30 [0] [0] 29 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)25 (0.740) 22/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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913_3_hsa_tpiA/1979133‑1979421
                                                                                                                                                  |                                                                                                                                              
aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggagggaggaaaagacctctgataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                              >  1:327986‑M2/147‑149 (MQ=255)
         tttacacttatctaagatttttcctaaatcgacgtaactgtactcgtaactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagagcacacaaaacaatcaattacttaacataaatgTATAAGTCAtaa                                                                                                                                     >  2:544181‑M2/138‑148 (MQ=255)
                    cgaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:484409‑M2/23‑1 (MQ=255)
                                ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                              >  1:617292‑M2/115‑149 (MQ=255)
                                 taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                             >  1:634480‑M2/114‑149 (MQ=255)
                                               tgtactcgtcactacacgcacatacaaaggaggggggctgcgattgtcaataatgcgtgatgcagatcacacaaaacactcaattacttaaaataaatgTATAAGTCAAACTCTTGTTTTGGGTGTATTTCAATCTGTAAAAAAGcttt                                                                                               >  2:394626‑M2/100‑149 (MQ=255)
                                                             cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                                                                 <  1:660888‑M2/64‑1 (MQ=255)
                                                               cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                               <  1:566712‑M2/66‑1 (MQ=255)
                                                                  acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                            >  1:316721‑M2/81‑149 (MQ=255)
                                                                  acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                            >  1:316854‑M2/81‑149 (MQ=255)
                                                                  acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                            <  2:513742‑M2/69‑1 (MQ=255)
                                                                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                          >  1:233593‑M2/79‑149 (MQ=255)
                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTTCGCTAGCAAGCaaaaa                                                                         <  1:601179‑M2/72‑1 (MQ=255)
                                                                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                         <  1:610068‑M2/72‑1 (MQ=255)
                                                                        aacggaggtgggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                      <  1:106971‑M2/75‑1 (MQ=255)
                                                                                               aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                               >  1:410692‑M2/52‑149 (MQ=255)
                                                                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgGATAGGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATCATCGAAATGGGATGTTg                                           >  1:29227‑M2/48‑149 (MQ=255)
                                                                                                       gtgatgcagatcacataaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                       <  1:481572‑M2/106‑1 (MQ=255)
                                                                                                       gtgatgcacatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                       <  2:46992‑M2/106‑1 (MQ=255)
                                                                                                          atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGGTAAAAAGTTTTTCGCTATGCTAGCAAGCAAAAATGAAACAGGAATGATCGAAATGGGATGTTGCGCACAg                                    >  1:634272‑M2/41‑149 (MQ=255)
                                                                                                                  cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTTTTTTTTTTTTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGTTAGCAAGTAAAAATGAAATAGGAATATTCGAAATGGGATGTTGCACACATAAAAAATa                            >  2:379320‑M2/33‑149 (MQ=255)
                                                                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTTTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                         >  1:364501‑M2/30‑149 (MQ=255)
                                                                                                                       aaaacactgcattacataacacaaatgGATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                       <  1:177498‑M2/122‑1 (MQ=255)
                                                                                                                       aaaacactcaattacttaacataaaggTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGTTTTTCGCTACGCTAGCAAGCAAAAATGCAACAGGAGTAATCGAAATGGGATGTTGCACACAGTCAAAATAACTCa                       >  1:521357‑M2/28‑149 (MQ=255)
                                                                                                                              tcaattaattaacataaatgTATAATTTCTACTTTTGTTTTGGTTCTATTTTAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa                <  1:214565‑M2/129‑1 (MQ=255)
                                                                                                                                 attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataa             <  1:212377‑M2/132‑1 (MQ=255)
                                                                                                                                      ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc        >  1:93999‑M2/13‑149 (MQ=255)
                                                                                                                                           ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAGGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta   >  1:568839‑M2/8‑149 (MQ=255)
                                                                                                                                            taaatgTATAAGTCATACTTTTGTATTGGGTGTATTTCCATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  >  2:176649‑M2/7‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913_3_hsa_tpiA/1979133‑1979421

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ 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.