Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913_3_hsa_pgi	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_pgi	1978503	1979278	776	24 [0]	[0] 24	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_hsa_pgi	= 1978502	0 (0.000)	22 (0.670)	21/258	NT	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_hsa_pgi	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_pgi/1979137‑1979420                                                                                                                                               |                                                                                                                                              ggtgctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacaataaacactcaattacttaacataaatgTATAAGt                                                                                                                                         <  1:660341‑M2/7‑1 (MQ=255)     ctgtacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                     <  2:460168‑M2/11‑1 (MQ=255)                 ctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                         >  1:359868‑M2/127‑149 (MQ=255)                    agctgtttcctaaatcgacgcaacggtactcgtcactacgcgcacatacaacggaggggggctgcgattttaaacaatgagtgatgcagatcacacaacacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                      <  2:1446‑M2/26‑1 (MQ=255)                             ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTCAAt                                                                                                             >  1:379650‑M2/115‑149 (MQ=255)                               aaatcgacggcactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                           <  2:655771‑M2/37‑1 (MQ=255)                                       tcagctgtactcgtcactacacgcacattaaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                   <  2:637393‑M2/45‑1 (MQ=255)                                        caactgtactcgtcactacacgcacatactacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                                                                  <  1:572651‑M2/46‑1 (MQ=255)                                            tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                              >  2:575573‑M2/100‑149 (MQ=255)                                                 tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAACGTTTTTccct                                                                                         >  2:120184‑M2/95‑149 (MQ=255)                                                 tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTATTTTGGGTGTATTTCAATCTGTTTAAAAGTTTTTcgct                                                                                         >  2:27710‑M2/95‑149 (MQ=255)                                                       ctacacgcacatacaacggaggggggccacgattttcaataatgcctgatgcagatcacacaaaacactcaattactcaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                   <  2:319698‑M2/61‑1 (MQ=255)                                                               acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                           >  1:265854‑M2/81‑149 (MQ=255)                                                                    caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                      >  1:596840‑M2/76‑149 (MQ=255)                                                                      acggaggggggctgcgattttcaatcatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                                                                    <  1:602782‑M2/76‑1 (MQ=255)                                                                             ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                             >  1:261344‑M2/67‑149 (MQ=255)                                                                                 ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                         >  1:148999‑M2/63‑149 (MQ=255)                                                                                  tgcgattttcaagaaggagtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTAGTAGAACATTTTTAGCTACGCTAGCAAGCATAACTGAAAAACgaaaaatg                                                        >  2:136770‑M2/62‑148 (MQ=255)                                                                                       ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataagtgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                   >  1:438165‑M2/57‑149 (MQ=255)                                                                                                      gatgcagatcacacaaaactctcaattacttaacattaatgTCTAAGTCATACTTTTTTTTTGGGTGGATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                    <  2:369336‑M2/108‑1 (MQ=255)                                                                                                                acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                          >  1:77917‑M2/32‑149 (MQ=255)                                                                                                                          ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCCCCGTaa                >  2:249747‑M2/22‑149 (MQ=255)                                                                                                                                  cttaacataaatgTATAAGTGATACTTTTGTTTTGGGTGTATTTCAATTTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATGGAAATGGGTTGTTGCGCACAGTACACATACCttaccgtaaaaattaat        >  2:289236‑M2/14‑142 (MQ=255)                                                                                                                                        ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGGCAAAATAGCTCACCGTAAATCATCAACTGCta  >  1:217017‑M2/8‑149 (MQ=255)                                                                                                                                               |                                                                                                                                              GTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_pgi/1979137‑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: ‑