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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_hsa_tpiA	1978503	1979278	776	40 [0]	[0] 40	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)	36 (0.600)	30/264	1.5	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

AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_hsa_tpiA/1978357‑1978650                                                                                                                                                  |                                                                                                                                                     aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACCTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                   >  2:555308‑M1/1‑146 (MQ=255) aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCATCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTTATGCAGATCACACAAAACACTCAATTACTTAAAATAAATgaat                                                                                                                                                   >  1:949102/1‑146 (MQ=255)         cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACGCGCACATACAACGGTGGGGGGCTGCGATTTTCAACAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatc                                                                                                                                           <  2:860290‑M1/149‑12 (MQ=255)           ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAGATAAATGtataagtcatact                                                                                                                                         >  1:807255‑M1/1‑136 (MQ=255)                     cTAAAAATTTTCCTAAATCGACGCAACTATACCCCTCACCACACGCACAAACAACGAAGGGAGGCTGCGATTTTCAATAATGCGTGATACAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgt                                                                                                                               <  2:263341‑M1/149‑24 (MQ=255)                      tAAGATTTTTCCCAATTCGACGATACCCGACCCGTCACTACACACACATACAACGGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACCCAATTACTTAACATAAATTtataagtcatacttttgttttggg                                                                                                                              <  1:1055411‑M1/149‑25 (MQ=255)                         gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                                                                           >  1:363326‑M1/1‑122 (MQ=255)                              ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                                                                      <  1:987610‑M1/149‑33 (MQ=255)                                ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                                    >  1:265229‑M1/1‑115 (MQ=255)                                                    cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                                                                                <  2:157428‑M1/149‑55 (MQ=255)                                                     tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATATTTACTTTTACGAGATGCAAATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagct                                                                                               >  1:268450‑M1/1‑94 (MQ=255)                                                               acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttagctacgctagcaa                                                                                     >  1:420074‑M1/1‑84 (MQ=255)                                                                 gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                   >  2:679491‑M1/1‑82 (MQ=255)                                                                   acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGGAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcggtaagctagcaagcaa                                                                                 >  1:14002‑M1/1‑80 (MQ=255)                                                                   acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagaaagaaa                                                                                 >  2:506665‑M1/1‑80 (MQ=255)                                                                     aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                      >  2:4395‑M1/1‑78 (MQ=255)                                                                     aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                      <  1:4395‑M1/94‑17 (MQ=255)                                                                      tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                              >  1:511683‑M1/1‑77 (MQ=255)                                                                             gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGAACAAAAAAAACAAACAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacaaaaaaaaaaaaaaatgaaaca                                                                       <  2:686832‑M1/149‑80 (MQ=255)                                                                                     cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                               >  1:463250‑M1/1‑62 (MQ=255)                                                                                                  taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatacaagaaggataatcgaattgggatgt                                                  >  1:458938‑M1/1‑49 (MQ=255)                                                                                                  aaaagcgtgAAGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaaaaaaaaaaaaaaaaatgaaacaggaataatcgaaatgggatgt                                                  <  2:463250‑M1/145‑101 (MQ=255)                                                                                                          gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggagtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgtaatgggatgttgtgcacg                                          >  2:345261‑M1/1‑41 (MQ=255)                                                                                                           aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                         >  2:173557‑M1/1‑40 (MQ=255)                                                                                                                 attacacaaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcactctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgatatgggatgttgcgcacagtcaaaa                                   <  2:981644‑M1/146‑116 (MQ=255)                                                                                                                  tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaat                                  >  2:967200‑M1/1‑33 (MQ=255)                                                                                                                      acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggaatgttgcgcacagtcaaaataact                              >  1:589495‑M1/1‑29 (MQ=255)                                                                                                                        aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcgagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                            >  2:1008434‑M1/1‑27 (MQ=255)                                                                                                                           acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                         >  1:509213‑M1/1‑24 (MQ=255)                                                                                                                             acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaactaaatcaccgta                       >  2:853686‑M1/1‑22 (MQ=255)                                                                                                                             acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta                       >  1:440453‑M1/1‑22 (MQ=255)                                                                                                                              cTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaa                      >  2:991659‑M1/1‑21 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                   >  2:883730‑M1/1‑18 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                   >  1:414507‑M1/1‑18 (MQ=255)                                                                                                                                 aaTTACTTAACATAAATGtataactcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgttacgctagcaagacaaagtgaaacaggaataatcgaaataggatgttgcgcacagtcaaacaaacacacagtaaata                   >  2:1046368‑M1/1‑18 (MQ=255)                                                                                                                                       ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc             >  1:1088898‑M1/1‑12 (MQ=255)                                                                                                                                          aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc          >  1:938971‑M1/1‑9 (MQ=255)                                                                                                                                            aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaactaatcatctgcta        >  1:230188‑M1/1‑7 (MQ=255)                                                                                                                                                aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa    >  2:779604‑M1/1‑3 (MQ=255)                                                                                                                                                  gtataattcatacttttgttttgtgtgtatttcaatctgttaaaaagtttttcgctacgttagcaagtaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaata  <  1:70131‑M1/149‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     AAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATT  >  NC_000913_3_hsa_tpiA/1978357‑1978650

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