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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_hsa_tpiA/1979133‑1979427                                                                                                                                                   |                                                                                                                                                     aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacagcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    >  1:28094‑M2/147‑149 (MQ=255) aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    >  1:28356‑M2/147‑149 (MQ=255) aaatgggtctttacacttatctaagattttccctaaatcgacgcaactatactcgtaactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    <  2:235952‑M2/3‑1 (MQ=255)   atggttctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgctattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                                  <  2:383770‑M2/5‑1 (MQ=255)             atccttatctaacatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                                                                        <  1:221461‑M2/15‑1 (MQ=255)              cacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgctgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACtttt                                                                                                                                       >  1:248763‑M2/134‑149 (MQ=255)                  tacccaagatttctcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggagctgcgatttccaacaatgcgtgatgcagaccacacaaaacactcaattacttaacataaatgTATAAGGCACacccttgtct                                                                                                                                   <  2:335899‑M2/20‑3 (MQ=255)                      taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggagctgcgattttcaataatgcgtgatgcagatcacacacaacagtcagttacttaacataaatgTATAAGTCATACTTTTGTTTtgga                                                                                                                               >  1:618041‑M2/126‑148 (MQ=255)                          atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaagcactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                           >  1:239006‑M2/122‑149 (MQ=255)                           tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                          <  2:452132‑M2/29‑1 (MQ=255)                                 ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcacttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                    >  2:565408‑M2/115‑149 (MQ=255)                                  ttaatcgacgctactgtattcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                   <  1:509832‑M2/36‑1 (MQ=255)                                          cgcaactgtactcgtcactacacgcacatacaacggagcggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaa                                                                                                           <  2:384398‑M2/44‑1 (MQ=255)                                           gcaactgtactcgtcactacacgcacatacaacggatgggggctgctattttcaataatgcgtgatgcatatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGttaaaca                                                                                                          <  2:439556‑M2/45‑3 (MQ=255)                                              actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                                                                       >  1:467359‑M2/102‑149 (MQ=255)                                                tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                     >  1:12187‑M2/100‑149 (MQ=255)                                                        tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGATTTTCTTTTATTACAGTCTGTTAAAACCTTTTttgctaca                                                                                             >  2:519016‑M2/92‑148 (MQ=255)                                                           ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                          <  1:576387‑M2/61‑1 (MQ=255)                                                             acacgcaaatacaacggagggaggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                        <  2:339361‑M2/63‑1 (MQ=255)                                                                            ggaggggcgctgcaattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                                                                         <  2:643257‑M2/78‑1 (MQ=255)                                                                                   ggctgcgattttcaataatgcgggatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                                  <  1:352169‑M2/85‑1 (MQ=255)                                                                                     ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                                <  2:86489‑M2/87‑1 (MQ=255)                                                                                           ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt                                                          >  1:384221‑M2/57‑149 (MQ=255)                                                                                                  taatgcgtgatgcagatcacacaaaaccctcaattacttaacataaatgTATAAGTCATACTTTTGTGTTGGGTGTATAGCAATCTGTTAAAAAGTGTTTCGCTACGCTAGCAAGCAAAAATGAAACAGCAATAATCGAAATGGGATGt                                                   >  1:460309‑M2/50‑149 (MQ=255)                                                                                                        gtgatgcagatcactcaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTACAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                             <  2:137480‑M2/106‑1 (MQ=255)                                                                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                           <  1:591365‑M2/108‑1 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                 >  1:24650‑M2/32‑149 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTACAATCTGTGAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                                 >  1:364532‑M2/32‑149 (MQ=255)                                                                                                                        aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                             >  1:352903‑M2/28‑149 (MQ=255)                                                                                                                        aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                             >  1:307761‑M2/28‑149 (MQ=255)                                                                                                                           acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                          >  1:528605‑M2/25‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAGAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATAGGATGTTGCGCACAGTCAAAATAACTACCCgtaaaca                    >  2:108087‑M2/19‑147 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAacta                    >  1:108001‑M2/19‑149 (MQ=255)                                                                                                                                          acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGGTAAAAAGTTTTTCGCTACGCTAGCAAGTAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc           <  1:472324‑M2/140‑1 (MQ=255)                                                                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAACCAGGATTAATGGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     >  1:483490‑M2/4‑149 (MQ=255)                                                                                                                                                   tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGGTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  >  1:432923/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAA  >  NC_000913_3_hsa_tpiA/1979133‑1979427

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