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	1979278	776	27 [0]	[0] 27	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)	26 (0.810)	25/584	0.2	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

CCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGACAGCAGGCTCATAAGACGCCCCAGCGTCGCCATAGTGCGTTCACC  >  NC_000913/1978269‑1978796                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       ccGGGAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGTGTCTTTACACTTATCTAAGATTTTTCCCAAATCGACGCAACTGTACTCGTAACTACACGAACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagc                                                                                                                                                                                                                                     <  2:265555‑M1/301‑68 (MQ=255)      aaaaaaTAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACCCAATTACTTAACATAAATGtataagtcatacttttgttttcggtgtatttcaatctgttaaaaagttcttcgctacgctagcaagcaaaaa                                                                                                                                                                                                                                >  1:145693‑M1/1‑229 (MQ=255)               ttAGAATTAAAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTAGTTTAATATGTCCTTACAAATAGAAATGGGTCCTTACACTTATCTAAGAATTTTCCTAAATCGACGCTACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTAATGGAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacagg                                                                                                                                                                                                                       <  2:279175‑M1/301‑82 (MQ=255)                               gtgcataagCTGCAAAAAGCAAAACCACCTTTTTGGTTTAATATGTCCTTAAAAATATAAATGGGTCTTTACACTTATCTAAGATTGTTCCTAAATAGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg                                                                                                                                                                                                        <  2:259978‑M1/296‑97 (MQ=255)                                      aGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACAGTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctaccctagcaagcaaaaatgaaacaggaataatcaaaatggatagttgcg                                                                                                                                                                                                >  2:120754‑M1/1‑197 (MQ=255)                                                             ttttGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacattcaaaataactcaccgta                                                                                                                                                                         >  1:3509‑M1/1‑174 (MQ=255)                                                                tGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacgggaataatcgaaatgggatgttgcgcacagtaaaaataaacaaccgtaaaa                                                                                                                                                                      >  2:184975‑M1/1‑171 (MQ=255)                                                                tGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaaa                                                                                                                                                                      >  2:60888‑M1/1‑171 (MQ=255)                                                                                       aGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaaccgtaaataatcatctgctataaataatcaa                                                                                                                                               >  2:6752‑M1/1‑148 (MQ=255)                                                                                                                        ccTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgtgatgttgcgcacagtcaaaataaatcaacgtaaagaatcatctgctataaataatcactttcatgcaataccagaaaagctatttttaaaa                                                                                                              >  2:162005‑M1/1‑115 (MQ=255)                                                                                                                           aaaTCGACGCAACTGTACTCGTCACTACATGCACATACAACTGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataattatctgctataaataaacactttcatgcaataacagataagctatttttaaacaga                                                                                                           >  2:89916‑M1/1‑112 (MQ=255)                                                                                                                                                            caTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactaacagtaaataatcatctgctataaataatcactttcatgaaataccagaaaagctatttttaaacagacaattaccgcacaacaaactaattaaataaaaa                                                                          >  2:59030‑M1/1‑79 (MQ=255)                                                                                                                                                             aTACAACGGAGGGGGGCTGCGATTTTCAATAATGCGGGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacgggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaag                                                                         <  1:184985‑M1/301‑224 (MQ=255)                                                                                                                                                                cAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttccgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaa                                                                      >  1:158097‑M1/1‑75 (MQ=255)                                                                                                                                                                             ttgCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttgttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaac                                                         <  2:3509‑M1/300‑240 (MQ=255)                                                                                                                                                                              tGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaact                                                        <  1:60892‑M1/301‑241 (MQ=255)                                                                                                                                                                                           aatGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgctttgtgtgtatttcaatctgttaaaaagtttttcgccacgctagcaagcataaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaatactcactttcatgtaataccagatcagttatttttaaacagacacttaccgcaaaacaaactaattaaataccactacaggcaacaattaacttaatgttaatgaa                                           >  1:247778‑M1/1‑48 (MQ=255)                                                                                                                                                                                           aaatcgagatgcagaTAACAAAAAATAAGCAATTTCTTAACATAAATGtataaggcatacttttgttttggatgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaaggaaacaggaataatctaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatg                                             <  2:247767‑M1/292‑252 (MQ=255)                                                                                                                                                                                               cGTGATGCAGATCACACAAAACACTCAATTCCTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaataccaagaaaggcaacaatcaacttaatgttaatgaaatgt                                       >  1:287946‑M1/1‑44 (MQ=255)                                                                                                                                                                                                 tGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcatcgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgaacaaaaaactaattaaatacaaaaaaaggaaacaataaacttaatgttaataaaatattt                                     >  2:257887‑M1/1‑42 (MQ=255)                                                                                                                                                                                                   aTGCAGATCACACAAAACACTCAATTACATAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatgaaatgtttt                                    <  1:120763‑M1/300‑261 (MQ=255)                                                                                                                                                                                                    tGCAGATCACACAAAACACTCAATTACTTAACATACATGtataagtcatacttttgtgttgggtttatttcaatctgttaaaaagtttttcgctacgctagcaagcaaacatgaaacaggaataatcgcaatgggatgttgcgcacagtccaaagaactacccgtaaataatcatatgctataaatcatcaatttcatgcaatacaagagaatctctttttaaaaagaaccttcctgcacaacaaacttattaaatagatagaaatccaacaagcaacttaatgttaatcaaattcttaca                                  >  2:207412‑M1/1‑39 (MQ=255)                                                                                                                                                                                                            acacaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatagaaatgggatgttgcgcgcagtcaaaataactaaccgtaaataatcatctgctataaataatcactttcatgcaataccacaaaagatatttttaaacaagacattatctctcaaaaacataattaaataaaaagaaaagaaaacatcaaattaatgttaaagaaaagctttaaatatatat                          >  2:244677‑M1/1‑31 (MQ=255)                                                                                                                                                                                                              acaAAACACTCAATTACTTAACATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatgaaatgttttacatatattaaa                        >  1:26603‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                    cacTCAATTACTTAAAAAAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagttttccgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatgaaatgttttacatatattaaccacttg                  <  2:158087‑M1/301‑279 (MQ=255)                                                                                                                                                                                                                                  aCATAAATGtataagtcatacttttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatgaaatgttttacatatattaaccacttgtttaatgcggagaa    <  1:162015‑M1/301‑293 (MQ=255)                                                                                                                                                                                                                                    aTAAATGtataagtcataattttgttttgggtgtatttcaatctgttaaaaagtttttcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgttgcgcacagtcaaaataactcaccgtaaataatcatctgctataaataatcactttcatgcaataccagataagctatttttaaacagacacttaccgcacaacaaactaattaaatacaaagaaaggcaacaatcaacttaatgttaatgaaatgttttacatatattaaccacttgtttaatgcggagaatc  <  1:244688‑M1/301‑295 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                       CCGGTAAAAAATAATTAGCATTAGAATAGTTGCGATAAGCTGCAATAAGCAGAACCACCTTTTTGGTTTAATATGTCCTTACAAATAGAAATGGGTCTTTACACTTATCTAAGATTTTTCCTAAATCGACGCAACTGTACTCGTCACTACACGCACATACAACGGAGGGGGGCTGCGATTTTCAATAATGCGTGATGCAGATCACACAAAACACTCAATTACTTAACATAAATGGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCATACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTGACAGCAGGCTCATAAGACGCCCCAGCGTCGCCATAGTGCGTTCACC  >  NC_000913/1978269‑1978796

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑