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	26 [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)	24 (0.730)	22/584	0.4	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

ACGCCCATAATGCGGGCTGTTGCCCGGCATCCAACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACA  >  NC_000913/1978992‑1979578                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            ccaactttatttttatgcggtttcaccgcaccccgtgatgtcgcagggaagccccgggaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatggtcttataaatagaaatgggtctttacacttagctaagatttttcctaaatcgacgaaactgtacgcgtcactacacgcaaaacaaacggaggagggctgcaattttcactaatgagtgatgcaaaacacaaaaaacactcaaatacttaaaataaattATTAATTCAtat                                                                                                                                                                                                                                                                                                  >  1:231294‑M2/288‑298 (MQ=255)          ttatgattctgtttcaccgccccccgtgatgtcgcctggaagccccggtaaaaaataattagcaatagaatagttgcgataagctacaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttgtctaatatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatacgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                                                                                                                                                                                       <  2:231302‑M2/23‑1 (MQ=255)                     tttcaccgcaccccgtgatgtcgccgggaagccccggtaaaaaataatttgcattagaatagttgcgataagctgtaataagcagaatcaccttgttggtttcatatgaccttacaaaaagaaatgggtctttacacttatctaagattttccttaaatcgacgcaactgtattcgtcactacatgcacattcaatggaggggggctgtgcttttcattaatgagtgatgtagattacacaaaacacacaattacttaacataaatgTATAAGTAATACTTTTGTTTTGGGTGTATTTCaa                                                                                                                                                                                                                                                                            >  2:25779‑M2/268‑301 (MQ=255)                          ccgcaccccgtgatgtcgccgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAAt                                                                                                                                                                                                                                                                           >  2:35456‑M2/263‑297 (MQ=255)                             caccccgtgatgtcgccgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgttcttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactatacgaacatacaacggaggggggctgcgatcttcaataatgcgtgatgcagacaacacaaaacactcaattacttaacataaatgTATAATTCATACTTTTTTTTTGGGTGTATTTCAATCTATtat                                                                                                                                                                                                                                                                    >  1:239675‑M2/260‑300 (MQ=255)                                cccgtgatgtcgccgggaagccccggtaaaaaataattagcattagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                                                                                                                                                                                 <  2:215161‑M2/45‑1 (MQ=255)                                                                          ttagaatagttgcgataagctgcaataagcagaaccacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTAGCTACGCTAGCAAGCAAAAATGAAACAGgaataaa                                                                                                                                                                                                                       >  2:284823‑M2/215‑300 (MQ=255)                                                                          ttaaaatagttgctacgaatgtcaataagcagaatcacctttttgtgttaataattcctaacaaatcgaaatgggtctttacacttatctaaaattttgcctaaatcgacgcaactgtactcgtcactacacgcacatacaactgaggggggctgcgagtttcaataatgcgtgatgcagaacacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGGTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                                                                                                                                                                                                                       <  1:201915‑M2/87‑1 (MQ=255)                                                                                                  ataagcagaaacacctttttggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTc                                                                                                                                                                                               <  2:158604‑M2/111‑1 (MQ=255)                                                                                                                      ggtttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATAGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAaata                                                                                                                                                                           >  2:274270‑M2/171‑301 (MQ=255)                                                                                                                         ttaatatgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGTAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCAACGTaaataata                                                                                                                                                                        >  2:314974‑M2/168‑300 (MQ=255)                                                                                                                         ttaaaatgtctatactaatagacatgggtctttacaattatatagaatttttcctaaagcgaagcatctatagtcgtcactgcaaacacagacaacggatgagggctgcgattttcaatgatgcgtgatgaagatagcacaaaacactcaattacttaacatgaatgTATAAGTCCTACTTTTGTTGTGGATGTGTTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAAGAATAATCGAAATGGGATGTTGCGCACAATCAAAATAACTCACCGTAAATAatc                                                                                                                                                                        <  1:25778‑M2/134‑1 (MQ=255)                                                                                                                         cgatcttgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacactcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc                                                                                                                                                                        <  1:75295‑M2/134‑1 (MQ=255)                                                                                                                               tgtccttacaaatagaaatgggtctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAAAGGAATAATCGAAATGGGATGTTGCGCACAGGAAAAATAACTCACCGTaaataatcagatcg                                                                                                                                                                  >  2:75297‑M2/162‑296 (MQ=255)                                                                                                                                        tccctccacgacgctcttccgatcttctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTataaataa                                                                                                                                                         <  1:9417‑M2/149‑1 (MQ=255)                                                                                                                                                             cttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGa                                                                                                                                    >  2:228868‑M2/132‑301 (MQ=255)                                                                                                                                                                 tctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAactcaccgtaaataatcatctgctataaataaagatcggaagagcacacgtctgaa                                                                                                                                 >  2:9418‑M2/128‑276 (MQ=255)                                                                                                                                                                                    ccacgccgctttactcttacatacaagcaaagacaacggagtgggaatgcgattttcaataaggcagaattcagatcacacaaaacacttgataacttaacataaatgTATAAGTCATACTTATGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATACTCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACt                                                                                                             <  1:108109‑M2/193‑1 (MQ=255)                                                                                                                                                                                          aactgtaatcggtactacacgaacatacaacgaaggggggcttcaattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaaagTATAAGTCATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGc                                                                                                       <  1:274262‑M2/199‑1 (MQ=255)                                                                                                                                                                                                   cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttgacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAAc                                                                                              <  2:239683‑M2/208‑1 (MQ=255)                                                                                                                                                                                                                    acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaagaa                                                                             >  2:229680‑M2/77‑301 (MQ=255)                                                                                                                                                                                                                     caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACaaagaaa                                                                            >  2:261737‑M2/76‑301 (MQ=255)                                                                                                                                                                                                                                   tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTTGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCGAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACt                                                              <  1:35455‑M2/240‑1 (MQ=255)                                                                                                                                                                                                                                                    cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGAAACAATCAACTTAATGTTAATGAAATGt                                             >  2:148536‑M2/45‑301 (MQ=255)                                                                                                                                                                                                                                                                                         ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATc        <  1:228860‑M2/294‑1 (MQ=255)                                                                                                                                                                                                                                                                                              tttataAGTTATACTGTTGGTTTGGTTGTATTTCAAACGGTTAAAAAGTTTGTCGCTACGCGAGCAAGCAAAAATGAAACAGGAATTATCGAAATGGGATGTTGCTCACAGTCAGAATAACTCACCGTAAATAAGCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTac   <  1:284815/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGTAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAAATAATTAAATACAAAGAAAGGAAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAgtaaa  >  2:231180‑M2/2‑299 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                            ACGCCCATAATGCGGGCTGTTGCCCGGCATCCAACGCCATTCATGGCCATATCAATGATTTTCTGGTGCGTACCGGGTTGAGAAGCGGTGTAAGTGAACTGCAGTTGCCATGTTTTACGGCAGTGAGAGCAGAGATAGCGCTGATGTCCGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTATAAATAATCACTTTCATGCAATACCAGATAAGCTATTTTTAAACAGACACTTACCGCACAACAAACTAATTAAATACAAAGAAAGGCAACAATCAACTTAATGTTAATGAAATGTTTTACATATATTAACCACTTGTTTAATGCGGAGAATCAGTACA  >  NC_000913/1978992‑1979578

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