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–1979270	1979278	9–776	51 [0]	[0] 51	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)	49 (1.150)	42/268	0.0	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

TGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979152‑1979417                                                                                                                                |                                                                                                                                           tctaagacttcgcctacatcgacgcagctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                             <  2:238738‑M2/16‑1 (MQ=255) gctaagatttttcctaattcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                             <  2:639728‑M2/16‑1 (MQ=255) gctaagatgtttcctaaatcgacgcgactgtactcgtcactacaggcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                             <  2:467948‑M2/16‑1 (MQ=255)   taagatttttcctagatcgacgcaactgttctcgtcactacacgcacatacagcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                           <  2:455694‑M2/18‑1 (MQ=255)   taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                           >  2:738514‑M2/126‑143 (MQ=255)    aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                                                          >  1:380437‑M2/125‑143 (MQ=255)       ttttttcccaagtcgacgcaactgctctcgtcactacacgcacatacaacggagggtggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                       <  2:141420‑M2/22‑1 (MQ=255)       atttttcctaaatcgacgcaactgtactcgtcactacacgcacataaaacggaggggggctgcgattttcaagaatacgtgaagcagatcacacaaaacactcaattactcaacataactgTCGAAGacaaacctttcatttc                                                                                                                       >  2:47109‑M2/122‑137 (MQ=255)       aattttcctatatcgacgcaactgtactcgtcactacacgcacctacaacggaggggggctgccattttcaataatgcgtgatgcagatcacacaaaacatttaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                       <  2:52574‑M2/22‑1 (MQ=255)              ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                <  1:500216‑M2/29‑1 (MQ=255)               taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTGGTTTTCGGTGTAtt                                                                                                               >  2:287146‑M2/114‑143 (MQ=255)                        gcaactgtgctcgtcactacccgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                                                                      <  2:272559‑M2/39‑1 (MQ=255)                              gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                <  1:732543‑M2/45‑1 (MQ=255)                                    ctcaccccacgcacctacaacggcgcggggctgccattttcaatactgcgtgatgcagatcacacaaaacactcaattacttaacacaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                          <  1:186950‑M2/51‑1 (MQ=255)                                     tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagagcacacaaaacactcaattacttaacatacatgTATCAGTCATACTTTTGATCAGAGTGTCTTTCAATATGTTAAACAGTTATTc                                                                                         >  2:419599‑M2/92‑143 (MQ=255)                                      cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatgatacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAACAGTTTTTcg                                                                                        <  2:470991‑M2/53‑1 (MQ=255)                                               cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTATTGTTTTGGGTGTATTTCAATCTGTTCAAAAGTTTTTCACTACACTag                                                                               >  2:43664‑M2/82‑143 (MQ=255)                                                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTGTTTGTCTACGCTagcaag                                                                           >  1:20447‑M2/78‑143 (MQ=255)                                                    acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGTTTTTCGCTACGCGagcacgc                                                                          >  2:499289‑M2/77‑143 (MQ=255)                                                      aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCtagcaagcat                                                                        >  1:490681‑M2/75‑142 (MQ=255)                                                      aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                        >  1:490806‑M2/75‑143 (MQ=255)                                                         ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                     <  1:678453‑M2/72‑1 (MQ=255)                                                           aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                   <  1:645643‑M2/74‑1 (MQ=255)                                                              ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                                <  1:500396‑M2/77‑1 (MQ=255)                                                               gggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                                                               >  2:479938‑M2/66‑143 (MQ=255)                                                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattactttacataaatgTATCAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTATTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                           >  1:216601‑M2/62‑143 (MQ=255)                                                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                           >  1:594452‑M2/62‑143 (MQ=255)                                                                      gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                        >  2:256365‑M2/59‑143 (MQ=255)                                                                            caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa                                                  <  1:459799‑M2/91‑1 (MQ=255)                                                                                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg                                              >  1:478009‑M2/49‑143 (MQ=255)                                                                                  tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                            <  1:79488‑M2/97‑1 (MQ=255)                                                                                  tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                            >  1:666797‑M2/47‑143 (MQ=255)                                                                                   gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                           >  2:485921‑M2/46‑143 (MQ=255)                                                                                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAGGAATAATCGAAATGGGATGTTgcgc                                    >  1:756044‑M2/39‑143 (MQ=255)                                                                                           cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                   <  1:9550‑M2/106‑1 (MQ=255)                                                                                             gatcacacaaatctctgaattacttaacttagatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                                 <  2:540123‑M2/108‑1 (MQ=255)                                                                                              agcacacaaaccactcaattacttatcataaatgTATAAGTCATACTTTTGTATTGAATGAAATTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAg                                <  1:780704‑M2/109‑1 (MQ=255)                                                                                                cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                              >  1:453769‑M2/33‑143 (MQ=255)                                                                                                 caacaaaacattcaattacttagcatatctcTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                             <  2:780120‑M2/112‑1 (MQ=255)                                                                                                 acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAACTGGGCTGTTGCGCACAGTCa                             >  1:695396‑M2/32‑143 (MQ=255)                                                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                        >  2:28993‑M2/28‑128 (MQ=255)                                                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt                                        <  1:28993‑M2/101‑1 (MQ=255)                                                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAt                         >  1:654820‑M2/28‑143 (MQ=255)                                                                                                            tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc                  >  1:657978‑M2/21‑143 (MQ=255)                                                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAGCAGGAAGAAGCGAAATGGGATGTTGCGGACAGTCAAAATAACTAACCg                >  1:418051‑M2/19‑143 (MQ=255)                                                                                                              aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATGTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCg                <  2:288178‑M2/125‑1 (MQ=255)                                                                                                                  acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaaa            <  1:738514‑M2/129‑1 (MQ=255)                                                                                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAACGTAaata          >  1:100664‑M2/13‑143 (MQ=255)                                                                                                                    ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata          >  2:524119‑M2/13‑143 (MQ=255)                                                                                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat     >  1:107515‑M2/8‑143 (MQ=255)                                                                                                                            aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTg  <  1:521516‑M2/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                           TGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTG  >  NC_000913/1979152‑1979417

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