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	42 [0]	[0] 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	0 (0.000)	39 (0.870)	32/266	0.5	100%	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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_bme_pgi/1978363‑1978644                                                                                                                                            |                                                                                                                                               gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                       >  1:801547‑M1/1‑140 (MQ=255)     tgctcgttatctaaCATTTTGCCTAAATCGACGCATCTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACACAGTCATTTAACATAAATGtataagtcatact                                                                                                                                   <  2:182586‑M1/143‑14 (MQ=255)          cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATACTGAATGATGCAGATCACAAATAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                              <  1:120170‑M1/149‑19 (MQ=255)          cTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgt                                                                                                                              <  1:516907‑M1/149‑19 (MQ=255)              tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGCGGGCTGCGATTTTCAATAATGCGTGATGAAGATCACACAAAGCCCTCAATTACTTCACATAAATGtataagtcatacttctgttttg                                                                                                                          >  2:328727‑M1/1‑127 (MQ=255)               cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                         <  1:415331‑M1/149‑24 (MQ=255)                 aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataaatcatacttctgttttgggt                                                                                                                       >  2:6319‑M1/1‑124 (MQ=255)                 aaGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGAGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATCAATAtatatgtaattctcttgttttttgt                                                                                                                       >  2:290157‑M1/1‑124 (MQ=255)                    aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgta                                                                                                                    >  2:689276‑M1/1‑121 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             >  2:405055‑M1/1‑114 (MQ=255)                           cTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaat                                                                                                             >  1:552987‑M1/1‑114 (MQ=255)                            tAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatc                                                                                                            <  1:337054‑M1/149‑37 (MQ=255)                               aTCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                         <  1:582439‑M1/149‑40 (MQ=255)                                tCGACGCCACTGTACGCGTCAATACAAGCACATACAAGTGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATTAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                        <  2:58834‑M1/149‑41 (MQ=255)                                tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgtt                                                                                                        >  1:35728‑M1/1‑109 (MQ=255)                                            tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                            >  1:346093‑M1/1‑97 (MQ=255)                                            tACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttc                                                                                            >  1:288837‑M1/1‑97 (MQ=255)                                               tCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTAAATAATGCGTGATGCAGATCACACAAAACACTAAATTACTTAACATAAATGtataagtaataattttgttttgggtgtatgtcaatctgttaaaaaggttgtcgct                                                                                         >  1:91698‑M1/1‑94 (MQ=255)                                                   cACTACACGCACATACAACGGAGGGTGGCTGCGATTTTTAATAATGCTTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                     <  2:146049‑M1/149‑60 (MQ=255)                                                      tACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctag                                                                                  >  2:322622‑M1/1‑87 (MQ=255)                                                                              gctactattttcaATAATACGAGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggggtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  2:494886‑M1/143‑87 (MQ=255)                                                                              gCTGCGACTTTCAATAAGGCGTGATGCAGTTTACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataa                                                          <  1:83200‑M1/149‑87 (MQ=255)                                                                                tGCGATTTTTAATAAAGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                        <  2:336309‑M1/149‑89 (MQ=255)                                                                                tGCGATTTTCAATAATGTGTGATGCAGATCACATAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctaacaagcaaaaatgaaacaggaataatc                                                        <  2:266806‑M1/149‑89 (MQ=255)                                                                                  cGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcga                                                      >  1:704305‑M1/1‑59 (MQ=255)                                                                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                   <  1:550958‑M1/149‑94 (MQ=255)                                                                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                   <  2:304630‑M1/149‑94 (MQ=255)                                                                                         cAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcagtctgttacaaagtttttcgctacgctagtaagcaaaaatgaaacaggaataatcgaaatggga                                               >  1:743008‑M1/1‑52 (MQ=255)                                                                                               tGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgc                                         <  2:346093‑M1/149‑104 (MQ=255)                                                                                                 cGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgc                                       <  2:229513‑M1/149‑106 (MQ=255)                                                                                                  gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                      >  1:361600‑M1/1‑43 (MQ=255)                                                                                                   tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcac                                     >  2:117374‑M1/1‑42 (MQ=255)                                                                                                      tGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagt                                  <  1:689276‑M1/149‑111 (MQ=255)                                                                                                                acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                        >  1:458393‑M1/1‑29 (MQ=255)                                                                                                                  aaaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgcaatggcatgttgcgcacagtcaaaataacttc                      >  2:669549‑M1/1‑27 (MQ=255)                                                                                                                   aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatctcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   <  1:782740‑M1/135‑110 (MQ=255)                                                                                                                   aaaCACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatctcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacag                                   >  2:782740‑M1/1‑26 (MQ=255)                                                                                                                         tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa               <  1:405055‑M1/149‑130 (MQ=255)                                                                                                                                cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat        <  1:585852‑M1/149‑137 (MQ=255)                                                                                                                                     cATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgct   <  2:35728‑M1/149‑142 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:680958‑M1/1‑7 (MQ=255)                                                                                                                                      aTAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgcta  >  1:410930‑M1/1‑7 (MQ=255)                                                                                                                                            |                                                                                                                                               GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGC  >  NC_000913_3_bme_pgi/1978363‑1978644

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