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–1979270	1979278	9–776	52 [0]	[0] 52	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)	48 (0.870)	37/268	0.5	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

CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979140‑1979421                                                                                                                                            |                                                                                                                                               tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                             >  1:616198‑M2/140‑143 (MQ=255) tctttacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                                                             <  2:65066‑M2/4‑1 (MQ=255)  ctgtacacttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                                                            <  1:442384‑M2/5‑1 (MQ=255)        gctgatccaagacttctgctaaagcgaggcaactgtactcgtcactacacgcacatacaacggaggggggctgcgcgttgcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      <  1:323417‑M2/11‑1 (MQ=255)        acttatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaccactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      >  2:551330‑M2/133‑143 (MQ=255)        acctatctaagatttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaactgaggtgggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                                                                      <  1:620204‑M2/11‑1 (MQ=255)                   agttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggagggagggagggatgtgcaagaatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                                                           <  2:467853‑M2/22‑1 (MQ=255)                    tttttcctacatcgacgccactgtactcgtcactacacgcacatacagcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:120884‑M2/23‑1 (MQ=255)                    tttttcctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgccgatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                                                                          <  2:331549‑M2/23‑1 (MQ=255)                          ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAt                                                                                                                    >  2:397627‑M2/115‑143 (MQ=255)                          ctaaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTGTGTTATGGGTGTAt                                                                                                                    >  1:731673‑M2/115‑143 (MQ=255)                           taaatcgacgcaactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaacacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                                                                   >  1:103875‑M2/114‑143 (MQ=255)                                      aactgtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                                                                        >  2:334658‑M2/103‑143 (MQ=255)                                          gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                                                                    <  2:564918‑M2/45‑1 (MQ=255)                                              tcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCAGACTATTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                                                                >  1:522489‑M2/95‑143 (MQ=255)                                                gtcactgcacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagaacacacattgccatcaatccctcaccataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttttt                                                                                              <  2:469653‑M2/51‑1 (MQ=255)                                                   actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTcgc                                                                                           >  2:14973‑M2/90‑143 (MQ=255)                                                     tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgata                                                                                         >  2:206319‑M2/88‑143 (MQ=255)                                                       cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                                                                       <  2:451070‑M2/58‑1 (MQ=255)                                                           cacatacaacgtaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                                                                   <  1:369113‑M2/62‑1 (MQ=255)                                                            acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                                                                  <  2:766755‑M2/63‑1 (MQ=255)                                                              atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTAGAAGTCAAACTTTTGTATTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                                                                >  1:761077‑M2/79‑143 (MQ=255)                                                                   acggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                                                           >  2:901309‑M2/74‑143 (MQ=255)                                                                     ggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                                                                         <  2:174461‑M2/72‑1 (MQ=255)                                                                       aggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                                                                       >  2:334445‑M2/70‑143 (MQ=255)                                                                            ggctgcgtgtgtcaataatgcgtgatgccgatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACa                                                                  <  1:632798‑M2/79‑1 (MQ=255)                                                                            ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAATTTTTCGCTACGCTAGCAAGCAAAAatcaacaa                                                                  >  1:594833‑M2/65‑140 (MQ=255)                                                                            ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTAGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTAAGCTAGCATGCAAAAATGAAACa                                                                  >  1:595247‑M2/65‑143 (MQ=255)                                                                             gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  1:222203‑M2/80‑1 (MQ=255)                                                                             gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGAGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                                                                 <  1:276697‑M2/80‑1 (MQ=255)                                                                              ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGATGTATTTCAATCTGGTAAAAAGTTCTTCGCCACgctagccagcaaaacggcaacaag                                                                >  1:336008‑M2/63‑133 (MQ=255)                                                                               tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                                                               <  1:274690‑M2/82‑1 (MQ=255)                                                                                 cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                                                             >  2:102544‑M2/60‑143 (MQ=255)                                                                                         aatactccgtgatgcagctcacacaaaacacacaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                     <  2:350949‑M2/92‑1 (MQ=255)                                                                                         aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa                                                     <  2:62163‑M2/92‑1 (MQ=255)                                                                                              tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa                                                <  2:616198‑M2/97‑1 (MQ=255)                                                                                               gcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATCCGTACTCTCCTAGTTTTGGGTGTATTTCTAACAGTTAAAAAGTTTTTCGCTACGTCAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt                                               >  2:470160‑M2/46‑143 (MQ=255)                                                                                                 gtgatgcagatcacacaaaaaactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt                                             >  1:811335‑M2/44‑143 (MQ=255)                                                                                                         gatcacagaaatcgcgcaattagtgaacataaatgTATAAGTCATACTGTTGTTTTGGGTGTATTTCAATCTGTGAAAAAGGTTTTCGCTACGCGAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTGGCGcaca                                     <  1:644213‑M2/108‑1 (MQ=255)                                                                                                             acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                 >  2:321708‑M2/32‑143 (MQ=255)                                                                                                             acacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTTTGGGTATATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCa                                 >  1:230139‑M2/32‑143 (MQ=255)                                                                                                               acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  2:30599‑M2/30‑143 (MQ=255)                                                                                                               acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  2:642268‑M2/30‑143 (MQ=255)                                                                                                               acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCaaa                               >  1:339343‑M2/30‑143 (MQ=255)                                                                                                               acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATAAAACAGGAATAATCGAAATAGTATGTTACGCACAGTcaac                               >  1:333671‑M2/30‑142 (MQ=255)                                                                                                                 aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATATCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAAACGAAATGGGAAGTTGACCACAATCAAAAt                             >  1:421362‑M2/28‑143 (MQ=255)                                                                                                                   aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATTTTGCGCACAGTAAAAATaa                           >  1:352458‑M2/26‑143 (MQ=255)                                                                                                                                   acataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatc           >  1:543162‑M2/10‑143 (MQ=255)                                                                                                                                     ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTGTTTACTACTCTCGAATGCCAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAatcat         >  1:214276‑M2/8‑143 (MQ=255)                                                                                                                                           gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTa   >  2:736413‑M2/2‑143 (MQ=255)                                                                                                                                           gTATAAGTAATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTa   <  1:743177‑M2/142‑1 (MQ=255)                                                                                                                                            tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCtat  <  2:118448/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                               CGGCGGTGCTTTTGCCGTTACGCACCACCCCGTCAGTAGCTGAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTTGCGCACAGTCAAAATAACTCACCGTAAATAATCATCTGCTAT  >  NC_000913/1979140‑1979421

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