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	35 [0]	[1] 36	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)	33 (0.960)	29/260	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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTAT  >  NC_000913_3_hsa_tpiA/1978363‑1978639                                                                                                                                            |                                                                                                                                          gTCTTTACCATTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                  <  2:568581‑M1/149‑10 (MQ=255)  tCTTTAAACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                 <  1:486189‑M1/149‑11 (MQ=255)    tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACCCGCACATACAACGGAGGGGGGCTGCGAATTTCAATAATGCGTGATGCAGATCACACAAAAAAATCAATTACTTAAAATAAATTtataattcatac                                                                                                                               >  2:207032‑M1/1‑137 (MQ=255)    tttACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAAAACTCAATTACTTAACATAAATGtataagtcatac                                                                                                                               >  2:673166‑M1/1‑137 (MQ=255)         acTTATCTAAGATTTTTTTTTTTTTTTCGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCATATCACACAAAACTCTCAATTACTTAACATAAATGtataagtcatacttttt                                                                                                                          >  2:393341‑M1/1‑132 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACCCAAAACACTCAATTACTTAACATAAATGtataagtcataattttgttttggt                                                                                                                   >  2:334650‑M1/1‑125 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCACTTACTTAACATAAATGtataagtcatacttttgttttggg                                                                                                                   >  1:477439‑M1/1‑125 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACAGACAACGGAGGGGGGATGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtctaagtcatacttttgttttggg                                                                                                                   >  1:607526‑M1/1‑125 (MQ=255)                 aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATCCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttttttttgggt                                                                                                                  >  1:498590‑M1/1‑124 (MQ=255)                   gATTTTTCCTAAATCAACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                                                                <  1:111495‑M1/149‑28 (MQ=255)                     tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATATCTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                              >  1:472574‑M1/1‑120 (MQ=255)                        ttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACCTCCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcctacttttgttttgggtgtatttc                                                                                                           >  2:307724‑M1/1‑117 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                        >  1:163286‑M1/1‑114 (MQ=255)                            tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                       >  2:510839‑M1/1‑113 (MQ=255)                            tAAAGCGGCGCAAATGTGCTTGTTACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgtgtgtatttcaatc                                                                                                       <  2:307002‑M1/149‑37 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAACTTACTTAGCATATATGtataagtcatacttttgttttgggtgtctttcaatctgttaaaaactttttcgct                                                                                    >  2:76350‑M1/1‑94 (MQ=255)                                                 gTCACTACACGCACATACAACAGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttcaatctgttaaaaagtttttccctac                                                                                  >  2:307170‑M1/1‑92 (MQ=255)                                                         acGCTCATACGACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                          <  1:374971‑M1/149‑66 (MQ=255)                                                                tACAACGGAGGGGGGCTGCGATTTGCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                   <  2:368594‑M1/149‑73 (MQ=255)                                                                     cGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                                                              <  2:346203‑M1/149‑78 (MQ=255)                                                                           ggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaa                                                        >  2:658294‑M1/1‑66 (MQ=255)                                                                                    tgttttaaTGATGCGTGATGCAGATCACACAAAACACGCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                               <  1:334650‑M1/147‑93 (MQ=255)                                                                                        tCAAGAATGGGTGAGGCAGATGACCCAAAACTGTCTATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttacatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                           <  2:410725‑M1/149‑97 (MQ=255)                                                                                         cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga                                          <  1:547777‑M1/149‑98 (MQ=255)                                                                                                            tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatcggatgttgcccccattcaacat                       >  2:652236‑M1/1‑33 (MQ=255)                                                                                                             cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                      >  1:354865‑M1/1‑32 (MQ=255)                                                                                                             cacacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaata                      >  1:354775‑M1/1‑32 (MQ=255)                                                                                                              acataAAACTCTTAATTACTTAACACAAATGtataagtcatacttttgctttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                     <  1:170155‑M1/149‑119 (MQ=255)                                                                                                               cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgagaccagtcaacataac                    >  2:578161‑M1/1‑30 (MQ=255)                                                                                                                  aaaaCACTCTATTACTTAACTTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                 <  2:159504‑M1/149‑123 (MQ=255)                                                                                                                    aaCCCTCAATTACTTAACATAGATGtataaggcatagttttgttttgggggtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcacc               <  1:395754‑M1/149‑125 (MQ=255)                                                                                                                     acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg              >  1:228775‑M1/1‑24 (MQ=255)                                                                                                                            aatacTAAACATAAATGtataagtcatacttttgtgttgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataa       <  2:269518‑M1/147‑133 (MQ=255)                                                                                                                                cTTAACATAAATGtataagtcatacttttgtttggggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggattaatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat   <  1:449902‑M1/149‑137 (MQ=255)                                                                                                                                 tgaacaTAAATGgattagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcgccgtaaataatcgtc  <  2:478455‑M1/147‑138 (MQ=255)                                                                                                                                            |                                                                                                                                          GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTAT  >  NC_000913_3_hsa_tpiA/1978363‑1978639

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