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	44 [0]	[0] 46	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)	43 (1.070)	32/270	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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913_3_hsa_tpiA/1978359‑1978645                                                                                                                                                |                                                                                                                                                aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                            <  2:512851‑M1/149‑6 (MQ=255)      tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCGCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACGCTCAATTACTTAACATAGATGtataagtcat                                                                                                                                       >  2:191743‑M1/1‑139 (MQ=255)         ttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                    >  1:573474‑M1/1‑136 (MQ=255)             acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTCCTCGTCACTATACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATTtataagtcatacttttt                                                                                                                                >  2:463273‑M1/1‑132 (MQ=255)                 ctcgaagattTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCCCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                            <  1:465694‑M1/145‑22 (MQ=255)                 aTCAAAGAGTGTTCCGAAATGGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                            <  2:256401‑M1/149‑22 (MQ=255)                    tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatacgtcatactcttgttctggg                                                                                                                         >  2:5369‑M1/1‑125 (MQ=255)                    tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                         >  1:386967‑M1/1‑125 (MQ=255)                    tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAACTGcataagtcatacttttgtttttgc                                                                                                                         >  2:671503‑M1/1‑125 (MQ=255)                    tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCAAATACTACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCAGATCCCACAAAACACTCAATTACTTAACATAAATGtataagttatacttttgttttggg                                                                                                                         >  2:124650‑M1/1‑125 (MQ=255)                    tAAGATTTTTCCTAAATCGACGCAACTGACGATATTCCAAGACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATAACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                         >  2:162738‑M1/1‑125 (MQ=255)                     aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                                                                        >  1:74029‑M1/1‑124 (MQ=255)                            ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttc                                                                                                                 >  1:93494‑M1/1‑117 (MQ=255)                               cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggtggtatttcaat                                                                                                              >  2:265477‑M1/1‑114 (MQ=255)                               cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                              >  2:584280‑M1/1‑114 (MQ=255)                               cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATCAATTACTTAACATAAATTtataagtcatacttttgttttgtttgtatttcaat                                                                                                              >  2:11816‑M1/1‑114 (MQ=255)                                tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                             >  2:674423‑M1/1‑113 (MQ=255)                                   aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAAACTCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaatctgt                                                                                                          >  1:664496‑M1/1‑110 (MQ=255)                                    tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctttt                                                                                                         >  1:129164‑M1/1‑109 (MQ=255)                                           aaCTGTACTCGTCACTACACGCACATACAACGGGGGGGGGCTGCGATATTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgttaaaaagt                                                                                                  >  1:218404‑M1/1‑102 (MQ=255)                                                   tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttggttttatttcaatctgttaaaaagtttttcgct                                                                                          >  2:149595‑M1/1‑94 (MQ=255)                                                          tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtaaaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                   >  1:111487‑M1/1‑87 (MQ=255)                                                             acGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                                >  2:248084‑M1/1‑84 (MQ=255)                                                                    tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                         <  2:93494‑M1/149‑73 (MQ=255)                                                                     aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcggttaaaaagtttttcgctccgctagccatcaaaaat                                                                        >  2:25381‑M1/1‑76 (MQ=255)                                                                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                     <  2:380858‑M1/149‑77 (MQ=255)                                                                          ggAGGGCGGCTGCGATTTTCAATAATGAGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttagaaagttgttcgctacgctagcacgcacacatgaaca                                                                   >  2:281067‑M1/1‑71 (MQ=255)                                                                               ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                              <  2:74029‑M1/149‑84 (MQ=255)                                                                                   cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                          >  1:604448‑M1/1‑62 (MQ=255)                                                                                   cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                          >  1:31892‑M1/1‑62 (MQ=255)                                                                                    tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcataattttgttttgggtgtatttcaatctgttaaaacgtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                         >  1:588460‑M1/1‑61 (MQ=255)                                                                                          tatcaATAATGCGTGACGCAGATCCAACAAAAAACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                   <  1:625878‑M1/147‑95 (MQ=255)                                                                                                  atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgatagcaagcaaaaatgaaacaggaataatcgaaaggggatgttg                                           >  1:635002‑M1/1‑47 (MQ=255)                                                                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                       <  1:134311‑M1/149‑107 (MQ=255)                                                                                                      gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                       <  1:248084‑M1/149‑107 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaccgacaaaa                              >  1:620384‑M1/1‑34 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                              >  2:407046‑M1/1‑34 (MQ=255)                                                                                                                tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcccagtcaaaat                             >  1:418829‑M1/1‑33 (MQ=255)                                                                                                                 cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaacaatgaaacaggagtaatcgaaatgaacagttgcgcaatgtcataaaa                            >  2:21944‑M1/1‑32 (MQ=255)                                                                                                                    acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                         >  1:201930‑M1/1‑29 (MQ=255)                                                                                                                    acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                         >  1:695449‑M1/1‑29 (MQ=255)                                                                                                                      aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                       <  1:481607‑M1/149‑123 (MQ=255)                                                                                                                                        aCATAAATGtataagtcacacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc     <  1:276782‑M1/149‑141 (MQ=255)                                                                                                                                           tAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctat  <  1:72593‑M1/149‑144 (MQ=255)                                                                                                                                                |                                                                                                                                                ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCT  >  NC_000913_3_hsa_tpiA/1978359‑1978645

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