Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913_3_hsa_tpiA	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_tpiA	1978503	1979278	776	46 [1]	[1] 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)	39 (0.930)	34/258	NT	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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_hsa_tpiA/1978363‑1978646                                                                                                                                            |                                                                                                                                                 gTCTCTACACTTATCTAACATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                         <  1:358386‑M1/149‑10 (MQ=255)     ttACACTTATCTAAGATTTTTCCTAAATCTACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATCCGTGATTCAAATCACACAAAACACCTAATTACTTAACATAAATGtataagtcatact                                                                                                                                     >  2:215286‑M1/1‑136 (MQ=255)             aTCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                                                             <  1:427456‑M1/149‑22 (MQ=255)               cTAAGATTTTTCCAAAATCGACGCATCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                           <  2:651805‑M1/149‑24 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttgtggg                                                                                                                          >  1:490383‑M1/1‑125 (MQ=255)                    aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                      >  2:815342‑M1/1‑121 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                               >  1:345411‑M1/1‑114 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAAAGGAGGGGGGCTGAGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATAAATTCCTTAACATAAATGtataagtcataattttgttttgggtgtctttaaat                                                                                                               >  2:184955‑M1/1‑114 (MQ=255)                              aaTTGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtatacgtcatacttttgttttgggtgtatttcaatctg                                                                                                            <  2:366800‑M1/149‑39 (MQ=255)                                tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGAGGGCTGCGATTTTCAATAATGCGTGATGCGGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttggtttgggtgtatttcaatctgtt                                                                                                          >  1:201586‑M1/1‑109 (MQ=255)                                     gCAACTGTTCTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaa                                                                                                     <  1:746166‑M1/149‑46 (MQ=255)                                        aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtc                                                                                                  >  2:276963‑M1/1‑101 (MQ=255)                                          tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                >  1:788733‑M1/1‑99 (MQ=255)                                           gTACTCGTTAATACACGGACATACAACGGTGGGGGGCTGCGATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               <  2:235925‑M1/149‑52 (MQ=255)                                           gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               <  1:495797‑M1/149‑52 (MQ=255)                                           gTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttt                                                                                               >  1:679437‑M1/1‑98 (MQ=255)                                            tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                              <  1:309466‑M1/149‑53 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCACTAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                           >  2:33777‑M1/1‑94 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACCAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaatgttttcacg                                                                                           >  1:59598‑M1/1‑94 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                           >  2:534806‑M1/1‑94 (MQ=255)                                                cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                          <  1:201915‑M1/149‑57 (MQ=255)                                                      tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctactctac                                                                                    >  2:380264‑M1/1‑87 (MQ=255)                                                                tACAACGGAGGGGGGCTGCGATGTTCAATAATGCGTCATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttccgctacgctagcaagcaaaaa                                                                          <  2:262345‑M1/149‑73 (MQ=255)                                                                 aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTGCTTAACATAAATGtatacgtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaat                                                                         >  1:327737‑M1/1‑76 (MQ=255)                                                                   aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatcttttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       <  2:561070‑M1/149‑76 (MQ=255)                                                                   aaCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatga                                                                       >  1:657776‑M1/1‑74 (MQ=255)                                                                    acggctgggggctgcgaTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                      >  2:312045‑M1/7‑73 (MQ=255)                                                                        aGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacag                                                                  >  1:321791‑M1/1‑69 (MQ=255)                                                                          gggggCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagaaaaaatgaaacagga                                                                >  2:814861‑M1/1‑67 (MQ=255)                                                                                             aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactcttgttttgggtgtatttcaatctgttaaaaagttcttcgctacgctagcaagcaaaaatgaaacaggaataatctaaatggcataca                                             >  2:838400‑M1/1‑48 (MQ=255)                                                                                               tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtgtaagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                           >  2:788785‑M1/1‑46 (MQ=255)                                                                                               tGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                           <  2:528605‑M1/149‑104 (MQ=255)                                                                                                gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                          >  2:605975‑M1/1‑45 (MQ=255)                                                                                                        cAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtca                                  >  1:471751‑M1/1‑37 (MQ=255)                                                                                                          gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttctctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacggtcaaa                                >  2:190527‑M1/1‑35 (MQ=255)                                                                                                           aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                               >  1:634904‑M1/1‑34 (MQ=255)                                                                                                                   aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcac                       <  2:56262‑M1/149‑124 (MQ=255)                                                                                                                    aaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc                      >  2:658118‑M1/1‑25 (MQ=255)                                                                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaata               >  1:631591‑M1/1‑18 (MQ=255)                                                                                                                                 ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatc         <  1:215286‑M1/149‑138 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaataatcatctg       >  1:262117‑M1/1‑10 (MQ=255)                                                                                                                                     cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaaacaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     >  2:215725‑M1/1‑8 (MQ=255)                                                                                                                                     cATAAACGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct     <  1:17441‑M1/149‑142 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta    >  1:519171‑M1/1‑7 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgctaaaaagtttttcgctacgctagcaagcacaaatgaaacagcaataatcgaaatgggaagttgcacacagtcaaaatacctcaacgtaaataatcatcctcta    >  2:682799‑M1/1‑7 (MQ=255)                                                                                                                                        attggggTAATGACTCCAACTTATTGACAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTTTGCCGCTc  >  1:177056/5‑149 (MQ=2)                                                                                                                                            |                                                                                                                                                 GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTC  >  NC_000913_3_hsa_tpiA/1978363‑1978646

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