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	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	= 1978502	0 (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

CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_tpiA/1979139‑1979420                                                                                                                                             |                                                                                                                                               gtctctacacttatctaacatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                                                                        <  1:358386‑M2/9‑1 (MQ=255)     ttacacttatctaagatttttcctaaatctacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatccgtgattcaaatcacacaaaacacctaattacttaacataaatgTATAAGTCATACt                                                                                                                                    >  2:215286‑M2/137‑149 (MQ=255)             atctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                            <  1:427456‑M2/21‑1 (MQ=255)               ctaagatttttccaaaatcgacgcatctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:651805‑M2/23‑1 (MQ=255)                taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTGTggg                                                                                                                         >  1:490383‑M2/126‑149 (MQ=255)                    atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                     >  2:815342‑M2/122‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAAt                                                                                                             >  1:345411‑M2/115‑149 (MQ=255)                           ctaaatcgacgcaactgtactcgtcactacacgcacatacaaaggaggggggctgagattttcaataatgcgtgatgcagatcacacaaaacaataaattccttaacataaatgTATAAGTCATAATTTTGTTTTGGGTGTCT‑TTAAAt                                                                                                             >  2:184955‑M2/115‑149 (MQ=255)                              aattgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATACGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTg                                                                                                          <  2:366800‑M2/38‑1 (MQ=255)                                tcgacgcaactgtactcgtcactacacgcacatacaacggaggagggctgcgattttcaataatgcgtgatgcggatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGGTTTGGGTGTAT‑TTCAATCTGtt                                                                                                        >  1:201586‑M2/110‑149 (MQ=255)                                     gcaactgttctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacataaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTaaaaa                                                                                                   <  1:746166‑M2/45‑1 (MQ=255)                                        actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAgtc                                                                                                >  2:276963‑M2/102‑148 (MQ=255)                                          tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGtttt                                                                                              >  1:788733‑M2/100‑149 (MQ=255)                                           gtactcgttaatacacggacatacaacggtggggggctgcgatttttaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGttttt                                                                                             <  2:235925‑M2/51‑1 (MQ=255)                                           gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGttttt                                                                                             <  1:495797‑M2/51‑1 (MQ=255)                                           gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGttttt                                                                                             >  1:679437‑M2/99‑149 (MQ=255)                                            tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTc                                                                                            <  1:309466‑M2/52‑1 (MQ=255)                                               tcgtcactacacgcacatacaacggaggggggctgcgattttcactaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTcgct                                                                                         >  2:33777‑M2/95‑149 (MQ=255)                                               tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacaccaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAATGTtttcacg                                                                                         >  1:59598‑M2/95‑146 (MQ=255)                                               tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTcgct                                                                                         >  2:534806‑M2/95‑149 (MQ=255)                                                cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTcgcta                                                                                        <  1:201915‑M2/56‑1 (MQ=255)                                                      tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTtcgctactctac                                                                                  >  2:380264‑M2/88‑148 (MQ=255)                                                                tacaacggaggggggctgcgatgttcaataatgcgtcatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTCCGCTACGCTAGCAAGCaaaaa                                                                        <  2:262345‑M2/72‑1 (MQ=255)                                                                 acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattgcttaacataaatgTATACGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                                                                       >  1:327737‑M2/77‑149 (MQ=255)                                                                   aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTTTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                     <  2:561070‑M2/75‑1 (MQ=255)                                                                   aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                     >  1:657776‑M2/75‑149 (MQ=255)                                                                    acggctgggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                    >  2:312045‑M2/74‑149 (MQ=255)                                                                        aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                >  1:321791‑M2/70‑149 (MQ=255)                                                                          gggggctgctattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGAAAAAATGAAACAGGa                                                              >  2:814861‑M2/68‑149 (MQ=255)                                                                                             aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTCTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCTAAATGgcataca                                           >  2:838400‑M2/49‑146 (MQ=255)                                                                                               tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTGTAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                         >  2:788785‑M2/47‑149 (MQ=255)                                                                                               tgcgtaatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgc                                         <  2:528605‑M2/103‑1 (MQ=255)                                                                                                gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                        >  2:605975‑M2/46‑149 (MQ=255)                                                                                                        cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                >  1:471751‑M2/38‑149 (MQ=255)                                                                                                          gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCTCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACGGTCaaa                              >  2:190527‑M2/36‑149 (MQ=255)                                                                                                           atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                             >  1:634904‑M2/35‑149 (MQ=255)                                                                                                                   aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                     <  2:56262‑M2/123‑1 (MQ=255)                                                                                                                    aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc                    >  2:658118‑M2/26‑149 (MQ=255)                                                                                                                           aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata             >  1:631591‑M2/19‑149 (MQ=255)                                                                                                                              tacttaacATAAATGTATAAGTAATACTTTTCGGGTCGGTGTATCATGTCGATG‑TAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc          <  1:397352/142‑1 (MQ=255)                                                                                                                                 ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc       <  1:215286‑M2/137‑1 (MQ=255)                                                                                                                                   aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTAACCGTAAATAATCATCTg     >  1:262117‑M2/11‑149 (MQ=255)                                                                                                                                     cataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAACAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt   >  2:215725‑M2/9‑149 (MQ=255)                                                                                                                                     cataaacgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCt   <  1:17441‑M2/141‑1 (MQ=255)                                                                                                                                      ataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta  >  1:519171‑M2/8‑149 (MQ=255)                                                                                                                                      ataaatgTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGCTAAAAAGTTTTTCGCTACGCTAGCAAGCACAAATGAAACAGCAATAATCGAAATGGGAAGTTGCACACAGTCAAAATACCTCAACGTaaataatcatcctcta  >  2:682799‑M2/8‑144 (MQ=255)                                                                                                                                             |                                                                                                                                               CCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTAT‑TTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_tpiA/1979139‑1979420

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