Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	1,978,503	Δ776 bp	100%		insB1–insA	insB1, insA

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	32 [0]	[0] 32	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	0 (0.000)	28 (0.680)	20/260	NT	100%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913/1979135‑1979424                                                                                                                                                 |                                                                                                                                                  atgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                               <  2:454808‑M2/5‑1 (MQ=255) atgggcctttacacttatctaacatttttcctaaattgacgcaactttacacgtcactacactcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcatttacttaacataaaggTATaa                                                                                                                                               <  2:424160‑M2/5‑1 (MQ=255)      tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcaaacacaacactcaattacttaacataaatgTATAAGTCAt                                                                                                                                          >  1:299680‑M2/140‑149 (MQ=255)        tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                                                                        <  2:502721‑M2/12‑1 (MQ=255)                     aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                           >  1:708243‑M2/125‑149 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  2:763531‑M2/115‑149 (MQ=255)                               ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                 >  2:178606‑M2/115‑149 (MQ=255)                                taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                                >  1:232955‑M2/114‑149 (MQ=255)                                       acgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGGGTATTTCAATCTGTTaaa                                                                                                         >  1:118767‑M2/107‑149 (MQ=255)                                              tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                                                                  >  1:722156‑M2/100‑149 (MQ=255)                                                     tgcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactctattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                                                                           <  1:108331‑M2/57‑1 (MQ=255)                                                       cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                         <  1:446072‑M2/59‑1 (MQ=255)                                                         ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAATTTTTctctacgcta                                                                                       >  2:638507‑M2/89‑149 (MQ=255)                                                          tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                      >  2:97343‑M2/88‑149 (MQ=255)                                                          tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTag                                                                                      >  2:275364‑M2/88‑149 (MQ=255)                                                             acgcacatacaacggagtggtgctgcgattttcaatattgcttgatgcagatcacacaaaacactcaattacttaacataattgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAACTTTTTCGCTACGctagcta                                                                                   <  2:404504‑M2/65‑3 (MQ=255)                                                                            gggggggctgctattttcaataatgcgtgatgcagatcacaccacacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTGTTTCAATCTGTTAAAAAGGTTTTCGCTCCGCTAGCAAGCAAAAATGAAACAg                                                                    <  2:718014‑M2/80‑1 (MQ=255)                                                                                           ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCGAGCAAAAATGAAACAGGAATAATCGAAATgg                                                     >  2:624356‑M2/55‑149 (MQ=255)                                                                                           ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                                     <  1:5840‑M2/95‑1 (MQ=255)                                                                                                         atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                       <  2:708243‑M2/109‑1 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGGACAGTCaaaa                                 >  1:629482‑M2/35‑149 (MQ=255)                                                                                                               atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                                 >  1:466440‑M2/35‑149 (MQ=255)                                                                                                                      aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTCCGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACCCa                          >  1:367608‑M2/28‑149 (MQ=255)                                                                                                                      aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGGATTTCAATCTGTTAAAAAGTTTTTCTTTTCTTTTGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGGGCCCAGCCAAAATAACtcc                          >  2:419178‑M2/28‑148 (MQ=255)                                                                                                                       aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTCTTTTTTTTTTAGCAAAAAAGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                         >  2:458849‑M2/27‑149 (MQ=255)                                                                                                                       aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                         >  1:325498‑M2/27‑149 (MQ=255)                                                                                                                               aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAACCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                 >  2:76687‑M2/19‑149 (MQ=255)                                                                                                                               aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAACCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                 >  2:73523‑M2/19‑149 (MQ=255)                                                                                                                                  tacttaacataaatgTATAATTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  2:333636‑M2/134‑1 (MQ=255)                                                                                                                                  tacttaacataaacgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc              <  2:463584‑M2/134‑1 (MQ=255)                                                                                                                                            aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtata    >  2:336240‑M2/6‑149 (MQ=255)                                                                                                                                              atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaa  >  1:57332‑M2/4‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  ATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAA  >  NC_000913/1979135‑1979424

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