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–1979269	1979278	10–776	59 [0]	[2] 63	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)	53 (0.880)	42/182	0.5	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

AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979182‑1979378                                                                                                  |                                                                                                    tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                   >  2:908618‑M2/98‑100 (MQ=255)  nctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                  >  2:325860‑M2/97‑100 (MQ=255)   ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                 >  2:902288‑M2/96‑100 (MQ=255)     cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                               >  2:445850‑M2/94‑100 (MQ=255)      gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                              <  1:1591807‑M2/8‑1 (MQ=255)       tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                             >  2:1346120‑M2/92‑100 (MQ=255)        cactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAt                                                                                            >  2:515346‑M2/91‑100 (MQ=255)          ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                          <  2:576421‑M2/12‑1 (MQ=255)             cacgcacatacaacggaggggggctgcgattttcaataatgcgtgaggcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                       >  1:890263‑M2/86‑100 (MQ=255)               cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                     <  2:144694‑M2/17‑1 (MQ=255)                gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgt                                                                                    >  2:1553752‑M2/83‑100 (MQ=255)                  acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATActtttgttn                                                                                  <  1:325857‑M2/20‑2 (MQ=255)                   catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                 >  2:311605‑M2/80‑100 (MQ=255)                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                <  1:210095‑M2/22‑1 (MQ=255)                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                >  2:1414233‑M2/79‑100 (MQ=255)                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  2:1156845‑M2/78‑100 (MQ=255)                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  1:1432732‑M2/78‑100 (MQ=255)                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  2:36876‑M2/24‑1 (MQ=255)                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  1:517978‑M2/24‑1 (MQ=255)                      acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  2:1531048‑M2/24‑1 (MQ=255)                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              >  2:870996‑M2/76‑99 (MQ=255)                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  1:38897‑M2/76‑100 (MQ=255)                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             >  2:70705‑M2/76‑100 (MQ=255)                        cgatctggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttgacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                            <  2:413187‑M2/26‑1 (MQ=255)                         acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtgt                                                                           <  1:146281‑M2/27‑1 (MQ=255)                              ggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttgacataaatgTATAAGTCATACTTTTGttttgggtgagatcg                                                                      >  1:413183‑M2/69‑94 (MQ=255)                               gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                     <  1:1064170‑M2/33‑1 (MQ=255)                                    ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                  <  1:1433721‑M2/36‑1 (MQ=255)                                    ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                <  2:357323‑M2/38‑1 (MQ=255)                                    ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                <  2:752246‑M2/38‑1 (MQ=255)                                    ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                >  1:1541973‑M2/63‑100 (MQ=255)                                        gatttttaataatgcgtgatgcagatcacacaaaacactcatttatttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                             <  1:1376697‑M2/41‑1 (MQ=255)                                            ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                        >  2:669320‑M2/55‑100 (MQ=255)                                                ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                    >  1:547736‑M2/51‑100 (MQ=255)                                                      cgtgatgcagatcacacaaaacactcaattacttaacatatatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                              <  1:1266301‑M2/56‑1 (MQ=255)                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgc                                           <  1:1370695‑M2/59‑1 (MQ=255)                                                         gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                               >  1:91536‑M2/42‑96 (MQ=255)                                                          atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                          >  2:235557‑M2/41‑100 (MQ=255)                                                          atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTCAAAAAGTTTTTcgctacgct                                          <  1:977675‑M2/60‑1 (MQ=255)                                                              agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                      <  2:890270‑M2/64‑1 (MQ=255)                                                                 tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                   <  1:1159085‑M2/67‑1 (MQ=255)                                                                 tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                   <  1:814147‑M2/67‑1 (MQ=255)                                                                  catacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                  <  1:1414219‑M2/68‑1 (MQ=255)                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:619600‑M2/31‑100 (MQ=255)                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                               >  1:926516‑M2/30‑100 (MQ=255)                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            <  2:1536617‑M2/74‑1 (MQ=255)                                                                          acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                          >  1:557559‑M2/25‑100 (MQ=255)                                                                           cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                         >  1:418442‑M2/24‑100 (MQ=255)                                                                            acttaattacttaacataaatgTATAAGTCATACTTTGGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                          <  1:678082‑M2/76‑1 (MQ=255)                                                                             ctcaattacttaacataaatgTATAAGTTATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        <  2:790804‑M2/78‑1 (MQ=255)                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  1:511623‑M2/81‑1 (MQ=255)                                                                                  ttacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                  <  2:414468‑M2/84‑1 (MQ=255)                                                                                      ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              <  2:848899‑M2/88‑1 (MQ=255)                                                                                        aacATAAATGTATAAGTCATACTTTAGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAAAAAAAAAAAAAAAATGAAACAGGAATAATCg             <  1:790546/97‑1 (MQ=255)                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTTTTTATGTTATCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:354181‑M2/9‑100 (MQ=255)                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  2:1154652‑M2/9‑100 (MQ=255)                                                                                          cATAAATGTATAAGTCATACTTTTTTTTTTTTTTTATTTCAATATGTTAAAAAGTTTTTCTCTACGTTAGCAATCATAAATGAAACAGTAATATTCGaa           >  1:257110/1‑99 (MQ=255)                                                                                              aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg      <  1:1346106‑M2/96‑1 (MQ=255)                                                                                              aatgTATAAGTCATACTTTTGTTTTGAGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg      <  1:1508502‑M2/96‑1 (MQ=255)                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   >  1:176460‑M2/2‑100 (MQ=255)                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   <  2:1522684‑M2/99‑1 (MQ=255)                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  2:802303/100‑1 (MQ=255)                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  >  2:18691/1‑100 (MQ=255)                                                                                                  |                                                                                                    AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979182‑1979378

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