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	84 [0]	[0] 83	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)	77 (1.120)	60/182	0.1	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

AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979182‑1979378                                                                                                  |                                                                                                    tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtat                                                                                                   >  2:1744931‑M2/98‑100 (MQ=255) tactcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgta                                                                                                    >  2:1604359‑M2/98‑99 (MQ=255)  actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                  >  1:1385137‑M2/97‑100 (MQ=255)  actcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgtata                                                                                                  >  2:178587‑M2/97‑100 (MQ=255)   ctcgtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATaa                                                                                                 >  2:1043842‑M2/96‑100 (MQ=255)      gtcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTc                                                                                              <  2:550741‑M2/8‑1 (MQ=255)       tcactacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCa                                                                                             <  1:169152‑M2/9‑1 (MQ=255)         actacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATa                                                                                           >  2:1558791‑M2/90‑100 (MQ=255)           tacacgcacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACt                                                                                         <  2:1424392‑M2/13‑1 (MQ=255)                 cacatacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtt                                                                                   >  2:1071031‑M2/82‑100 (MQ=255)                   catacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgtttt                                                                                 <  1:1431033‑M2/21‑1 (MQ=255)                    atacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACttttgttttg                                                                                >  2:249335‑M2/79‑100 (MQ=255)                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  1:1592752‑M2/78‑100 (MQ=255)                     tacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  1:1127889‑M2/78‑100 (MQ=255)                     nacaacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTgg                                                                               >  2:1575105‑M2/78‑100 (MQ=255)                      acaactgagggggtctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTggg                                                                              <  2:1418400‑M2/24‑1 (MQ=255)                       catcgcaggggggctgccattttcaataatgcgtgatgcagatcacacaaaacactccattacttaacataaatgTATAAGTCATACttatcttttg                                                                                >  1:906423‑M2/76‑97 (MQ=255)                       caacggaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGgt                                                                             <  2:836323‑M2/25‑1 (MQ=255)                            gaggggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTAtt                                                                        <  1:1478384‑M2/30‑1 (MQ=255)                              ggggggctgcgattgtcaataatgcgtgatgcagatcacacaaaacactcaatgtcttgacataaatgTATAAGTCATACTTTTGTTGTGGGTGTAGTTc                                                                      <  1:1580444‑M2/32‑1 (MQ=255)                               gggggctgcgattatcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                     >  2:1232068‑M2/68‑100 (MQ=255)                                ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                    >  1:1701399‑M2/67‑100 (MQ=255)                                ggggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCaa                                                                    <  1:1325068‑M2/34‑1 (MQ=255)                                ggggctgcgaatttcaataatgcgtgatgcaaatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCa                                                                     <  1:959465‑M2/33‑1 (MQ=255)                                  ggctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGACATACTTTTGTTTTGGGTGTATTTCAATc                                                                  >  2:613047‑M2/65‑100 (MQ=255)                                   gctgcgatttttaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATc                                                                  <  2:611443‑M2/36‑1 (MQ=255)                                   gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                 <  1:1748795‑M2/37‑1 (MQ=255)                                   gctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCt                                                                 >  1:225926‑M2/64‑100 (MQ=255)                                    ctgcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTg                                                                >  2:661846‑M2/63‑100 (MQ=255)                                      gcgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              >  1:1581696‑M2/61‑100 (MQ=255)                                       cgattttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGtt                                                              <  1:1593341‑M2/40‑1 (MQ=255)                                           tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                         <  2:622274‑M2/45‑1 (MQ=255)                                           tttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTaaaaa                                                         >  1:665611‑M2/56‑100 (MQ=255)                                            ttcaataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAg                                                        <  2:1152744‑M2/46‑1 (MQ=255)                                              caataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtt                                                      >  1:590490‑M2/53‑100 (MQ=255)                                                ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                    >  1:1775792‑M2/51‑100 (MQ=255)                                                ataatgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGtttt                                                    >  1:1479118‑M2/51‑100 (MQ=255)                                                  aatgcgtgctgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTCGGCGTATTTCAATCTGTTAAACATtttt                                                    >  1:1785430‑M2/49‑98 (MQ=255)                                                   atgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcg                                                 <  2:1316169‑M2/53‑1 (MQ=255)                                                    tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                >  1:208612‑M2/47‑100 (MQ=255)                                                    tgcgtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgc                                                <  1:6035‑M2/54‑1 (MQ=255)                                                       gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctac                                             >  2:273566‑M2/44‑100 (MQ=255)                                                       gtgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgcta                                              >  2:889703‑M2/44‑99 (MQ=255)                                                        tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                            >  1:1682246‑M2/43‑100 (MQ=255)                                                        tgatgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacg                                            >  2:95509‑M2/43‑100 (MQ=255)                                                          atgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgct                                          <  1:186387‑M2/60‑1 (MQ=255)                                                           tgcagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTcgctacgcta                                         <  1:628964‑M2/61‑1 (MQ=255)                                                             cagatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTag                                        <  2:896697‑M2/62‑1 (MQ=255)                                                               gatcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaa                                     >  1:488153‑M2/36‑100 (MQ=255)                                                                atcacacaaaatactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                    <  2:308702‑M2/66‑1 (MQ=255)                                                                atcacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGTTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaag                                    <  1:395359‑M2/66‑1 (MQ=255)                                                                  cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                  <  1:31340‑M2/68‑1 (MQ=255)                                                                  cacacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTagcaagca                                  >  1:1060831‑M2/33‑100 (MQ=255)                                                                   acacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTCTTTCAATCTCTTAAAAAGTTTTTCGCTACGCTAgcaagcaa                                 >  2:410607‑M2/32‑100 (MQ=255)                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                <  2:1127900‑M2/70‑1 (MQ=255)                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  2:1125161‑M2/31‑100 (MQ=255)                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaa                                >  1:995943‑M2/31‑100 (MQ=255)                                                                    cacaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGATACGCTAGCAAGCaaa                                >  1:1064743‑M2/31‑100 (MQ=255)                                                                     acaaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCaaaa                               <  2:148091‑M2/71‑1 (MQ=255)                                                                      caaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAATCaaaaa                              >  1:627869‑M2/29‑100 (MQ=255)                                                                       aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  1:1686882‑M2/73‑1 (MQ=255)                                                                       aaaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAAt                             <  2:1069545‑M2/73‑1 (MQ=255)                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            <  2:1592763‑M2/74‑1 (MQ=255)                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            >  2:979943‑M2/27‑100 (MQ=255)                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            >  2:808270‑M2/27‑100 (MQ=255)                                                                        aaacactcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATg                            >  1:325264‑M2/27‑100 (MQ=255)                                                                            actcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAAc                        >  2:431803‑M2/23‑100 (MQ=255)                                                                              tcaattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:343688‑M2/21‑100 (MQ=255)                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     >  2:1675714‑M2/20‑100 (MQ=255)                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     >  1:224437‑M2/20‑100 (MQ=255)                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  2:1338456‑M2/81‑1 (MQ=255)                                                                               caattacttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  2:695117‑M2/20‑99 (MQ=255)                                                                                     cttaacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  2:831334‑M2/14‑100 (MQ=255)                                                                                       taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  1:252284‑M2/12‑98 (MQ=255)                                                                                       taacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg             <  1:1405380‑M2/89‑1 (MQ=255)                                                                                       taacataaatgTATAAGTCATACTTTGGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg             >  2:1102939‑M2/12‑100 (MQ=255)                                                                                        aacataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa            >  2:658824‑M2/11‑100 (MQ=255)                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:1370810‑M2/9‑100 (MQ=255)                                                                                          cataaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGAAATAATCGaaa          >  2:969968‑M2/9‑100 (MQ=255)                                                                                            taaatgTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg        <  1:339950‑M2/94‑1 (MQ=255)                                                                                                 gTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   >  2:314913‑M2/2‑100 (MQ=255)                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGGTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  >  1:772206/1‑100 (MQ=255)                                                                                                  tataAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  2:225927/100‑1 (MQ=255)                                                                                                  |                                                                                                    AACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCCATAAATGTATAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979182‑1979378

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