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	30 [0]	[0] 30	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)	29 (1.030)	25/268	0.1	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

GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCT  >  NC_000913/1979141‑1979416                                                                                                                                           |                                                                                                                                          ctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaatcatgcgtgatgcagatcacacaaaacactcaattccttaacataaatgTATaa                                                                                                                                       <  1:244082‑M2/5‑1 (MQ=255)            tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggtgctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                            <  1:70365‑M2/16‑1 (MQ=255)            gctaagattttgcccaaatcgacgcgacggcactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                            <  1:209985‑M2/16‑1 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggggggggactgccattttcaataatgcgtgatgcagatcacacaaaacactcacttacttcacataaatgTATAAgccataccatcgc                                                                                                                          >  1:442196‑M2/126‑137 (MQ=255)                 tatttttcctgagtcggcgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                       <  1:261510‑M2/21‑1 (MQ=255)                   tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                     <  2:81979‑M2/23‑1 (MQ=255)                          taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagttcacacaaaactcccaattacttagcataaataTATAAGTCATACGTTTGTTTTGGGTGTAtt                                                                                                              >  2:241524‑M2/114‑143 (MQ=255)                            aatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTAGTTTTGGGTGTATTTc                                                                                                            >  1:441404‑M2/112‑143 (MQ=255)                                   gcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                     <  1:445077‑M2/39‑1 (MQ=255)                                       ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaacctgaacaa                                                                                                 >  2:321164‑M2/101‑138 (MQ=255)                                                 cactacacgcacatacaacggaggggggctgcgattttcaatagtgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                                                       >  1:69945‑M2/91‑143 (MQ=255)                                                   ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                     >  1:286544‑M2/89‑143 (MQ=255)                                                   ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                                                                     >  1:286032‑M2/89‑143 (MQ=255)                                                             atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                           >  2:310626‑M2/79‑143 (MQ=255)                                                                         ggggcttcgctttcaattaatgcgtgatgcatatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                               <  1:262235‑M2/77‑1 (MQ=255)                                                                         ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                               <  2:31897‑M2/77‑1 (MQ=255)                                                                                 gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGTTTTTTGCTTCGCTAGGAGACACAaatgagaccgaaaga                                                       >  1:279788‑M2/59‑136 (MQ=255)                                                                                   ttttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                     >  1:454082‑M2/57‑143 (MQ=255)                                                                                            atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGCTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAACggg                                            >  2:439709‑M2/48‑143 (MQ=255)                                                                                               cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                         >  2:318042‑M2/45‑143 (MQ=255)                                                                                                      cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATATTCGAAATGggacgttttgca                                  >  1:177363‑M2/38‑138 (MQ=255)                                                                                                      cagatcacacaaaacactcaattacttaacacaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                  <  1:138130‑M2/106‑1 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGCTTTTGGTGTATTTAAATCTTTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGCTGTTGCCcaca                                >  1:372824‑M2/36‑143 (MQ=255)                                                                                                          tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                              <  2:286544‑M2/110‑1 (MQ=255)                                                                                                          tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                              <  2:286032‑M2/110‑1 (MQ=255)                                                                                                              acaaaacactcaattacctaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                   <  1:488932‑M2/89‑1 (MQ=255)                                                                                                              acaaaacactcaattacctaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                   >  2:488932‑M2/30‑118 (MQ=255)                                                                                                                aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAACTGGGATGTTGCGCACCGTCAAAAt                        >  2:479667‑M2/28‑143 (MQ=255)                                                                                                                                 aacataaatgTATAAGTCAGACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaataat       <  2:81553‑M2/133‑1 (MQ=255)                                                                                                                                      aaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCt  <  1:383870‑M2/138‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCT  >  NC_000913/1979141‑1979416

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