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	36 [0]	[0] 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)	29 (0.890)	24/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

GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_hsa_tpiA/1978356‑1978648                                                                                                                                                   |                                                                                                                                                   aaaaTGGGTTCTAACACTTATCCAAGATCTTACCTAAATCAACGCAACTGTACTCGTCACTACACGCACATACAACGTAGCGGTCCTCCGATTATCAACAATGCGTAATGCAGATCACACAAAACACTCAATTACTTAACATAAATgta                                                                                                                                                  <  2:2950/148‑3 (MQ=255)  aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAGCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 >  1:28094‑M1/1‑146 (MQ=255)  aaaTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 >  1:28356‑M1/1‑146 (MQ=255)  aaaTGGGTCTTTACACTTATCTAAGATTTTCCCTAAATCGACGCAACTATACTCGTAACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                                                                 <  2:235952‑M1/149‑4 (MQ=255)    aTGGTTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCTATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               <  2:383770‑M1/149‑6 (MQ=255)              atccttaTCTAACATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                                                                     <  1:221461‑M1/146‑16 (MQ=255)               cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGCTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatactttt                                                                                                                                    >  1:248763‑M1/1‑133 (MQ=255)                   tacccaagatttCTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGATTTCCAACAATGCGTGATGCAGACCACACAAAACACTCAATTACTTAACATAAATGtataaggcacacccttgtct                                                                                                                                <  2:335899‑M1/144‑21 (MQ=255)                       tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGATTTTCAATAATGCGTGATGCAGATCACACACAACAGTCAGTTACTTAACATAAATGtataagtcatacttttgttttgga                                                                                                                            >  1:618041‑M1/1‑125 (MQ=255)                           aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAGCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                        >  1:239006‑M1/1‑121 (MQ=255)                            tttttCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtat                                                                                                                       <  2:452132‑M1/149‑30 (MQ=255)                                  cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCACTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                 >  2:565408‑M1/1‑114 (MQ=255)                                   ttaatcGACGCTACTGTATTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                <  1:509832‑M1/147‑37 (MQ=255)                                           cGCAACTGTACTCGTCACTACACGCACATACAACGGAGCGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaa                                                                                                        <  2:384398‑M1/149‑45 (MQ=255)                                            gCAACTGTACTCGTCACTACACGCACATACAACGGATGGGGGCTGCTATTTTCAATAATGCGTGATGCATATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaca                                                                                                       <  2:439556‑M1/149‑46 (MQ=255)                                               aCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtt                                                                                                    >  1:467359‑M1/1‑101 (MQ=255)                                                 tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  >  1:12187‑M1/1‑99 (MQ=255)                                                         tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgattttcttttattacagtctgttaaaaccttttttgctaca                                                                                          >  2:519016‑M1/1‑91 (MQ=255)                                                            cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgcta                                                                                       <  1:576387‑M1/149‑62 (MQ=255)                                                              acacGCAAATACAACGGAGGGAGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                     <  2:339361‑M1/149‑64 (MQ=255)                                                                             ggAGGGGCGCTGCAATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                                                                      <  2:643257‑M1/149‑79 (MQ=255)                                                                                    ggCTGCGATTTTCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaata                                                               <  1:352169‑M1/149‑86 (MQ=255)                                                                                      cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                             <  2:86489‑M1/149‑88 (MQ=255)                                                                                            ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                       >  1:384221‑M1/1‑56 (MQ=255)                                                                                                   taatGCGTGATGCAGATCACACAAAACCCTCAATTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatagcaatctgttaaaaagtgtttcgctacgctagcaagcaaaaatgaaacagcaataatcgaaatgggatgt                                                >  1:460309‑M1/1‑49 (MQ=255)                                                                                                         gTGATGCAGATCACTCAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttacaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                          <  2:137480‑M1/149‑107 (MQ=255)                                                                                                           gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                        <  1:591365‑M1/149‑109 (MQ=255)                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:24650‑M1/1‑31 (MQ=255)                                                                                                                     acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattacaatctgtgaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              >  1:364532‑M1/1‑31 (MQ=255)                                                                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                          >  1:352903‑M1/1‑27 (MQ=255)                                                                                                                         aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactca                          >  1:307761‑M1/1‑27 (MQ=255)                                                                                                                            acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccg                       >  1:528605‑M1/1‑24 (MQ=255)                                                                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttagaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaataggatgttgcgcacagtcaaaataactacccgtaaaca                 >  2:108087‑M1/1‑18 (MQ=255)                                                                                                                                  aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaacta                 >  1:108001‑M1/1‑18 (MQ=255)                                                                                                                                           aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctggtaaaaagtttttcgctacgctagcaagtaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgc        <  1:472324‑M1/149‑141 (MQ=255)                                                                                                                                                 aTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaccaggattaatggaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  >  1:483490‑M1/1‑3 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_hsa_tpiA/1978356‑1978648

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