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–1979269	1979278	10–776	75 [0]	[1] 74	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)	65 (1.100)	50/182	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

GAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAT  >  NC_000913/1979181‑1979376                                                                                                   |                                                                                                  gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                                                  >  1:1104828‑M2/99‑100 (MQ=255)    ctcgtcactacacgcacatacaacggaggggggctgctattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                               >  2:844723‑M2/96‑100 (MQ=255)      cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                               <  2:338040‑M2/5‑1 (MQ=255)      cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                             <  2:345883‑M2/7‑1 (MQ=255)      cgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGt                                                                                             >  2:711503‑M2/94‑100 (MQ=255)           ctacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatgacacaaaacactcaattacttaacataaatgTATAAGTCATAc                                                                                        >  2:1007950‑M2/89‑100 (MQ=255)            tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                       >  2:148733‑M2/88‑100 (MQ=255)            tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                       <  1:1307217‑M2/13‑1 (MQ=255)              cacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttt                                                                                     <  2:1058878‑M2/15‑1 (MQ=255)               acgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACtttt                                                                                    <  1:564872‑M2/16‑1 (MQ=255)                cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  2:348585‑M2/17‑1 (MQ=255)                cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  1:81751‑M2/17‑1 (MQ=255)                 gcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  2:348750‑M2/17‑1 (MQ=255)                  cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 <  2:4153‑M2/19‑1 (MQ=255)                   acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 >  2:153625‑M2/81‑99 (MQ=255)                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                              <  2:1113936‑M2/22‑1 (MQ=255)                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                             >  2:813644‑M2/78‑100 (MQ=255)                       acaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                            <  2:933006‑M2/24‑1 (MQ=255)                        caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                           >  2:348402‑M2/76‑100 (MQ=255)                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                          >  1:1057794‑M2/75‑100 (MQ=255)                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                          >  1:1454640‑M2/75‑100 (MQ=255)                                        cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                           <  2:652524‑M2/41‑1 (MQ=255)                                        cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                           <  2:689390‑M2/41‑1 (MQ=255)                                             ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGTGTGTATTTCAATCTGTTAAAAAg                                                      >  2:151489‑M2/55‑100 (MQ=255)                                             ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                      <  1:882251‑M2/46‑1 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     <  1:657107‑M2/47‑1 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     >  1:635952‑M2/54‑100 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     <  1:717370‑M2/47‑1 (MQ=255)                                              tcaataatgcgtgatgcagatcaaacaaaacactcaattacttaacataaatgTATAAGTCATACTCTTGTTATGGGTGTCTCTCAATCTGTTAAACAg                                                      <  2:1319920‑M2/46‑1 (MQ=255)                                               caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                    >  2:874261‑M2/53‑100 (MQ=255)                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                   <  1:556397‑M2/49‑1 (MQ=255)                                                 ccgatctgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGcttt                                                  <  1:326779‑M2/50‑1 (MQ=255)                                                   aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                <  1:1224023‑M2/52‑1 (MQ=255)                                                   aatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTc                                                >  2:40978‑M2/49‑100 (MQ=255)                                                    atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                               <  2:445361‑M2/53‑1 (MQ=255)                                                    atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                               <  1:348402‑M2/53‑1 (MQ=255)                                                       cgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                            >  1:383025‑M2/45‑100 (MQ=255)                                                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTtaaaaagctttagatcgg                                           >  2:326779‑M2/44‑93 (MQ=255)                                                          gatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgct                                             >  2:1233755‑M2/42‑96 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCtgttaaaaagttttttttttttttt                                       >  2:1515671‑M2/40‑90 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  1:903087‑M2/40‑100 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       <  1:844728‑M2/61‑1 (MQ=255)                                                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                      >  2:1131330‑M2/39‑100 (MQ=255)                                                              cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagc                                     >  2:1044429‑M2/38‑100 (MQ=255)                                                               agatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagca                                    >  2:1101821‑M2/37‑100 (MQ=255)                                                                gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                   >  2:989217‑M2/36‑100 (MQ=255)                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                  <  1:190448‑M2/66‑1 (MQ=255)                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                  >  1:1381439‑M2/35‑100 (MQ=255)                                                                 atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                  <  1:131066‑M2/66‑1 (MQ=255)                                                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                 <  1:148733‑M2/67‑1 (MQ=255)                                                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                 <  2:1125030‑M2/67‑1 (MQ=255)                                                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagc                                 <  1:604044‑M2/67‑1 (MQ=255)                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                             >  2:249051‑M2/30‑100 (MQ=255)                                                                       caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                            >  2:368661‑M2/29‑100 (MQ=255)                                                                           acactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaa                        <  1:1333742‑M2/76‑1 (MQ=255)                                                                            cactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGaaa                       <  2:320251‑M2/77‑1 (MQ=255)                                                                                caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   >  1:611464‑M2/20‑100 (MQ=255)                                                                                caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   >  1:1082580‑M2/20‑100 (MQ=255)                                                                                  attacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                 >  1:353205‑M2/18‑100 (MQ=255)                                                                                    tacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaata               <  2:203861‑M2/85‑1 (MQ=255)                                                                                     acttaacataaatgTATAAGTAATACTTTTGTTTTGGGTGTATTTCAATATGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa              >  2:612194‑M2/15‑100 (MQ=255)                                                                                      cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTGAATTTGTTAAAATGTTTTTCGCTATGCTAGTAAGCAAAAATGAAACAGGgataat             >  2:1108285‑M2/14‑100 (MQ=255)                                                                                        taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg           >  2:1360491‑M2/12‑100 (MQ=255)                                                                                         aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa          >  2:372219‑M2/11‑100 (MQ=255)                                                                                         aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa          >  2:1282603‑M2/11‑100 (MQ=255)                                                                                            ataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt       >  2:1343112‑M2/8‑100 (MQ=255)                                                                                             taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg      <  2:729060‑M2/94‑1 (MQ=255)                                                                                             taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg      <  2:611462‑M2/94‑1 (MQ=255)                                                                                             taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg      >  1:114865‑M2/7‑100 (MQ=255)                                                                                              aaaTGTATAAGTCATACTTTTGTTTTGAGAATATTTCAATCTGTTAAAAAGTTTTTCGCTAAAAAAAAAAAAAAAAATGAAACAGGAATAATCGAAATg      <  2:941135/99‑1 (MQ=255)                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa   <  1:217787‑M2/97‑1 (MQ=255)                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa   <  1:379395‑M2/97‑1 (MQ=255)                                                                                                atgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa   <  1:967894‑M2/97‑1 (MQ=255)                                                                                                 tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt  <  2:1310075‑M2/98‑1 (MQ=255)                                                                                                   |                                                                                                  GAACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAT  >  NC_000913/1979181‑1979376

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑