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–1979269	1979278	10–776	38 [0]	[0] 38	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)	37 (0.940)	33/274	0.2	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

ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_bme_pgi/1978359‑1978648                                                                                                                                                |                                                                                                                                                   aTGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACGCAAAACACTCAATTACTTAACATAAATGtataa                                                                                                                                               >  2:383266‑M1/1‑144 (MQ=255)  tGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACCCAAATAAACTCGTCACTACACGCACATACAACGAAGCGGCGCTACGATTTTCAATAATACGTGATTCAGATCACACAAAACACTCAATTCCTTTCCATAAATGtataaa                                                                                                                                              >  2:42322‑M1/1‑143 (MQ=255)       cTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcata                                                                                                                                         >  1:588186‑M1/1‑138 (MQ=255)             acTTATCTAGGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACCGAGGGGGGCTGCGATTTTCAATAATGCGTCATGCAGAGAACACAAAGAAGTCCATTACTTAACATAAATGtataacacgtccttttg                                                                                                                                   >  2:2799‑M1/1‑132 (MQ=255)                 aTCAAAGAGTGTTCCGAAATCGACGCAACTGTACTCATCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACTCTCAATTTCTTAACAGAAATGtataagtcatacttttgtttt                                                                                                                               <  2:123775‑M1/149‑22 (MQ=255)                   cTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgg                                                                                                                             <  1:58037‑M1/149‑24 (MQ=255)                                    tCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtttttggtgtatttaaacttttt                                                                                                            >  2:384305‑M1/1‑109 (MQ=255)                                             cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                   >  2:290491‑M1/1‑100 (MQ=255)                                             cTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttt                                                                                                   >  2:659145‑M1/1‑100 (MQ=255)                                              tGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttt                                                                                                  >  1:591693‑M1/1‑99 (MQ=255)                                                    cGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgcta                                                                                            <  1:226742‑M1/149‑57 (MQ=255)                                                     gTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcattctgttaaaaagtttttcgctac                                                                                           >  1:400868‑M1/1‑92 (MQ=255)                                                           acacGCAGATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagc                                                                                     <  2:668676‑M1/149‑64 (MQ=255)                                                              cGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacatttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaag                                                                                  <  1:481427‑M1/149‑67 (MQ=255)                                                                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagccaaaatgaa                                                                        >  1:649588‑M1/1‑73 (MQ=255)                                                                        aCGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaa                                                                        >  2:392474‑M1/1‑73 (MQ=255)                                                                                   cgtcgatTTTTAATAATGCGGGAGGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtttatttcaatctgttaaaaagtttttcgctacgcgagcaagcaaaaatgaaacaggaataat                                                             <  2:47989‑M1/146‑88 (MQ=255)                                                                                   cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                             >  1:519986‑M1/1‑62 (MQ=255)                                                                                   cTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataat                                                             <  2:69505‑M1/149‑88 (MQ=255)                                                                                   cTGCGATTTTCAACAATGCGTGATGCAGAACACACAAAACACTCAATTACGTAACATAAATGtataagtcatacttttgttttgggtgtatgtcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataag                                                             <  1:679398‑M1/149‑88 (MQ=255)                                                                                    tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttttttttgttttatttcaatctgttaaaaagtttttcgttacgctagtaagcaaacgtgaaacaggagtaata                                                            >  2:93740‑M1/1‑61 (MQ=255)                                                                                       gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                     >  2:279785‑M1/1‑58 (MQ=255)                                                                                       gATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtattt                                                                                                                     <  1:279785‑M1/89‑32 (MQ=255)                                                                                        aTTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcatacaaaaatgaaacaggaataatctcaa                                                        >  2:85510‑M1/1‑57 (MQ=255)                                                                                         ttttCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgtgttgggtgtatttcaatctgttaaaaagtttttggctacgctagcaagcaaaaatgaaacaggaataatcgaaat                                                       >  2:384979‑M1/1‑56 (MQ=255)                                                                                            tCAATAAGTCGTGATGAAGATCATACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                    <  1:383266‑M1/149‑97 (MQ=255)                                                                                               ataatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                                 >  2:241925‑M1/1‑50 (MQ=255)                                                                                               ataatGACTGATCCAGAACAAACAAAACCCTCAATTACTTAACACAAAGGtataagtcctacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctccgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg                                                 <  1:42322‑M1/149‑100 (MQ=255)                                                                                                    gCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcg                                            <  1:275472‑M1/149‑105 (MQ=255)                                                                                                              gATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaa                                  <  1:91890‑M1/149‑115 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 >  1:712875‑M1/1‑34 (MQ=255)                                                                                                               aTCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaa                                 <  2:712875‑M1/149‑116 (MQ=255)                                                                                                                  acaaaAACCACGCAATTCCTTACCATAAATGtataagtcatacttttgtttggggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataa                              <  2:505603‑M1/149‑119 (MQ=255)                                                                                                                    acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataact                            >  1:703599‑M1/1‑29 (MQ=255)                                                                                                                     caAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactc                           >  1:471023‑M1/1‑28 (MQ=255)                                                                                                                             tCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaa                   >  1:307585‑M1/1‑20 (MQ=255)                                                                                                                                  tACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatc              <  2:588186‑M1/149‑135 (MQ=255)                                                                                                                                              aTGtataaggtagacttttgttttgcgtgtagttcaatcggttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaatactcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaa  <  2:620897‑M1/149‑147 (MQ=255)                                                                                                                                                |                                                                                                                                                   ATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAA  >  NC_000913_3_bme_pgi/1978359‑1978648

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