Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,121,521	Δ9 bp	coding (219‑227/1332 nt)	hslU ←	molecular chaperone and ATPase component of HslUV protease

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4121521	4121529	9	59 [0]	[0] 60	hslU	molecular chaperone and ATPase component of HslUV protease

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4121520	0 (0.000)	57 (0.670)	48/196	1.0	100%	coding (228/1332 nt)	hslU	molecular chaperone and ATPase component of HslUV protease
	?	NC_000913	4121530 =	0 (0.000)					coding (218/1332 nt)	hslU	molecular chaperone and ATPase component of HslUV protease

GCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGA  >  NC_000913/4121431‑4121629                                                                                                    |                                                                                                    ttcacggcggcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggca                                                                                                     <  1:1650755‑M2/1‑1 (MQ=255)   cacggcggcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGc                                                                                                   >  1:1943590‑M2/98‑100 (MQ=255)     cggcggcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCtt                                                                                                 <  1:1140726‑M2/5‑1 (MQ=255)     cggcggcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCtt                                                                                                 >  2:278239‑M2/96‑100 (MQ=255)          ncatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCa                                                                                            >  1:618316‑M2/91‑100 (MQ=255)          gcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCa                                                                                            <  2:1559137‑M2/10‑1 (MQ=255)          gcatcggtcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCa                                                                                            <  2:1859383‑M2/10‑1 (MQ=255)                 tcagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGacgacg                                                                                     >  1:769333‑M2/84‑100 (MQ=255)                  cagatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACgg                                                                                    <  1:368284‑M2/18‑1 (MQ=255)                   agatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACggg                                                                                   <  2:403482‑M2/19‑1 (MQ=255)                   agatcgcgaataatagaatccacttccttaccgacgtagcccactgcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACggg                                                                                   <  2:1074084‑M2/19‑1 (MQ=255)                    gatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGc                                                                                  >  2:1871750‑M2/81‑100 (MQ=255)                    gatcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGc                                                                                  >  1:846319‑M2/81‑100 (MQ=255)                     atcgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCg                                                                                 >  1:566513‑M2/80‑100 (MQ=255)                       cgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGAt                                                                               >  1:1785673‑M2/78‑100 (MQ=255)                       cgcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGAt                                                                               >  1:2143216‑M2/78‑100 (MQ=255)                        gcgaataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGAtt                                                                              >  1:2027721‑M2/77‑100 (MQ=255)                            ataatagaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAg                                                                          <  1:1388035‑M2/28‑1 (MQ=255)                                  gaatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTAc                                                                    >  2:673770‑M2/67‑100 (MQ=255)                                   aatccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTAcc                                                                   <  2:1329546‑M2/35‑1 (MQ=255)                                     tccacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCg                                                                  <  1:893923‑M2/36‑1 (MQ=255)                                      ccacttccttaccgacgtagcccacttcggtgaatttgatcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGac                                                                >  2:1419939‑M2/63‑100 (MQ=255)                                       cacttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGaca                                                               >  2:620639‑M2/62‑100 (MQ=255)                                        acttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGacac                                                              >  1:1370010‑M2/61‑100 (MQ=255)                                         cttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACAcc                                                             >  1:2138172‑M2/60‑100 (MQ=255)                                          ttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCg                                                            <  1:2003267‑M2/42‑1 (MQ=255)                                          ttccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCg                                                            <  1:848630‑M2/42‑1 (MQ=255)                                           tccttaccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCgg                                                           >  1:670555‑M2/58‑100 (MQ=255)                                                 ccgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGc                                                     >  2:929272‑M2/52‑100 (MQ=255)                                                  cgacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGcc                                                    <  1:929274‑M2/50‑1 (MQ=255)                                                   gacgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCg                                                   >  1:511733‑M2/50‑100 (MQ=255)                                                     cgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGat                                                 <  2:618316‑M2/53‑1 (MQ=255)                                                     cgtagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGa                                                  <  2:602894‑M2/52‑1 (MQ=255)                                                     atctgcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGat                                                 <  1:622009‑M2/53‑1 (MQ=255)                                                       tagcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGatca                                               >  2:1527657‑M2/46‑100 (MQ=255)                                                         gcccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGgccgatagat                                             >  2:622009‑M2/44‑96 (MQ=255)                                                          cccacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGAtcatca                                            >  1:1779744‑M2/43‑100 (MQ=255)                                                            cacttcggtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAgg                                          <  1:1871757‑M2/60‑1 (MQ=255)                                                                   gtgaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                   >  1:184422‑M2/34‑100 (MQ=255)                                                                     gaatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTc                                 >  2:1334310‑M2/32‑100 (MQ=255)                                                                      aatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCg                                <  1:1713775‑M2/70‑1 (MQ=255)                                                                      aatttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCg                                >  2:1538581‑M2/31‑100 (MQ=255)                                                                       atttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCgg                               <  2:1370007‑M2/71‑1 (MQ=255)                                                                       atttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGACGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCgg                               <  2:1943582‑M2/71‑1 (MQ=255)                                                                        tttggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCggg                              <  1:160431‑M2/72‑1 (MQ=255)                                                                           ggtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTc                           >  2:1682998‑M2/26‑100 (MQ=255)                                                                            gtcgcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCa                          >  1:1974661‑M2/25‑100 (MQ=255)                                                                               gcttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACtt                       >  1:892218‑M2/22‑100 (MQ=255)                                                                                cttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTc                      >  2:1852422‑M2/21‑100 (MQ=255)                                                                                 ttcaactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCa                     <  1:1419942‑M2/81‑1 (MQ=255)                                                                                   caactttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATg                   >  1:653711‑M2/18‑100 (MQ=255)                                                                                     actttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGgc                 >  2:553915‑M2/16‑100 (MQ=255)                                                                                       tttgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGgcgc               <  2:670554‑M2/87‑1 (MQ=255)                                                                                         tgatgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAg             <  2:566514‑M2/89‑1 (MQ=255)                                                                                           atgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGct           >  2:1925187‑M2/10‑100 (MQ=255)                                                                                           atgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGct           >  2:1435700‑M2/10‑100 (MQ=255)                                                                                            tgaacggcAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGctc          >  1:1201326‑M2/9‑100 (MQ=255)                                                                                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg   <  1:334877‑M2/99‑1 (MQ=255)                                                                                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg   <  1:620639‑M2/99‑1 (MQ=255)                                                                                                    aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa  <  1:673771/100‑1 (MQ=255)                                                                                                    |                                                                                                    GCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGA  >  NC_000913/4121431‑4121629

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