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	36 [0]	[0] 37	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)	36 (0.850)	29/270	0.6	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

GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_bme_pgi/1978363‑1978641                                                                                                                                            |                                                                                                                                            gTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtca                                                                                                                                    <  2:79097‑M1/149‑10 (MQ=255)  tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                   >  1:201783‑M1/1‑139 (MQ=255)  tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATCCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcat                                                                                                                                   >  1:543318‑M1/1‑139 (MQ=255)  tCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATAAAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAATAAATTACTTAACATAAATAtataagtcat                                                                                                                                   >  2:63561‑M1/1‑139 (MQ=255)     ttACACTTATCTAAGTTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGACGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatact                                                                                                                                <  2:266900‑M1/149‑14 (MQ=255)         acTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttg                                                                                                                            >  1:91769‑M1/1‑132 (MQ=255)              tCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttg                                                                                                                       <  2:121732‑M1/149‑23 (MQ=255)               cgaagaTGTTCGCTAAATCGACGCAACTGGACTCGTCACCACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATTCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                      <  2:646538‑M1/147‑24 (MQ=255)                tAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataactcatacatttgtttttgg                                                                                                                     >  1:594861‑M1/1‑125 (MQ=255)                    aTTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtaatacttttgttttgggtgta                                                                                                                 >  1:604096‑M1/1‑121 (MQ=255)                               aTCGAAGGAACTGTACTCGGCACTACACGCACATACACCGGAGGGGGGCTGCGATTTTCAATTATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAAGGtataagtcatacttttgttttgggtgtatttcaatctgt                                                                                                      <  2:135570‑M1/149‑40 (MQ=255)                                         cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                            >  1:669046‑M1/1‑100 (MQ=255)                                              cTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgc                                                                                       >  1:173630‑M1/1‑95 (MQ=255)                                                   caccacACGCACATACAACGGAGTGGGTCTGCGAGTTGCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttcgctacgc                                                                                  <  2:436786‑M1/149‑60 (MQ=255)                                                         atcaacataaaACAGAAGGGGCCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCCATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctactctagcaa                                                                            <  2:357226‑M1/145‑66 (MQ=255)                                                         aagaacagacaACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaa                                                                            <  1:142779‑M1/145‑66 (MQ=255)                                                           gCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                          >  2:23831‑M1/1‑82 (MQ=255)                                                                      ggAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACAAAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaac                                                               <  2:610060‑M1/149‑79 (MQ=255)                                                                            gggCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttttgctacgctagcaagcaaaaatgaaacaggaat                                                         >  2:421776‑M1/1‑65 (MQ=255)                                                                                    aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa                                                 <  1:23831‑M1/149‑93 (MQ=255)                                                                                     ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                >  1:261960‑M1/1‑56 (MQ=255)                                                                                       ttgaataatgcgTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatttgttataaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatgaaaatgt                                              >  2:531951‑M1/4‑54 (MQ=255)                                                                                              atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                     >  2:293662‑M1/1‑47 (MQ=255)                                                                                              atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatc                                                     <  1:293662‑M1/135‑89 (MQ=255)                                                                                              atGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcataattttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaaccggaataatcgaactgggatggtg                                       >  1:91487‑M1/1‑47 (MQ=255)                                                                                                  gTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgca                                   >  1:671537‑M1/1‑43 (MQ=255)                                                                                                          gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                           >  2:750293‑M1/1‑35 (MQ=255)                                                                                                          gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                           >  1:351515‑M1/1‑35 (MQ=255)                                                                                                              acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaatgggatgttgcgctaagtcgaaataa                       >  1:71886‑M1/1‑31 (MQ=255)                                                                                                                acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcaacgtcaaaataact                     >  1:299299‑M1/1‑29 (MQ=255)                                                                                                                              tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc       >  1:224541‑M1/1‑15 (MQ=255)                                                                                                                               aCTTAACATAAATGtataagtcatacttttgttgtgggtgtatttcaatctgttaaaaagtttttcgctgcgctaggaggcaaaaatgaaacgggaaaaatcgaaatgggatgttgcgcacagtcaaagtaacggaccgtaaataatca      >  2:403982‑M1/1‑14 (MQ=255)                                                                                                                                cTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcat     >  1:641132‑M1/1‑13 (MQ=255)                                                                                                                                  tAACATAAATGtataaatcatacttttgttgtgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatct   <  1:63561‑M1/149‑139 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  <  2:224541‑M1/149‑140 (MQ=255)                                                                                                                                   aaCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctg  <  2:201783‑M1/149‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            GTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGC  >  NC_000913_3_bme_pgi/1978363‑1978641

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