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	51 [0]	[0] 50	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)	47 (0.800)	46/264	0.3	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

TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_tpiA/1979133‑1979420                                                                                                                                                   |                                                                                                                                              aaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactaaacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                                                             <  1:533263‑M2/3‑1 (MQ=255)   attgttcttgacactagtcgaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                           <  2:825348‑M2/5‑1 (MQ=255)    tgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                          >  2:407545‑M2/144‑149 (MQ=255)          tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                    <  1:688652‑M2/12‑1 (MQ=255)                 ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgatttttaataatgcgtgatgcagatcacacaaaacactcaattacttaacatagatgTATAAGTCATACttttgtt                                                                                                                             <  1:81835‑M2/19‑1 (MQ=255)                 ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                             >  1:127206‑M2/131‑149 (MQ=255)                   atctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacgacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                           <  1:606087‑M2/21‑1 (MQ=255)                       aacatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                       <  1:507805‑M2/25‑1 (MQ=255)                        agatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                                                                      >  2:385698‑M2/124‑149 (MQ=255)                          atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataactgTATAAGTCATAATTTTgttttggtttta                                                                                                                    >  2:205073‑M2/122‑144 (MQ=255)                          atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                                                                    >  2:774495‑M2/122‑149 (MQ=255)                                 ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                             >  2:321211‑M2/115‑149 (MQ=255)                                   aaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                           >  2:749478‑M2/113‑149 (MQ=255)                                           gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                   <  1:777288‑M2/45‑1 (MQ=255)                                             aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                                                                 >  2:797986‑M2/103‑149 (MQ=255)                                                     tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatcTATAAGTCATACTTTTCTTTTGGGTGTATTTCAATCTGTTAAAAAATTTTTtgct                                                                                         >  2:273420‑M2/95‑149 (MQ=255)                                                       gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                       >  2:889840‑M2/93‑149 (MQ=255)                                                        tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                      >  1:390692‑M2/92‑149 (MQ=255)                                                         caatacaagcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                     <  2:311988‑M2/59‑1 (MQ=255)                                                              cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                                                                >  1:850860‑M2/86‑149 (MQ=255)                                                               acgaccatacaacggatgagggctccgattttcaataatgcgttatgcagatcacacaaaacactcaattacttaacctaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                               <  1:734220‑M2/65‑1 (MQ=255)                                                                cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                              <  1:774495‑M2/66‑1 (MQ=255)                                                                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTACAATCGTTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                         >  1:853983‑M2/79‑149 (MQ=255)                                                                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                                                                     >  1:359132‑M2/75‑149 (MQ=255)                                                                           cggaggggggctgcgattttcaataaggcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATCTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                   <  2:753174‑M2/77‑1 (MQ=255)                                                                                  gggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                            <  1:22987‑M2/84‑1 (MQ=255)                                                                                   ggctgcgattttcaataaggctggattcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                           <  1:419417‑M2/85‑1 (MQ=255)                                                                                      tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                        >  2:686883‑M2/62‑149 (MQ=255)                                                                                         gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                     >  1:890251‑M2/59‑149 (MQ=255)                                                                                               caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                               <  1:889840‑M2/97‑1 (MQ=255)                                                                                                  taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                               <  2:147985‑M2/49‑1 (MQ=255)                                                                                                  taatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                               >  1:147985‑M2/50‑98 (MQ=255)                                                                                                   aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                           <  2:198792‑M2/101‑1 (MQ=255)                                                                                                    atgcgtgatgcagatcacacaaaacactcaattacttaacatgaatgTATAAGTCATACTTTTGTTTTGGGTGTACTTCAATCTGTTAAAAAGTTTTGCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                          >  2:1039383‑M2/48‑149 (MQ=255)                                                                                                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                      <  1:797986‑M2/106‑1 (MQ=255)                                                                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTATTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCAAAAAGGTATGTTGCGCACAGt                                  >  2:467057‑M2/40‑149 (MQ=255)                                                                                                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGAGATGTTTGGCACAGTc                                 >  2:763025‑M2/39‑149 (MQ=255)                                                                                                               agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATAGAAATGGAAATTTTCGTACAGTCaa                               >  2:722659‑M2/37‑149 (MQ=255)                                                                                                                gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATAGAAATGGGATGTTGCGCACAGTCaaa                              >  2:726643‑M2/36‑149 (MQ=255)                                                                                                                 ataacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                             <  1:749478‑M2/115‑1 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                             <  1:878452‑M2/35‑1 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                             >  2:878452‑M2/32‑66 (MQ=255)                                                                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                          <  1:321211‑M2/118‑1 (MQ=255)                                                                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACt                        >  2:576596‑M2/30‑149 (MQ=255)                                                                                                                          aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAGGAATAATCGAAATTGGATGTTGAGCACAGTCAAAACAACACAcc                    >  2:374531‑M2/26‑149 (MQ=255)                                                                                                                                      cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat        >  1:865096‑M2/14‑149 (MQ=255)                                                                                                                                         aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg     >  1:1009699‑M2/11‑149 (MQ=255)                                                                                                                                           cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAACAATCATCTGCt   <  1:691968‑M2/141‑1 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCCAGCAAAAATGAAACAGGAATAATAGAAATGGGATGTTGCGCACAGTCAGAATAACTCACCGTGAATAATCATCTGCta  >  1:548774‑M2/8‑149 (MQ=255)                                                                                                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCta  >  2:920862‑M2/8‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                              TGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTA  >  NC_000913_3_hsa_tpiA/1979133‑1979420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑