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

TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978406‑1978601                                                                                                 |                                                                                                    tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtat                                                                                                  >  2:908618‑M1/1‑97 (MQ=255)  nctCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtata                                                                                                 >  2:325860‑M1/2‑96 (MQ=255)   cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataa                                                                                                >  2:902288‑M1/1‑95 (MQ=255)     cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagt                                                                                              >  2:445850‑M1/1‑93 (MQ=255)      gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtc                                                                                             <  1:1591807‑M1/100‑9 (MQ=255)       tCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                            >  2:1346120‑M1/1‑91 (MQ=255)        cACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                           >  2:515346‑M1/1‑90 (MQ=255)          cTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatac                                                                                         <  2:576421‑M1/100‑13 (MQ=255)             cacGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttt                                                                                      >  1:890263‑M1/1‑85 (MQ=255)               cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                    <  2:144694‑M1/100‑18 (MQ=255)                gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                   >  2:1553752‑M1/1‑82 (MQ=255)                  acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttn                                                                                 <  1:325857‑M1/100‑21 (MQ=255)                   caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttt                                                                                >  2:311605‑M1/1‑79 (MQ=255)                    aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                               >  2:1414233‑M1/1‑78 (MQ=255)                    aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                               <  1:210095‑M1/100‑23 (MQ=255)                     tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  1:1432732‑M1/1‑77 (MQ=255)                     tACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                              >  2:1156845‑M1/1‑77 (MQ=255)                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:1531048‑M1/100‑25 (MQ=255)                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  2:36876‑M1/100‑25 (MQ=255)                      aCAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             <  1:517978‑M1/100‑25 (MQ=255)                       cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  1:38897‑M1/1‑75 (MQ=255)                       cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggt                                                                            >  2:70705‑M1/1‑75 (MQ=255)                       cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttggg                                                                             >  2:870996‑M1/1‑75 (MQ=255)                        cgatctggggggctgcgaTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTGACATAAATGtataagtcatacttttgttttgggtg                                                                           <  2:413187‑M1/94‑27 (MQ=255)                         aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgt                                                                          <  1:146281‑M1/100‑28 (MQ=255)                              ggggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTGACATAAATGtataagtcatacttttgttttgggtgagatcg                                                                     >  1:413183‑M1/1‑68 (MQ=255)                               gggggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                    <  1:1064170‑M1/100‑34 (MQ=255)                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               <  2:752246‑M1/100‑39 (MQ=255)                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               >  1:1541973‑M1/1‑62 (MQ=255)                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctg                                                               <  2:357323‑M1/100‑39 (MQ=255)                                    cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                 <  1:1433721‑M1/98‑37 (MQ=255)                                        gATTTTTAATAATGCGTGATGCAGATCACACAAAACACTCATTTATTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtta                                                            <  1:1376697‑M1/99‑42 (MQ=255)                                            ttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaag                                                       >  2:669320‑M1/1‑54 (MQ=255)                                                ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                   >  1:547736‑M1/1‑50 (MQ=255)                                                      cGTGATGCAGATCACACAAAACACTCAATTACTTAACATATATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                             <  1:1266301‑M1/100‑57 (MQ=255)                                                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                          <  1:1370695‑M1/100‑60 (MQ=255)                                                         gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                              >  1:91536‑M1/1‑41 (MQ=255)                                                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgct                                         >  2:235557‑M1/1‑40 (MQ=255)                                                          aTGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtcaaaaagtttttcgctacgct                                         <  1:977675‑M1/100‑61 (MQ=255)                                                              aGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagca                                     <  2:890270‑M1/100‑65 (MQ=255)                                                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                  <  1:1159085‑M1/100‑68 (MQ=255)                                                                 tCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                  <  1:814147‑M1/100‑68 (MQ=255)                                                                  catacaaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagca                                 <  1:1414219‑M1/97‑69 (MQ=255)                                                                    cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaa                               >  1:619600‑M1/1‑30 (MQ=255)                                                                     acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaa                              >  1:926516‑M1/1‑29 (MQ=255)                                                                        aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatg                           <  2:1536617‑M1/100‑75 (MQ=255)                                                                          acacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                         >  1:557559‑M1/1‑24 (MQ=255)                                                                           cacTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaa                        >  1:418442‑M1/1‑23 (MQ=255)                                                                            acTTAATTACTTAACATAAATGtataagtcatactttggttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                         <  1:678082‑M1/98‑77 (MQ=255)                                                                             cTCAATTACTTAACATAAATGtataagttatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                       <  2:790804‑M1/99‑79 (MQ=255)                                                                               cAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                    <  1:511623‑M1/100‑82 (MQ=255)                                                                                  ttACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaat                 <  2:414468‑M1/100‑85 (MQ=255)                                                                                      ttAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc             <  2:848899‑M1/100‑89 (MQ=255)                                                                                          cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttttttatgttatcaagcaaaaatgaaacaggaataatcgaaa         >  1:354181‑M1/1‑8 (MQ=255)                                                                                          cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa         >  2:1154652‑M1/1‑8 (MQ=255)                                                                                              aaTGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  1:1346106‑M1/100‑97 (MQ=255)                                                                                              aaTGtataagtcatacttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg     <  1:1508502‑M1/100‑97 (MQ=255)                                                                                                 gtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  >  1:176460‑M1/1‑1 (MQ=255)                                                                                                 gtataagtcatacttttgttttgggtgtatttcaatatgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg  <  2:1522684‑M1/100‑100 (MQ=255)                                                                                                 |                                                                                                    TACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCG  >  NC_000913/1978406‑1978601

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