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	58 [0]	[0] 58	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)	56 (1.170)	46/266	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

CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979140‑1979419                                                                                                                                            |                                                                                                                                             tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                           >  2:483165‑M2/140‑143 (MQ=255)     tacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                                                                       <  2:833808‑M2/8‑1 (MQ=255)             tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                                                               <  1:862575‑M2/16‑1 (MQ=255)               taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttagcataactgTATAAGTCAtcctttggg                                                                                                                             >  1:314010‑M2/126‑140 (MQ=255)               taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaaaactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                             >  2:15217‑M2/126‑143 (MQ=255)               taagatgattcttaaatcgacgcaactgtactcgtcactacacgcacagacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaatctctcaattacttaacataaatgTATAAGTCATACtgttgt                                                                                                                             <  2:438060‑M2/18‑1 (MQ=255)                  gatttttcctaacgctacgctgctgtcctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                                                          <  2:868426‑M2/21‑1 (MQ=255)                   atttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                         <  1:870116‑M2/22‑1 (MQ=255)                    tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                        <  1:232017‑M2/23‑1 (MQ=255)                    tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                        <  1:102255‑M2/23‑1 (MQ=255)                          ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgggtgtgt                                                                                                                  >  2:26164‑M2/115‑141 (MQ=255)                            aagtcgacgcaacggtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                                                                <  1:499070‑M2/31‑1 (MQ=255)                                cgactcaacggttctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatacagatcacacaaaactcgtaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                            <  2:864090‑M2/35‑1 (MQ=255)                                 gacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                                                           <  2:810057‑M2/36‑1 (MQ=255)                                    gcacgtgtactcggcacgacacgcgcatacaaggggggggggctgcgattttcactaatgcgtgatgcagatcacacaaaactctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATCTCCATCTGt                                                                                                        <  2:712995‑M2/39‑1 (MQ=255)                                       actgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgccgatcacccaaaacactcaattactgaacataaatgTATAAGTGATAATTGTGTTTTGGATGTATCTccatctctcaa                                                                                                     >  2:198713‑M2/102‑138 (MQ=255)                                          gtactcgttacgacacggccatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcagttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                  <  2:27927‑M2/45‑1 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                              >  1:527694‑M2/95‑143 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                              >  1:527791‑M2/95‑143 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATATAAATCTGTTAAAAAGttt                                                                                              >  2:572605‑M2/95‑143 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTCTGGGTGTATTTCAATCTGTTAAAAAtttt                                                                                              >  2:251858‑M2/95‑143 (MQ=255)                                                gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTACAAGTCATACTTTTGTTTTGGGCCTATTTCCTTCTGCTAACAAGtcttt                                                                                            >  1:401370‑M2/93‑143 (MQ=255)                                                gtcaccacacgcacatacacagggggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                            <  1:580050‑M2/51‑1 (MQ=255)                                                     tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                                                                       <  1:287496‑M2/56‑1 (MQ=255)                                                          gcacttacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                  <  1:887203‑M2/61‑1 (MQ=255)                                                          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                  <  2:71249‑M2/61‑1 (MQ=255)                                                          gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                                                                  <  1:144596‑M2/61‑1 (MQ=255)                                                            acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                >  1:332502‑M2/81‑143 (MQ=255)                                                                  aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                                                          >  2:806024‑M2/75‑143 (MQ=255)                                                                   acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                                                         >  1:157943‑M2/74‑143 (MQ=255)                                                                    cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                                                                        <  1:125862‑M2/71‑1 (MQ=255)                                                                            ggctgcgattttcaataatgcgtgatgcagatcacacaaaactctctattacttaacatttatgTTTAATTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                <  1:883985‑M2/79‑1 (MQ=255)                                                                             gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                               >  1:425909‑M2/64‑143 (MQ=255)                                                                              ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAACAATGAAACAgg                                                              >  2:894015‑M2/63‑143 (MQ=255)                                                                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAATAGTTTTTCGCTATGCTAGCAAGCAAATATGAAACAGGa                                                             >  2:534237‑M2/62‑143 (MQ=255)                                                                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAACTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTCGCAAGCAAAAATGAAAcagac                                                             >  1:176163‑M2/62‑141 (MQ=255)                                                                                 cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTGAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                           >  2:626524‑M2/60‑143 (MQ=255)                                                                                    tttttattaatgcgtgatttagatcacacaaaacattcaatttcttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                        <  2:650762‑M2/87‑1 (MQ=255)                                                                                      ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg                                                      <  2:293079‑M2/89‑1 (MQ=255)                                                                                         aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGGTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:416606‑M2/52‑143 (MQ=255)                                                                                         aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                   >  1:641228‑M2/52‑143 (MQ=255)                                                                                             atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg                                               <  1:866264‑M2/96‑1 (MQ=255)                                                                                             atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCAAGCAAAACTGAAACAGGAATAATCGAAATggg                                               >  1:607380‑M2/48‑143 (MQ=255)                                                                                              tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAGTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCCAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                              >  1:533265‑M2/47‑143 (MQ=255)                                                                                                     tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAGCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTgcg                                       <  1:105138‑M2/104‑1 (MQ=255)                                                                                                       cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGCAACAGGAATAATCGAAATGGGATGTTGCGca                                     >  1:70473‑M2/38‑143 (MQ=255)                                                                                                          atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGCTGCGAACAg                                  >  2:850274‑M2/35‑143 (MQ=255)                                                                                                           tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAACGTTGTCCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGt                                 >  1:127839‑M2/34‑143 (MQ=255)                                                                                                               acaaaacactcgattacttaacagaaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAGCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                             <  1:200205‑M2/114‑1 (MQ=255)                                                                                                                caaaacactcaattacttaacattaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaaa                            <  1:483165‑M2/115‑1 (MQ=255)                                                                                                                 aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAGGGGATGTTGCGCACAGTCAAAAt                           >  1:290033‑M2/28‑143 (MQ=255)                                                                                                                               cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaat             <  1:367394‑M2/130‑1 (MQ=255)                                                                                                                                ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata            >  1:832726‑M2/13‑143 (MQ=255)                                                                                                                                ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTAACCTTAaata            >  2:435987‑M2/13‑143 (MQ=255)                                                                                                                                  aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAACACAACTCACCGTAaaaaat          >  2:238417‑M2/11‑143 (MQ=255)                                                                                                                                   acataaatgTATAAGTCATACTTTTGTTTTGGGGGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc         <  2:70400‑M2/134‑1 (MQ=255)                                                                                                                                         atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGc   <  1:240573‑M2/140‑1 (MQ=255)                                                                                                                                          tttataACTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAATAGGAATAATCGAAATGGGATGTTGCGCACATTCAAAATAACTCACCGTAAATAATCATCTGCt  <  1:747491/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCT  >  NC_000913/1979140‑1979419

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