Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	3,014,068	(AAC)3→2	coding (649‑651/1137 nt)	ECB_02815 →	Capsule polysaccharide export inner‑membrane protein kpsE

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	3014060	3014062	3	70 [0]	[0] 71	ECB_02815	Capsule polysaccharide export inner‑membrane protein kpsE

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	3014063 =	0 (0.000)	54 (0.950)	41/134	0.3	100%	coding (644/1137 nt)	ECB_02815	Capsule polysaccharide export inner‑membrane protein kpsE
	?	REL606	= 3014059	0 (0.000)					coding (640/1137 nt)	ECB_02815	Capsule polysaccharide export inner‑membrane protein kpsE

TTGCAGAAACGGAGATGGAAAAGGCACGCCAGCGTCTGGACGCCAGCAAAGCGGAATTGCTCTCTTATCAAGACAACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGCCAGAAG  >  REL606/3013988‑3014136                                                                            |                                                                          cctttgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaaga                                                                           >  2:896112‑M1/76‑76 (MQ=255) cctttgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaaga                                                                           <  2:614636‑M1/1‑1 (MQ=255) cctttgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaaga                                                                           <  1:1309853‑M1/1‑1 (MQ=255)    ttgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaa                                                                        >  2:463728‑M1/73‑76 (MQ=255)    ttgcagaaacggaagggaaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaa                                                                        >  1:905822‑M1/73‑76 (MQ=255)     tgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaac                                                                       >  1:1400330‑M1/72‑76 (MQ=255)     tgcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaac                                                                       >  1:1782237‑M1/72‑76 (MQ=255)      gcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaaca                                                                      >  2:59940‑M1/71‑76 (MQ=255)      gcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaaca                                                                      >  2:525553‑M1/71‑76 (MQ=255)      gcagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaaca                                                                      >  2:525535‑M1/71‑76 (MQ=255)       cagaaccggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagacaacan                                                                     >  1:856053‑M1/70‑75 (MQ=255)       cagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagAcaacaa                                                                     >  2:409930‑M1/70‑76 (MQ=255)       cagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagAcaacaa                                                                     >  2:1731181‑M1/70‑76 (MQ=255)       cagaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagAcaacaa                                                                     >  2:1981433‑M1/70‑76 (MQ=255)        agaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACaacaac                                                                    >  1:1516882‑M1/69‑76 (MQ=255)        agaaacggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACaacaac                                                                    >  1:1988008‑M1/69‑76 (MQ=255)          aaagggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGt                                                                  >  1:609240‑M1/67‑76 (MQ=255)            acggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTc                                                                >  1:1678926‑M1/65‑76 (MQ=255)              ggagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTg                                                              >  2:320375‑M1/63‑76 (MQ=255)               gagatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTgg                                                             >  1:1700616‑M1/62‑76 (MQ=255)                agatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGa                                                            <  1:1050593‑M1/16‑1 (MQ=255)                agatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGa                                                            <  2:592585‑M1/16‑1 (MQ=255)                agatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGa                                                            >  1:1402727‑M1/61‑76 (MQ=255)                 gatggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGCTCTGGAt                                                           >  2:436543‑M1/60‑76 (MQ=255)                 gatgaaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGAt                                                           >  2:567613‑M1/60‑76 (MQ=255)                  gtggaaaagccacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATc                                                          >  1:499599‑M1/59‑76 (MQ=255)                    ggaaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAAAGTTATggataca                                                        <  1:976332‑M1/20‑4 (MQ=255)                      aaaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcaca                                                      <  1:1075135‑M1/22‑1 (MQ=255)                       aaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacag                                                     >  2:1803374‑M1/54‑76 (MQ=255)                       aaaggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacag                                                     <  2:278138‑M1/23‑1 (MQ=255)                         aggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggc                                                   >  1:1452132‑M1/52‑76 (MQ=255)                         aggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggc                                                   <  1:1562711‑M1/25‑1 (MQ=255)                          ggcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggca                                                  >  2:1879261‑M1/51‑76 (MQ=255)                           gcacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcac                                                 <  1:495282‑M1/27‑1 (MQ=255)                            cacgccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcaca                                                <  2:742379‑M1/28‑1 (MQ=255)                                ccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcacaggca                                            >  2:235854‑M1/45‑76 (MQ=255)                                ccagcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcacaggca                                            <  2:1811005‑M1/32‑1 (MQ=255)                                  agcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcacaggcaca                                          >  2:1898059‑M1/43‑76 (MQ=255)                                  agcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcacaggcaca                                          <  2:397111‑M1/34‑1 (MQ=255)                                  agcgtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCcacaggcacaggcaca                                          >  2:1821440‑M1/43‑76 (MQ=255)                                     gtctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCacaggcacaggcacaggc                                       >  2:111643‑M1/40‑76 (MQ=255)                                      tctggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAggcg                                      <  2:1556618‑M1/38‑1 (MQ=255)                                         ggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGgcggcg                                   <  1:1792805‑M1/41‑1 (MQ=255)                                         ggacgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGgcggcg                                   >  2:415867‑M1/36‑76 (MQ=255)                                            cgccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGc                                >  1:764078‑M1/33‑76 (MQ=255)                                             gccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCa                               >  2:822444‑M1/32‑76 (MQ=255)                                             gccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCa                               >  2:495282‑M1/32‑76 (MQ=255)                                             gccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCa                               >  1:787554‑M1/32‑76 (MQ=255)                                             gccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCa                               >  1:1835978‑M1/32‑76 (MQ=255)                                              ccagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCAc                              >  1:1454691‑M1/31‑76 (MQ=255)                                               cagcaaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACg                             <  2:74920‑M1/47‑1 (MQ=255)                                                  caaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTa                          <  1:579633‑M1/50‑1 (MQ=255)                                                   acagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGCGCACGTTAg                         <  2:378597‑M1/51‑1 (MQ=255)                                                   aaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAg                         <  1:402121‑M1/51‑1 (MQ=255)                                                   aaagcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAg                         <  2:1636930‑M1/51‑1 (MQ=255)                                                    atncggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGt                        <  2:1420201‑M1/52‑1 (MQ=255)                                                     agcggaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTg                       >  1:14248‑M1/24‑76 (MQ=255)                                                         gaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATa                   <  2:634814‑M1/57‑1 (MQ=255)                                                         gaattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATa                   <  1:1300111‑M1/57‑1 (MQ=255)                                                          aattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATAc                  <  2:704710‑M1/58‑1 (MQ=255)                                                          aattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATAc                  <  1:1045451‑M1/58‑1 (MQ=255)                                                          aattgctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATAc                  <  1:1757445‑M1/58‑1 (MQ=255)                                                             tcctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCt               >  1:278138‑M1/16‑76 (MQ=255)                                                              gctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTg              >  1:1337972‑M1/15‑76 (MQ=255)                                                              gctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTg              <  1:823781‑M1/62‑1 (MQ=255)                                                               ctctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGa             <  1:422979‑M1/63‑1 (MQ=255)                                                                 ctcttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATg           <  2:82250‑M1/65‑1 (MQ=255)                                                                    ttatcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGc        >  2:316239‑M1/9‑76 (MQ=255)                                                                      atcaagACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGCCa      >  1:930316‑M1/7‑76 (MQ=255)                                                                         aacacagcaacGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGCCagaa   <  1:399205/73‑1 (MQ=255)                                                                          agACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGCCagaag  >  1:1903724‑M1/3‑76 (MQ=255)                                                                            |                                                                          TTGCAGAAACGGAGATGGAAAAGGCACGCCAGCGTCTGGACGCCAGCAAAGCGGAATTGCTCTCTTATCAAGACAACAACAACGTTCTGGATCCACAGGCACAGGCACAGGCGGCGAGCACGTTAGTGAATACGCTGATGGGCCAGAAG  >  REL606/3013988‑3014136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑