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

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913_3_bme_pgi	1978503	1979278	776	41 [2]	[1] 42	insB1–insA	insB1,insA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913_3_bme_pgi	= 1978502	2 (0.040)	36 (0.860)	32/270	0.3	97.4%	intergenic (‑305/+16)	flhD/insB1	flagellar class II regulon transcriptional activator, with FlhC/IS1 transposase B
	?	NC_000913_3_bme_pgi	1979279 =	0 (0.000)					intergenic (‑64/‑474)	insA/uspC	IS1 repressor TnpA/universal stress protein

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_bme_pgi/1978363‑1978649                                                                                                                                            |                                                                                                                                                    gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATCAATGtataagtca                                                                                                                                            >  2:422120‑M1/1‑140 (MQ=255)        cacTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttt                                                                                                                                     <  1:494976‑M1/149‑17 (MQ=255)               cTAAGATTGTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAATCACTTAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                              <  2:303636‑M1/149‑24 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatactttagatttggg                                                                                                                             >  1:559341‑M1/1‑125 (MQ=255)                         tCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttca                                                                                                                    <  2:699485‑M1/149‑34 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATCACGCGTGAGGCAGAGCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                  >  2:510034‑M1/1‑114 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGCTCAAACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                                  >  1:372055‑M1/1‑114 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTACAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAAAATAAATGtataagtcatagttttgttttggtggtatttcaat                                                                                                                  >  2:320260‑M1/1‑114 (MQ=255)                            tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                                 >  2:635684‑M1/1‑113 (MQ=255)                             aaaTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatct                                                                                                                <  1:411896‑M1/149‑38 (MQ=255)                               ctcGACGCAACTCTGCTCTCCACAACACGCACATACAACTGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                              <  2:411545‑M1/148‑40 (MQ=255)                                tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatacttttgttttgggtgtatctcaatctgtt                                                                                                             >  2:727059‑M1/1‑109 (MQ=255)                                            tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGAGTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                                 >  2:709219‑M1/1‑97 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgct                                                                                              >  1:554377‑M1/1‑94 (MQ=255)                                                cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                             <  1:24803‑M1/149‑57 (MQ=255)                                                 gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctac                                                                                            >  1:732335‑M1/1‑92 (MQ=255)                                                  tCACTACACGCCCATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttaaatctgttcaaacgtttttcgctccg                                                                                           >  2:271989‑M1/1‑91 (MQ=255)                                                      ttcacgcACATACAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTGAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                       <  2:207723‑M1/147‑63 (MQ=255)                                                      tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcactctgttgtaatgtttttcgctacgctag                                                                                       >  1:561036‑M1/1‑87 (MQ=255)                                                             acaTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaa                                                                                <  2:732335‑M1/149‑70 (MQ=255)                                                                tACAGCGGAGGGGGGCTGCTATTTACAATAATGCGTGCTGCAGATCACATAAAACACTTTATTACTTAACATAAATGtataagtcatatttttgttttgagtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaa                                                                             <  2:118075‑M1/149‑73 (MQ=255)                                                                    aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                         >  2:721164‑M1/1‑73 (MQ=255)                                                                                        tCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                     >  1:268570‑M1/1‑53 (MQ=255)                                                                                            taatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt                                                 >  1:633877‑M1/1‑49 (MQ=255)                                                                                              atGCGTGAAGCAGATCACAAAAAACGCTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttg                                               <  1:422120‑M1/149‑103 (MQ=255)                                                                                                   tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtctttcaatctgctgaaaagttgttcgctacgctagcaagcaaaaataaaacagtaataatcgaaatgggatgttgccgac                                          >  2:277675‑M1/1‑42 (MQ=255)                                                                                                    gATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaca                                         <  1:267839‑M1/149‑109 (MQ=255)                                                                                                      tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                       >  2:108923‑M1/1‑39 (MQ=255)                                                                                                           aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                  >  1:46961‑M1/1‑34 (MQ=255)                                                                                                           aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                  >  2:686564‑M1/1‑34 (MQ=255)                                                                                                               cacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattttaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                              >  2:741470‑M1/1‑30 (MQ=255)                                                                                                               cacaAAACACTCAATGACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataac                              <  2:414041‑M1/149‑120 (MQ=255)                                                                                                                       acTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggattttgcgcacagtcaaaataactccccttc                      >  2:308938‑M1/1‑22 (MQ=255)                                                                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcagtctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatggaaatgggatgttgcgcacagtcaaaataactcaccgtaaata                  >  1:156913‑M1/1‑18 (MQ=255)                                                                                                                           aaTTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaacataactcaccgtaaata                  >  1:425628‑M1/1‑18 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaataatcatctg          >  2:237358‑M1/1‑10 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                       >  2:388923‑M1/1‑7 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg                                                       <  1:388923‑M1/101‑95 (MQ=255)                                                                                                                                           tGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaacaaggaataatcaaaatgggatgttgcgcacagtcaaaatgactcaccgtaaataatcatctgctactaat  >  2:519595‑M1/1‑2 (MQ=255)                                                                                                                                           cgggTAATGCCTGCAACTTATGGATAGTGTTTTATGTCGAGATAGTGTGTGATGTCTTTGTCATGCAGCTCCACTGATTGTGAGAGCGACAGGGACATCCGGATCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTCTGCCGCTCaaa  <  2:448619/148‑2 (MQ=2)                                                                                                                                           cgggTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTAATCCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCCGGTTATGCCGCTCAAt  >  1:215926/2‑149 (MQ=2)                                                                                                                                            |                                                                                                                                                    GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAAT  >  NC_000913_3_bme_pgi/1978363‑1978649

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