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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	1978503–1979269	1979278	10–776	74 [0]	[0] 74	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)	59 (0.800)	44/162	NT	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

CAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979184‑1979378                                                                                                |                                                                                                    ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                 >  1:2159687‑M2/96‑100 (MQ=255)         tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          <  1:979711‑M2/12‑1 (MQ=255)         tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                         >  2:208106‑M2/88‑100 (MQ=255)              gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                    <  1:494589‑M2/18‑1 (MQ=255)              gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                    <  2:1088847‑M2/18‑1 (MQ=255)              gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                    <  2:404242‑M2/18‑1 (MQ=255)               cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   <  2:375802‑M2/19‑1 (MQ=255)                acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   <  2:454084‑M2/19‑1 (MQ=255)                  atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTAATACttttgttttg                                                                                >  2:1742851‑M2/79‑100 (MQ=255)                   tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  1:1636054‑M2/78‑100 (MQ=255)                    acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  1:1036906‑M2/24‑1 (MQ=255)                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  1:1176751‑M2/76‑100 (MQ=255)                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  1:1907740‑M2/76‑100 (MQ=255)                        cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                          >  1:1095658‑M2/73‑100 (MQ=255)                         ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                         <  2:482610‑M2/29‑1 (MQ=255)                           agaggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcatttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                       <  2:1915514‑M2/31‑1 (MQ=255)                            ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAttt                                                                       >  2:496666‑M2/69‑99 (MQ=255)                              ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                    >  2:407214‑M2/67‑100 (MQ=255)                                 gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                   >  1:1276736‑M2/64‑98 (MQ=255)                                 gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                 >  2:1428417‑M2/64‑100 (MQ=255)                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                               <  1:893445‑M2/39‑1 (MQ=255)                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGt                                                               <  1:1837355‑M2/39‑1 (MQ=255)                                    gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              <  2:1588906‑M2/40‑1 (MQ=255)                                       attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                           >  1:1177390‑M2/58‑100 (MQ=255)                                       attttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaa                                                           >  1:58988‑M2/58‑100 (MQ=255)                                        ttttcaataatgcttgaggcagatcacacaaaacacttattttcttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              <  2:990700‑M2/40‑1 (MQ=255)                                           tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                       <  1:1726287‑M2/47‑1 (MQ=255)                                                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  <  2:655728‑M2/52‑1 (MQ=255)                                                aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  <  2:1374554‑M2/52‑1 (MQ=255)                                                 atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 >  2:1387‑M2/48‑100 (MQ=255)                                                 atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  >  2:1963304‑M2/48‑99 (MQ=255)                                                    cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCCTTTTGTTAAACAGTTTTTcgct                                               <  2:1303780‑M2/55‑1 (MQ=255)                                                    cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                              >  1:1121035‑M2/45‑100 (MQ=255)                                                      tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATGTGTTAAAAAGTTTTTcgctacg                                            <  2:1302128‑M2/58‑1 (MQ=255)                                                      tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                  >  2:567155‑M2/43‑94 (MQ=255)                                                        atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                          <  1:1424182‑M2/60‑1 (MQ=255)                                                         tgcagatcacacaaaacactcaatttcttaacataaatgTATATGTAATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAATTATTcgctacact                                          <  2:967197‑M2/60‑1 (MQ=255)                                                         tgcagatcacacaaaacactaaattacttaacataaatgTATAAGTCATACTTTTGTTATTGGTGTATTACAATCTGTTAAAAAATTATTcgttacgcta                                         >  2:1563709‑M2/40‑100 (MQ=255)                                                          gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                        <  2:1988414‑M2/62‑1 (MQ=255)                                                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTagca                                      >  2:790222‑M2/37‑100 (MQ=255)                                                            agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGCATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                      <  1:150012‑M2/64‑1 (MQ=255)                                                              atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                    <  1:380495‑M2/66‑1 (MQ=255)                                                                cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                  <  2:360731‑M2/68‑1 (MQ=255)                                                                 acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                 <  1:1524882‑M2/69‑1 (MQ=255)                                                                 acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGATAAAAAGTTTTTCGCTACGCTGgcaagcaa                                 >  1:1276113‑M2/32‑100 (MQ=255)                                                                  cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:1139340‑M2/31‑100 (MQ=255)                                                                  cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                <  2:1992961‑M2/70‑1 (MQ=255)                                                                  cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:566269‑M2/31‑100 (MQ=255)                                                                   acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                               >  1:1776509‑M2/30‑100 (MQ=255)                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             >  2:717823‑M2/28‑100 (MQ=255)                                                                     aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  2:1342518‑M2/73‑1 (MQ=255)                                                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                              >  2:1082562‑M2/27‑98 (MQ=255)                                                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            >  2:594442‑M2/27‑100 (MQ=255)                                                                      aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            >  2:1930634‑M2/27‑100 (MQ=255)                                                                        acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                          <  1:1428427‑M2/76‑1 (MQ=255)                                                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         >  1:1300499‑M2/24‑100 (MQ=255)                                                                         cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTCTTTCGCTACGCTAGCAAGCAAAAATGaaa                         >  2:1346068‑M2/24‑100 (MQ=255)                                                                           ctcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTAAAAAAGTTTTTCGATACGCTAGCAAGCAAAAATGAAAc                        <  2:306234‑M2/78‑1 (MQ=255)                                                                             caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     >  2:1245952‑M2/20‑100 (MQ=255)                                                                             caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  1:1798329‑M2/81‑1 (MQ=255)                                                                                   cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  2:528302‑M2/14‑100 (MQ=255)                                                                                     taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg             <  2:2159672‑M2/89‑1 (MQ=255)                                                                                     taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg             <  1:1205743‑M2/89‑1 (MQ=255)                                                                                      aacataaatgTATAAGTCATATTTTTTTTTTGGGTGTATTTAAATCTATTAAAAAGTTTTTCGTTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa            <  2:1078549‑M2/90‑1 (MQ=255)                                                                                      aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa            <  1:2054304‑M2/90‑1 (MQ=255)                                                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:1489057‑M2/9‑100 (MQ=255)                                                                                        cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:956410‑M2/9‑100 (MQ=255)                                                                                          taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt         >  1:153316‑M2/7‑99 (MQ=255)                                                                                            aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg      >  2:1816726‑M2/5‑100 (MQ=255)                                                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    >  2:1139153‑M2/3‑100 (MQ=255)                                                                                              tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAGCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa     >  1:1579480‑M2/3‑99 (MQ=255)                                                                                               gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   >  2:1729961‑M2/2‑100 (MQ=255)                                                                                                tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  2:105008/100‑1 (MQ=255)                                                                                                tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  1:1970184/100‑1 (MQ=255)                                                                                                |                                                                                                    CAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979184‑1979378

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