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	42 [0]	[0] 44	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)	42 (0.960)	34/264	0.4	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

TTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGT  >  NC_000913_3_hsa_tpiA/1978367‑1978636                                                                                                                                        |                                                                                                                                       tgacacttATCTTAGATGTTTCCGAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                           >  2:117702‑M1/3‑136 (MQ=255)   acacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGTGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                         <  1:583616‑M1/149‑16 (MQ=255)           cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                 <  1:535730‑M1/149‑24 (MQ=255)            tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                >  1:133871‑M1/1‑125 (MQ=255)             aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                               <  2:448955‑M1/149‑26 (MQ=255)               gATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatacttttgttttgggtgt                                                                                                             >  2:362921‑M1/1‑122 (MQ=255)                 tttttCCTAAATCGACGCAACTGTACTCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAGAAGACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                           <  1:644559‑M1/149‑30 (MQ=255)                         aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                   >  2:233772‑M1/1‑112 (MQ=255)                            tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                >  1:157712‑M1/1‑109 (MQ=255)                                 gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTCACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                >  1:362482‑M1/1‑104 (MQ=255)                                 gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                           >  2:81534‑M1/1‑104 (MQ=255)                                 gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                <  2:362482‑M1/144‑41 (MQ=255)                                       gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                     >  1:554813‑M1/1‑98 (MQ=255)                                           tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                 >  1:74152‑M1/1‑94 (MQ=255)                                             gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctttcaatctgttaaaaagtttttcgctac                                                                               >  1:151351‑M1/1‑92 (MQ=255)                                             gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctag                                                                               >  1:151244‑M1/1‑92 (MQ=255)                                             gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                               >  2:594143‑M1/1‑92 (MQ=255)                                                  tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGGTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                          >  1:541333‑M1/1‑87 (MQ=255)                                                     acGGACATACAACGGAGGGTGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                       <  2:408490‑M1/149‑66 (MQ=255)                                                         acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                 >  2:385095‑M1/1‑80 (MQ=255)                                                         acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                 <  1:385095‑M1/119‑40 (MQ=255)                                                          gttccaccggaggggGGCTCCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCACTTACTTAACGTACATCgccgtgacatacttttgttttgggtggatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                                                                  <  1:362921‑M1/142‑71 (MQ=255)                                                             aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                               <  1:213829‑M1/149‑74 (MQ=255)                                                               aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATGGCACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttccatctgttaaaaagtttttcgctacgctagcaagaaaaaatga                                                             >  2:353810‑M1/1‑74 (MQ=255)                                                                      gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagga                                                      <  1:488269‑M1/149‑83 (MQ=255)                                                                          gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                  >  2:546854‑M1/1‑63 (MQ=255)                                                                            tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                >  1:175628‑M1/1‑61 (MQ=255)                                                                            tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                >  2:149227‑M1/1‑61 (MQ=255)                                                                                   ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                         <  2:74152‑M1/149‑96 (MQ=255)                                                                                                gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                            <  1:117702‑M1/149‑109 (MQ=255)                                                                                                          acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaagctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                  <  1:233772‑M1/149‑119 (MQ=255)                                                                                                          acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                      <  1:776928‑M1/65‑35 (MQ=255)                                                                                                          acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaa                                                                                                      >  2:776928‑M1/1‑31 (MQ=255)                                                                                                            acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctttttgtttttaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcccagtcaaaataact                >  2:73593‑M1/1‑29 (MQ=255)                                                                                                              aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                         <  1:729172‑M1/122‑96 (MQ=255)                                                                                                              aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg                                         >  2:729172‑M1/1‑27 (MQ=255)                                                                                                                   acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgta         >  2:602059‑M1/1‑22 (MQ=255)                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata     >  2:578444‑M1/1‑18 (MQ=255)                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                   <  2:360882‑M1/119‑102 (MQ=255)                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                   >  1:360882‑M1/1‑18 (MQ=255)                                                                                                                       aaTTACTTAACATAAATGtataagtcatacttttgtttagggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata     >  2:716940‑M1/1‑18 (MQ=255)                                                                                                                          tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc  <  1:149227‑M1/149‑135 (MQ=255)                                                                                                                                        |                                                                                                                                       TTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGT  >  NC_000913_3_hsa_tpiA/1978367‑1978636

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑