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	64 [1]	[0] 63	insB1–insA	insB1, insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1978502	1 (0.020)	61 (0.980)	44/268	0.3	99.2%	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

GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCAT  >  NC_000913/1979141‑1979414                                                                                                                                           |                                                                                                                                        ctttacagttctctaagatctttcctaactcgacgcaactgttctcttcactacacgcacatacaacggagggtggctgcgattttcaataatccgtgatgcagatcacataaaacactcaattacttaacataaatgTATaa                                                                                                                                     <  2:681611‑M2/5‑1 (MQ=255)  tttacacttatctaagatttttcctaaatcgacgcacctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaatactgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                    >  1:744244‑M2/138‑143 (MQ=255)  tttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataattcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAg                                                                                                                                    >  2:110002‑M2/138‑143 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                        >  1:788003‑M2/126‑143 (MQ=255)              taagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattgtcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                                                        >  1:1138496‑M2/126‑143 (MQ=255)               aagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaccataaatgTATAAGTCATACttttgtt                                                                                                                       >  1:367675‑M2/125‑143 (MQ=255)                   tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                   <  1:91595‑M2/23‑1 (MQ=255)                   tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                   <  2:140098‑M2/23‑1 (MQ=255)                   ttttgcctaaatcgacgcgactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                   <  2:1147082‑M2/23‑1 (MQ=255)                   tttctcctgaaccgacgcaactgtactcgtcactccacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                   <  2:1147038‑M2/23‑1 (MQ=255)                   ctcttcctagatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                   <  2:896261‑M2/23‑1 (MQ=255)                     tttcctaaatcgacgcaactgtgttcggcacgacacgcacatgccccgcatgggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                 <  1:1025873‑M2/25‑1 (MQ=255)                     tgtcctagagcggcgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                                                                 <  1:468576‑M2/25‑1 (MQ=255)                             atcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                                                         >  1:346154‑M2/111‑143 (MQ=255)                                 acgcaacggtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                                                     >  2:470529‑M2/107‑143 (MQ=255)                                  cgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                                                    <  2:271983‑M2/38‑1 (MQ=255)                                       ctgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                                                               <  2:964576‑M2/43‑1 (MQ=255)                                        tgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTAGCTTTGGGTATATTTCAATCTGTTaaaa                                                                                              >  2:67382‑M2/100‑143 (MQ=255)                                             tcgtcagtacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAACttt                                                                                         >  2:201050‑M2/95‑143 (MQ=255)                                             tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCAGTCATACTTTTGTTTTGGGAGTATTTCAGTCTTTTACAAAGttt                                                                                         >  2:436346‑M2/95‑143 (MQ=255)                                                tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                                                      <  1:880011‑M2/52‑1 (MQ=255)                                                  actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                                                    >  1:1092300‑M2/90‑143 (MQ=255)                                                           acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                           <  2:1084710‑M2/63‑1 (MQ=255)                                                           acatacaacggaggggggcggcgattttcaataatgcgtgatgcagatcacacaaaaccctcaataagtggacataagtgTTTAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                           <  1:266571‑M2/63‑1 (MQ=255)                                                             atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattccttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                         <  2:346154‑M2/65‑1 (MQ=255)                                                              tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaccactcaattacttaacataaatgTATAAGTCATACTTTTGTACTGGGTGTATTTCCACCAGTTAAAAAATTTTTCGCTACGCTagcaag                                                                        >  2:10544‑M2/78‑143 (MQ=255)                                                              tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                        <  1:504871‑M2/66‑1 (MQ=255)                                                              gacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                                                        <  2:217456‑M2/66‑1 (MQ=255)                                                                    ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCGATCTGTGAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                  >  1:1085962‑M2/72‑143 (MQ=255)                                                                      aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCACAAATg                                                                >  2:333589‑M2/70‑143 (MQ=255)                                                                         ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                                                             <  1:234984‑M2/77‑1 (MQ=255)                                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAACAAGTTTTTCGCTACGCTAGCAAGCAAACATGAAACAGGa                                                        >  1:97380‑M2/62‑143 (MQ=255)                                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                        >  1:623820‑M2/62‑143 (MQ=255)                                                                              tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATAATTTTGTTCTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCGAGAAAAAGTGAAACAGGa                                                        >  2:922468‑M2/62‑143 (MQ=255)                                                                                cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                      >  1:278779‑M2/60‑143 (MQ=255)                                                                                 gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                                                     >  2:275681‑M2/59‑143 (MQ=255)                                                                                  attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                    >  1:1128672‑M2/58‑143 (MQ=255)                                                                                  attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                                                    >  1:473783‑M2/58‑143 (MQ=255)                                                                                    tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTAGTTTGGGGTGTATTTCAACATGTTAATAAGTGGTGCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc                                                  >  1:314692‑M2/56‑143 (MQ=255)                                                                                      tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGAAAGCAAGCAAAAATGAAACAGGAATAATCGa                                                >  1:89329‑M2/54‑143 (MQ=255)                                                                                             tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATTGGa                                         >  2:833890‑M2/47‑143 (MQ=255)                                                                                               cgtgatgaagatcacacaaaagacttaattacttagcataaatgTATAAGTCATACTTTTGTTTTGGGGGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg                                       <  2:906389‑M2/99‑1 (MQ=255)                                                                                                gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                      >  2:253886‑M2/44‑143 (MQ=255)                                                                                                  gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTg                                    >  2:884327‑M2/42‑143 (MQ=255)                                                                                                      cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGca                                >  2:590812‑M2/38‑143 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                              <  1:248143‑M2/108‑1 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                              <  2:623820‑M2/108‑1 (MQ=255)                                                                                                        gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGcaca                              >  1:1153074‑M2/36‑143 (MQ=255)                                                                                                          tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGCAATGGGATGTTGCGCACAGt                            >  2:856918‑M2/34‑143 (MQ=255)                                                                                                              acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                        >  1:490895‑M2/30‑143 (MQ=255)                                                                                                               caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTCAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAGACGCGATGTTCAGCGCAGTCaaaa                       >  2:294629‑M2/29‑143 (MQ=255)                                                                                                                  aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATaa                    >  1:935345‑M2/26‑143 (MQ=255)                                                                                                                     actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTc                 >  2:357838‑M2/23‑143 (MQ=255)                                                                                                                       tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAAtaactaac               >  2:310799‑M2/21‑143 (MQ=255)                                                                                                                       tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGCTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAc               >  2:311276‑M2/21‑143 (MQ=255)                                                                                                                        caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAcc              <  1:982235‑M2/124‑1 (MQ=255)                                                                                                                         aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAACGTTTTTCGCTACGCTAGCACGCAAAACCGAAACCGGACTCATCGAAACGGGCTGTTGCGCACAGTCCAAATAACCCACCg             >  2:558379‑M2/19‑143 (MQ=255)                                                                                                                         aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGGAATGGGATGTTGCGCACAGTCAAAATAACTCACCg             >  2:760323‑M2/19‑143 (MQ=255)                                                                                                                           ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTa           >  2:361993‑M2/17‑143 (MQ=255)                                                                                                                            tacttagtgtaaatgTATGAGTCATAATTTTGTTTTGGGTATATTTCGATCTGTTAAAAAGTTTTTCGTGACGCTAGCAAGTAAAAATGAAACAGGAATAAGCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTaa          <  2:487842‑M2/128‑1 (MQ=255)                                                                                                                               ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATACTCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata       >  2:220496‑M2/13‑143 (MQ=255)                                                                                                                               ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAAATCACAgtaaaga       >  2:1105709‑M2/13‑141 (MQ=255)                                                                                                                                    ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat  >  2:417437‑M2/8‑143 (MQ=255)                                                                                                                                           |                                                                                                                                        GGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCAT  >  NC_000913/1979141‑1979414

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