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

TTCACGGCGGCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGC  >  NC_000913/4121422‑4121619                                                                                                   |                                                                                                    ttCACGGCGGCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCa                                                                                                    <  1:1650755‑M1/100‑2 (MQ=255)   cACGGCGGCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagc                                                                                                  >  1:1943590‑M1/1‑97 (MQ=255)     cggcggCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagctt                                                                                                <  1:1140726‑M1/100‑6 (MQ=255)     cggcggCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagctt                                                                                                >  2:278239‑M1/1‑95 (MQ=255)          ncaTCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagcca                                                                                           >  1:618316‑M1/2‑90 (MQ=255)          gcATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagcca                                                                                           <  2:1559137‑M1/100‑11 (MQ=255)          gcATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagcca                                                                                           <  2:1859383‑M1/100‑11 (MQ=255)                 tCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacg                                                                                    >  1:769333‑M1/1‑83 (MQ=255)                  cAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgg                                                                                   <  1:368284‑M1/100‑19 (MQ=255)                   aGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacggg                                                                                  <  2:403482‑M1/100‑20 (MQ=255)                   aGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTGCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacggg                                                                                  <  2:1074084‑M1/100‑20 (MQ=255)                    gATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggc                                                                                 >  1:846319‑M1/1‑80 (MQ=255)                    gATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggc                                                                                 >  2:1871750‑M1/1‑80 (MQ=255)                     aTCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcg                                                                                >  1:566513‑M1/1‑79 (MQ=255)                       cgcgAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgat                                                                              >  1:1785673‑M1/1‑77 (MQ=255)                       cgcgAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgat                                                                              >  1:2143216‑M1/1‑77 (MQ=255)                        gcgAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatt                                                                             >  1:2027721‑M1/1‑76 (MQ=255)                            ataataGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcag                                                                         <  1:1388035‑M1/100‑29 (MQ=255)                                  gAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttac                                                                   >  2:673770‑M1/1‑66 (MQ=255)                                   aaTCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttacc                                                                  <  2:1329546‑M1/100‑36 (MQ=255)                                     tCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccg                                                                 <  1:893923‑M1/99‑37 (MQ=255)                                      ccACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGATCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgac                                                               >  2:1419939‑M1/1‑62 (MQ=255)                                       cACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgaca                                                              >  2:620639‑M1/1‑61 (MQ=255)                                        aCTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacac                                                             >  1:1370010‑M1/1‑60 (MQ=255)                                         cTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacacc                                                            >  1:2138172‑M1/1‑59 (MQ=255)                                          ttCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccg                                                           <  1:2003267‑M1/100‑43 (MQ=255)                                          ttCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccg                                                           <  1:848630‑M1/100‑43 (MQ=255)                                           tCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccgg                                                          >  1:670555‑M1/1‑57 (MQ=255)                                                 ccGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggc                                                    >  2:929272‑M1/1‑51 (MQ=255)                                                  cGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggcc                                                   <  1:929274‑M1/100‑51 (MQ=255)                                                   gACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccg                                                  >  1:511733‑M1/1‑49 (MQ=255)                                                     cGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgat                                                <  2:618316‑M1/100‑54 (MQ=255)                                                     cGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccga                                                 <  2:602894‑M1/99‑53 (MQ=255)                                                     atctgcccaCTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgat                                                <  1:622009‑M1/96‑54 (MQ=255)                                                       tAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatca                                              >  2:1527657‑M1/1‑45 (MQ=255)                                                         gCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatagat                                            >  2:622009‑M1/1‑43 (MQ=255)                                                          cccACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatca                                           >  1:1779744‑M1/1‑42 (MQ=255)                                                            cACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcagg                                         <  1:1871757‑M1/100‑61 (MQ=255)                                                                   gTGAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttt                                  >  1:184422‑M1/1‑33 (MQ=255)                                                                     gAATTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttc                                >  2:1334310‑M1/1‑31 (MQ=255)                                                                      aaTTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcg                               >  2:1538581‑M1/1‑30 (MQ=255)                                                                      aaTTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcg                               <  1:1713775‑M1/100‑71 (MQ=255)                                                                       aTTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcgg                              <  2:1370007‑M1/100‑72 (MQ=255)                                                                       aTTTGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacggacgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcgg                              <  2:1943582‑M1/100‑72 (MQ=255)                                                                        tttGGTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggg                             <  1:160431‑M1/100‑73 (MQ=255)                                                                           ggTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtc                          >  2:1682998‑M1/1‑25 (MQ=255)                                                                            gTCGCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtca                         >  1:1974661‑M1/1‑24 (MQ=255)                                                                               gCTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcactt                      >  1:892218‑M1/1‑21 (MQ=255)                                                                                cTTCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttc                     >  2:1852422‑M1/1‑20 (MQ=255)                                                                                 ttCAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttca                    <  1:1419942‑M1/100‑82 (MQ=255)                                                                                   cAACTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatg                  >  1:653711‑M1/1‑17 (MQ=255)                                                                                     aCTTTGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggc                >  2:553915‑M1/1‑15 (MQ=255)                                                                                       tttGATGAACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgc              <  2:670554‑M1/100‑88 (MQ=255)                                                                                         tgatgaACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcag            <  2:566514‑M1/100‑90 (MQ=255)                                                                                           atgaACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcagct          >  2:1925187‑M1/1‑9 (MQ=255)                                                                                           atgaACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcagct          >  2:1435700‑M1/1‑9 (MQ=255)                                                                                            tgaACGGCagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcagctc         >  1:1201326‑M1/1‑8 (MQ=255)                                                                                                   cagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcagctcttcgttg  <  1:334877‑M1/100‑100 (MQ=255)                                                                                                   cagcttagccagacgacgggcgatttcagttttaccgacaccggtcgggccgatcatcaggatatttttcggggtcacttcatggcgcagctcttcgttg  <  1:620639‑M1/100‑100 (MQ=255)                                                                                                   |                                                                                                    TTCACGGCGGCATCGGTCAGATCGCGAATAATAGAATCCACTTCCTTACCGACGTAGCCCACTTCGGTGAATTTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGC  >  NC_000913/4121422‑4121619

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