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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503	1979278	776	34 [0]	[0] 34	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)	32 (0.890)	28/268	0.3	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

GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979141‑1979417                                                                                                                                           |                                                                                                                                           ttttagacttatctaagtggtttcctaaatcgacgcaactgtactcgtgactacacgcacatacaacggaggggggctgcgattttcaataatgagtgatgcagatcacacaaaatactcaattacttaacataaatgTATaa                                                                                                                                        <  2:262863‑M2/5‑1 (MQ=255) ctttacccttatcgaagattttttctaaatcgacgcaactgtactcttcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                        <  2:93678‑M2/5‑1 (MQ=255)  tttacacttatctaagatttttccgaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTAtaaa                                                                                                                                       >  1:102385‑M2/138‑142 (MQ=255)   ttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                                                                      <  1:579106‑M2/7‑1 (MQ=255)         ttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatacagatcacacaaaacactcaattacttaacataaatgTATAAGTCATacc                                                                                                                                >  2:284627‑M2/131‑142 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                           >  1:199133‑M2/126‑143 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                           >  2:334333‑M2/126‑143 (MQ=255)                  atctttcctaaagcgacgcaactgtacgcgtaactacacgcacatacgacggaggggggctgcgattttcattaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                       <  2:391811‑M2/22‑1 (MQ=255)                           aaatcgaggcagttgtactcgtccctacacgcgcatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                              <  2:627990‑M2/31‑1 (MQ=255)                              tcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                                                           >  1:446551‑M2/110‑143 (MQ=255)                                         gaactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacactctcaattacttaccataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTTAATCTGTTaaaaa                                                                                                <  2:138674‑M2/45‑1 (MQ=255)                                            ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTAATGTAGGGTATATATCaaaaagttaaaaagtt                                                                                             >  1:310184‑M2/96‑143 (MQ=255)                                             tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                            >  2:29080‑M2/95‑143 (MQ=255)                                               gtcactacacgcacatacaactgaggggggctgcgattttcaattatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                          <  2:298734‑M2/51‑1 (MQ=255)                                                 cactacacgcacataaaacggaggtgtgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaaggTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                        <  1:48726‑M2/53‑1 (MQ=255)                                                       acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                                                                  <  2:199133‑M2/59‑1 (MQ=255)                                                             agacaacggaggggggcggcgagtttcgataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                            <  1:55564‑M2/65‑1 (MQ=255)                                                               acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                              <  2:523242‑M2/31‑1 (MQ=255)                                                               acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                              >  1:523242‑M2/77‑107 (MQ=255)                                                                    ggaggggtgctgcgattttcaataatgcgttatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                     <  2:305463‑M2/72‑1 (MQ=255)                                                                    ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCACTACGCTAGCAGGCaaaaa                                                                     >  2:55025‑M2/72‑143 (MQ=255)                                                                      aggggggctgcgattttcaataatgcgtgttgctgattatacataatacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                   <  1:562836‑M2/74‑1 (MQ=255)                                                                         ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTCTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAACAAAAATGaaa                                                                >  2:145086‑M2/67‑143 (MQ=255)                                                                         ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                >  1:632451‑M2/67‑143 (MQ=255)                                                                             ctgcggttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGCAACAgg                                                            >  1:552081‑M2/63‑143 (MQ=255)                                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACAATAGCAAGCAAAAATGAAACAGGa                                                           >  1:656351‑M2/62‑143 (MQ=255)                                                                               gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                                                          >  1:715503‑M2/61‑143 (MQ=255)                                                                                                                aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCTCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                         >  1:559273‑M2/28‑143 (MQ=255)                                                                                                                  aacacttaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                       <  1:356651‑M2/118‑1 (MQ=255)                                                                                                                      ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAGGTTGCGCACAGGCCAAATAACTCa                   >  1:424391‑M2/22‑143 (MQ=255)                                                                                                                       tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTACGCACAATAAAAATAacacac                  >  2:411346‑M2/21‑143 (MQ=255)                                                                                                                              cttaacataattgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaat           <  2:102385‑M2/130‑1 (MQ=255)                                                                                                                                   cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatca      >  1:664172‑M2/9‑143 (MQ=255)                                                                                                                                       aacgTATAAGTCATACTTTTGTTTTGGGTGTATCTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg  <  2:145920‑M2/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                           GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979141‑1979417

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