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

GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_hsa_pgi/1978361‑1978644                                                                                                                                              |                                                                                                                                               ggtgctttaCACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAATAAACACTCAATTACTTAACATAAATGtataagt                                                                                                                                         <  1:660341‑M1/145‑8 (MQ=255)     ctgtacaCTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                     <  2:460168‑M1/146‑12 (MQ=255)                 cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         >  1:359868‑M1/1‑126 (MQ=255)                    agctgtttcctaaATCGACGCAACGGTACTCGTCACTACGCGCACATACAACGGAGGGGGGCTGCGATTTTAAACAATGAGTGATGCAGATCACACAACACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtg                                                                                                                      <  2:1446‑M1/144‑27 (MQ=255)                             cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctttcaat                                                                                                             >  1:379650‑M1/1‑114 (MQ=255)                               aaaTCGACGGCACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                           <  2:655771‑M1/149‑38 (MQ=255)                                       tcagctgtaCTCGTCACTACACGCACATTAAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                   <  2:637393‑M1/145‑46 (MQ=255)                                        cAACTGTACTCGTCACTACACGCACATACTACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                                  <  1:572651‑M1/149‑47 (MQ=255)                                            tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                              >  2:575573‑M1/1‑99 (MQ=255)                                                 tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttaaaacgtttttccct                                                                                         >  2:120184‑M1/1‑94 (MQ=255)                                                 tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttattttgggtgtatttcaatctgtttaaaagtttttcgct                                                                                         >  2:27710‑M1/1‑94 (MQ=255)                                                       cTACACGCACATACAACGGAGGGGGGCCACGATTTTCAATAATGCCTGATGCAGATCACACAAAACACTCAATTACTCAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                   <  2:319698‑M1/149‑62 (MQ=255)                                                               acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                           >  1:265854‑M1/1‑80 (MQ=255)                                                                    cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                                                                      >  1:596840‑M1/1‑75 (MQ=255)                                                                      aCGGAGGGGGGCTGCGATTTTCAATCATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                    <  1:602782‑M1/149‑77 (MQ=255)                                                                             ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                             >  1:261344‑M1/1‑66 (MQ=255)                                                                                 cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                         >  1:148999‑M1/1‑62 (MQ=255)                                                                                  tGCGATTTTCAAGAAGGAGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctagtagaacatttttagctacgctagcaagcataactgaaaaacgaaaaatg                                                        >  2:136770‑M1/1‑61 (MQ=255)                                                                                       ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAGTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                   >  1:438165‑M1/1‑56 (MQ=255)                                                                                                      gATGCAGATCACACAAAACTCTCAATTACTTAACATTAATGtctaagtcatactttttttttgggtggatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                    <  2:369336‑M1/149‑109 (MQ=255)                                                                                                                acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                          >  1:77917‑M1/1‑31 (MQ=255)                                                                                                                          cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactccccgtaa                >  2:249747‑M1/1‑21 (MQ=255)                                                                                                                                  cTTAACATAAATGtataagtgatacttttgttttgggtgtatttcaatttgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatggaaatgggttgttgcgcacagtacacataccttaccgtaaaaattaat        >  2:289236‑M1/1‑13 (MQ=255)                                                                                                                                        aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacaggcaaaatagctcaccgtaaatcatcaactgcta  >  1:217017‑M1/1‑7 (MQ=255)                                                                                                                                              |                                                                                                                                               GGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_hsa_pgi/1978361‑1978644

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