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	40 [0]	[0] 40	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)	36 (0.600)	30/264	1.5	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                                                                                                                                                   |                                                                                                                                                     aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacctacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                                    >  2:555308‑M2/147‑149 (MQ=255)         ctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacgcgcacatacaacggtggggggctgcgattttcaacaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCatc                                                                                                                                            <  2:860290‑M2/11‑2 (MQ=255)           ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaagataaatgTATAAGTCATACt                                                                                                                                          >  1:807255‑M2/137‑149 (MQ=255)                     ctaaaaattttcctaaatcgacgcaactatacccctcaccacacgcacaaacaacgaagggaggctgcgattttcaataatgcgtgatacagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttgttttgt                                                                                                                                <  2:263341‑M2/23‑2 (MQ=255)                      taagatttttcccaattcgacgatacccgacccgtcactacacacacatacaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacacccaattacttaacataaattTATAAGTCATACTTTTGTTTTggg                                                                                                                               <  1:1055411‑M2/24‑1 (MQ=255)                         gatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                                                                                                            >  1:363326‑M2/123‑149 (MQ=255)                              ttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                                                                       <  1:987610‑M2/32‑1 (MQ=255)                                cctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                                                                     >  1:265229‑M2/116‑149 (MQ=255)                                                    ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                                                                 <  2:157428‑M2/54‑1 (MQ=255)                                                     tcgtcactacacgcacatacaacggaggggggctgcgatatttacttttacgagatgcaaatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTagct                                                                                                >  1:268450‑M2/95‑149 (MQ=255)                                                               acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTagcaa                                                                                      >  1:420074‑M2/85‑149 (MQ=255)                                                                 gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                                                                    >  2:679491‑M2/83‑149 (MQ=255)                                                                   acatacaacggaggggggctgcgattttcaataatgcgtgatggagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGGTAAGCTAgcaagcaa                                                                                  >  1:14002‑M2/81‑149 (MQ=255)                                                                   acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgaaagaaa                                                                                  >  2:506665‑M2/81‑149 (MQ=255)                                                                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                       >  2:4395‑M2/79‑94 (MQ=255)                                                                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                                       <  1:4395‑M2/16‑1 (MQ=255)                                                                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                               >  1:511683‑M2/78‑149 (MQ=255)                                                                             gaggggggctgcgattttcaataatgcgtgatgcagaacaaaaaaaacaaacaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACAAAAAAAAAAAAAAATGAAACa                                                                        <  2:686832‑M2/79‑1 (MQ=255)                                                                                     ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                                >  1:463250‑M2/63‑149 (MQ=255)                                                                                                  taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATACAAGAAGGATAATCGAATTGGGATGt                                                   >  1:458938‑M2/50‑149 (MQ=255)                                                                                                  aaaagcgtgaagcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAAAAAAAAAAAAAAAATGAAACAGGAATAATCGAAATGGGATGt                                                   <  2:463250‑M2/100‑1 (MQ=255)                                                                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGAGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGTAATGGGATGTTGTgcacg                                           >  2:345261‑M2/42‑148 (MQ=255)                                                                                                           atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                          >  2:173557‑M2/41‑149 (MQ=255)                                                                                                                 attacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCACTCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAACAGGAATAATCGATATGGGATGTTGCGCACAGTCaaaa                                    <  2:981644‑M2/115‑1 (MQ=255)                                                                                                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                                   >  2:967200‑M2/34‑149 (MQ=255)                                                                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGAATGTTGCGCACAGTCAAAATAACt                               >  1:589495‑M2/30‑149 (MQ=255)                                                                                                                        aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCa                             >  2:1008434‑M2/28‑149 (MQ=255)                                                                                                                           acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                          >  1:509213‑M2/25‑149 (MQ=255)                                                                                                                             actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAACTAAATCACCGTa                        >  2:853686‑M2/23‑149 (MQ=255)                                                                                                                             actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTa                        >  1:440453‑M2/23‑149 (MQ=255)                                                                                                                              ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa                       >  2:991659‑M2/22‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                    >  2:883730‑M2/19‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                    >  1:414507‑M2/19‑149 (MQ=255)                                                                                                                                 aattacttaacataaatgTATAACTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGTTACGCTAGCAAGACAAAGTGAAACAGGAATAATCGAAATAGGATGTTGCGCACAGTCAAACAAACACACAGTAaata                    >  2:1046368‑M2/19‑149 (MQ=255)                                                                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcatc              >  1:1088898‑M2/13‑149 (MQ=255)                                                                                                                                          acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc           >  1:938971‑M2/10‑149 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAACTAATCATCTGCta         >  1:230188‑M2/8‑149 (MQ=255)                                                                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa     >  2:779604‑M2/4‑149 (MQ=255)                                                                                                                                                  gTATAATTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGTTAGCAAGTAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaata   <  1:70131‑M2/148‑1 (MQ=255)                                                                                                                                                   tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa  <  1:888403/149‑1 (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: ‑