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	71 [0]	[1] 72	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)	66 (1.030)	50/182	0.3	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:684247‑M2/2‑1 (MQ=255) gtactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                                                  <  1:684244‑M2/2‑1 (MQ=255)  tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                 <  1:1585409‑M2/3‑1 (MQ=255)  tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                 <  1:1567283‑M2/3‑1 (MQ=255)       gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                            <  2:484439‑M2/8‑1 (MQ=255)          actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                         >  1:1612046‑M2/90‑100 (MQ=255)                cgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttg                                                                                   <  1:1044522‑M2/17‑1 (MQ=255)                  ttccgatctcggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 <  1:1557001‑M2/19‑1 (MQ=255)                  cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 >  2:770736‑M2/82‑100 (MQ=255)                  cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                 <  2:130762‑M2/19‑1 (MQ=255)                   acatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttt                                                                                <  2:967098‑M2/20‑1 (MQ=255)                     atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                              >  2:1633928‑M2/79‑100 (MQ=255)                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                             <  2:1567843‑M2/23‑1 (MQ=255)                      tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                             <  1:217645‑M2/23‑1 (MQ=255)                        caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                           >  1:537239‑M2/76‑100 (MQ=255)                        caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                           <  1:316208‑M2/25‑1 (MQ=255)                        caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                           >  1:1674993‑M2/76‑100 (MQ=255)                         aacgggggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                            <  2:558772‑M2/24‑1 (MQ=255)                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                          >  1:796917‑M2/75‑100 (MQ=255)                         aacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgtg                                                                          >  1:438817‑M2/75‑100 (MQ=255)                           cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaagtcatacttttgttagatcggaa                                                                        >  2:1557008‑M2/73‑91 (MQ=255)                           cggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTa                                                                        <  2:1101366‑M2/28‑1 (MQ=255)                             gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                      <  1:27824‑M2/30‑1 (MQ=255)                             gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                      >  2:322814‑M2/71‑100 (MQ=255)                                gggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTc                                                                    <  1:208579‑M2/32‑1 (MQ=255)                                 ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                  <  2:1467341‑M2/34‑1 (MQ=255)                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                               <  2:950153‑M2/37‑1 (MQ=255)                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                               <  2:1458163‑M2/37‑1 (MQ=255)                                    gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                               <  2:1458288‑M2/37‑1 (MQ=255)                                      tgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAatcag                                                              >  1:1495044‑M2/62‑97 (MQ=255)                                        cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTa                                                           >  1:404954‑M2/60‑100 (MQ=255)                                         gattttcaataatgcgtgatgcagatcacacaagacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                          >  1:1020441‑M2/59‑100 (MQ=255)                                         gattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaa                                                          <  2:388036‑M2/42‑1 (MQ=255)                                            tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                       <  1:1650394‑M2/45‑1 (MQ=255)                                            tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                       >  2:739316‑M2/56‑100 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     >  1:819310‑M2/54‑100 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     >  1:1274088‑M2/54‑100 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGt                                                     <  2:387530‑M2/47‑1 (MQ=255)                                              tcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTCGGTGTATTTCAATCTGtt                                                            <  1:1377381‑M2/40‑1 (MQ=255)                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                   >  2:1340855‑M2/52‑100 (MQ=255)                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                   >  2:1419291‑M2/52‑100 (MQ=255)                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                   >  1:1519486‑M2/52‑100 (MQ=255)                                                aataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGttt                                                   <  1:1336517‑M2/49‑1 (MQ=255)                                                        gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                           <  1:631406‑M2/57‑1 (MQ=255)                                                           atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                        <  2:514436‑M2/60‑1 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  1:210297‑M2/40‑100 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  1:1299921‑M2/40‑100 (MQ=255)                                                            tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                       >  2:203612‑M2/40‑100 (MQ=255)                                                             gcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                      <  2:621321‑M2/62‑1 (MQ=255)                                                                  tcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGCTAAAAAGTTTTTCGCTACGCTagcaagc                                 >  1:449169‑M2/34‑100 (MQ=255)                                                                   cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                >  2:326338‑M2/33‑100 (MQ=255)                                                                    acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTATTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                               >  1:73152‑M2/32‑100 (MQ=255)                                                                      acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                             <  2:1251935‑M2/71‑1 (MQ=255)                                                                       caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaaa                            >  2:248154‑M2/29‑100 (MQ=255)                                                                        aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                           <  2:210297‑M2/73‑1 (MQ=255)                                                                         aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTAAATCTGTTAAAAAGTTTTTCGCTACGCTAGCATGCAAAAAt                           >  1:161741‑M2/27‑99 (MQ=255)                                                                          aacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGa                         <  1:524825‑M2/75‑1 (MQ=255)                                                                               ttaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                    <  1:248154‑M2/80‑1 (MQ=255)                                                                                caattacttaacataaatgTATAAGTTATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                   <  1:267257‑M2/81‑1 (MQ=255)                                                                                caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                    >  1:1385001‑M2/20‑99 (MQ=255)                                                                                 aattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGa                  >  1:675010‑M2/19‑100 (MQ=255)                                                                                     acttatcataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTATAAAGTTTTTCGCTACGCTTGCAAGCAAAGATGAAACAGGaat                >  1:674513‑M2/15‑98 (MQ=255)                                                                                     acttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa              >  2:795135‑M2/15‑100 (MQ=255)                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc            <  1:207515‑M2/88‑1 (MQ=255)                                                                                       ttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc            <  2:796537‑M2/88‑1 (MQ=255)                                                                                       ttaacATAAATGTATAAGTCATACTTTTGTTTTGTGTGTATAAAAAAAAAAAAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc            <  1:1142683/96‑1 (MQ=255)                                                                                        taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc            >  1:1064899‑M2/12‑99 (MQ=255)                                                                                           cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAatgaaacaggagatcggaaga        >  2:1134335‑M2/9‑90 (MQ=255)                                                                                           cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa        >  1:1269020‑M2/9‑100 (MQ=255)                                                                                             taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg      <  1:1402190‑M2/94‑1 (MQ=255)                                                                                               aatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg    <  2:1274094‑M2/96‑1 (MQ=255)                                                                                                 tgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt  >  2:348253‑M2/3‑100 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑