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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_hsa_tpiA	1978503–1979267	1979278	12–776	30 [0]	[0] 30	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	= 1978502	0 (0.000)	29 (0.950)	21/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

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_hsa_tpiA/1979139‑1979424                                                                                                                                             |                                                                                                                                                  gcctttacacttatctaagatttttcctaaaccgagtcaactgcactcggcacgacactcacatacaacggaggggggctgcgattttcaataatgcgtgctgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                                                                           >  2:456479‑M2/141‑149 (MQ=255)    tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                        >  1:10767‑M2/138‑149 (MQ=255)    tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAAAc                                                                                                                                        >  1:363195‑M2/138‑149 (MQ=255)    tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcaaacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                        >  1:37027‑M2/138‑149 (MQ=255)          cttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                                  <  1:591813‑M2/18‑1 (MQ=255)           ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                                 >  2:148294‑M2/131‑149 (MQ=255)             atctaaaatctttcctaaatcgacgcaactatactcgtcactacacccacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaccctcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                               <  1:49102‑M2/21‑1 (MQ=255)              cctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatttccaataatgcgtgatgcagaccacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                              <  1:157933‑M2/22‑1 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  1:428479‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  2:587258‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaataTATAATTCATACTTTTTTTTTTGGTGTATTTCAAt                                                                                                                 >  2:446450‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatttccaatcatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  2:562371‑M2/115‑149 (MQ=255)                            taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                <  2:340536‑M2/36‑1 (MQ=255)                                        actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                    >  2:6805‑M2/102‑149 (MQ=255)                                               tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcact                                                                                             >  1:342326‑M2/95‑149 (MQ=255)                                                    actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTTATTCTGTTAAAATGTTTTTtgctacgct                                                                                        >  1:315977‑M2/90‑149 (MQ=255)                                                     ctacacgcacatacaacggaggggggctgcgattttcaataatgcgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacccta                                                                                       >  2:179721‑M2/89‑149 (MQ=255)                                                         acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCGGTTAAAAAGTTTTTCGCTACGCTagaaa                                                                                   >  1:15194‑M2/85‑149 (MQ=255)                                                                       gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                     <  2:612407‑M2/79‑1 (MQ=255)                                                                               ctgggattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                             <  1:289876‑M2/87‑1 (MQ=255)                                                                                tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCCCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                            >  2:295733‑M2/62‑149 (MQ=255)                                                                                  ctgttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa                                                          <  1:472445‑M2/90‑1 (MQ=255)                                                                                          gataatgcgtgaggcagatctcacaaaacattcaattacttaccgtaaattTATAAGTCATTCTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTGTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                  <  2:149432‑M2/98‑1 (MQ=255)                                                                                          aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                                  <  1:322987‑M2/98‑1 (MQ=255)                                                                                           ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                                 >  1:550291‑M2/51‑149 (MQ=255)                                                                                                 cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTTTTCTTTCTCAAACAAAAATGAAACACGAATACTCGACATGGtatgtttagc                                           >  2:75056‑M2/45‑145 (MQ=255)                                                                                                           atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                 >  1:5175‑M2/35‑149 (MQ=255)                                                                                                            tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCTAAATGGGATGTTGCGCACAGTAAAAAt                                >  2:553088‑M2/34‑149 (MQ=255)                                                                                                            tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTAAAAAt                                >  1:101316‑M2/34‑149 (MQ=255)                                                                                                                                          atgTAAAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  <  2:41929‑M2/146‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913_3_hsa_tpiA/1979139‑1979424

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: ‑