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–1979269	1979278	10–776	35 [0]	[0] 35	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)	35 (0.950)	31/266	0.3	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

TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_hsa_tpiA/1978354‑1978643                                                                                                                                                     |                                                                                                                                              tAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTCAATCGACGCAACTGTACTCGTCCCTACACGCCCATACAACGGAGGGGGGCGGCGATTTTCAATAATGCGTGATGCACATCACACAAAACACCAACTTACTTAACTTACATGt                                                                                                                                              >  2:553985‑M1/1‑149 (MQ=255)     aaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTATCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                          <  2:349295‑M1/150‑6 (MQ=255)     aaTGGGCCTTTACACTTATCTAAGATTTTTCCTAAAGCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCTATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                          <  2:251604‑M1/150‑6 (MQ=255)          gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACCCAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                     >  2:639316‑M1/1‑140 (MQ=255)                  acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCCCACACAACACTCAATTTCTTAACATAAATGtataagtcatacttttgt                                                                                                                             >  1:190042‑M1/1‑132 (MQ=255)                   cTTATCTAAGATTTTTCCCAAACCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtt                                                                                                                            <  1:319659‑M1/150‑20 (MQ=255)                      aTCTAAGATCTTTACTAAATCGACGCAACTGTACCCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACACAAATGtataagtcatacttttgttttg                                                                                                                         <  1:206106‑M1/150‑23 (MQ=255)                       tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                        <  1:326944‑M1/150‑24 (MQ=255)                         tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                      <  2:207947‑M1/150‑26 (MQ=255)                          aaGATTTTTCCTAAATGGACGCAACTGTACTCGTCGCTACACGCATATAAAACGGAGGGGGGCTGCGATTTTCAATGATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAGATAACTGtataagtcattattttgctttgggtg                                                                                                                     >  2:96231‑M1/1‑124 (MQ=255)                           aGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                                                                    >  1:389424‑M1/1‑123 (MQ=255)                               ttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                <  2:385127‑M1/150‑32 (MQ=255)                                     tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                         <  1:447035‑M1/135‑23 (MQ=255)                                     tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                         >  2:447035‑M1/1‑113 (MQ=255)                                           gACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaa                                                                                                    >  1:445202‑M1/1‑107 (MQ=255)                                            aCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                   >  2:296233‑M1/1‑106 (MQ=255)                                             cGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                  <  1:36012‑M1/150‑46 (MQ=255)                                              gCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                                                                 >  1:503334‑M1/1‑104 (MQ=255)                                                aaCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                               >  2:164283‑M1/1‑102 (MQ=255)                                                                    gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                                                           >  1:107635‑M1/1‑82 (MQ=255)                                                                        aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                       >  1:62810‑M1/1‑78 (MQ=255)                                                                             aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                                  >  2:275921‑M1/1‑73 (MQ=255)                                                                                gAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                               <  2:503334‑M1/150‑81 (MQ=255)                                                                                   gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                            <  1:275921‑M1/150‑84 (MQ=255)                                                                                      ggttgcgATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                         <  2:189622‑M1/147‑87 (MQ=255)                                                                                       gCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcctacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                        <  1:633475‑M1/150‑88 (MQ=255)                                                                                                 tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                              <  1:164283‑M1/150‑98 (MQ=255)                                                                                                     taatGCGTGAAGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgtt                                          <  1:553037‑M1/150‑102 (MQ=255)                                                                                                       atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                        <  1:97449‑M1/150‑104 (MQ=255)                                                                                                       atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                        <  2:445202‑M1/150‑104 (MQ=255)                                                                                                              aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                 <  1:624430‑M1/150‑111 (MQ=255)                                                                                                                          caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                     >  2:276294‑M1/1‑28 (MQ=255)                                                                                                                          caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                     >  1:57545‑M1/1‑28 (MQ=255)                                                                                                                                    aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa           >  1:355401‑M1/1‑18 (MQ=255)                                                                                                                                             aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct  <  1:296233‑M1/150‑142 (MQ=255)                                                                                                                                                     |                                                                                                                                              TAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCG  >  NC_000913_3_hsa_tpiA/1978354‑1978643

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