Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	4,168,970	Δ3 bp	coding (4010‑4012/4224 nt)	rpoC →	DNA‑directed RNA polymerase subunit beta'

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	4168970	4168972	3	33 [0]	[0] 33	rpoC	DNA‑directed RNA polymerase subunit beta'

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 4168969	0 (0.000)	33 (0.850)	24/148	0.6	100%	coding (4009/4224 nt)	rpoC	DNA‑directed RNA polymerase subunit beta'
	?	REL606	4168973 =	0 (0.000)					coding (4013/4224 nt)	rpoC	DNA‑directed RNA polymerase subunit beta'

GTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGG  >  REL606/4168896‑4169039                                                                          |                                                                       gTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcg                                                                      >  1:72239‑M1/1‑74 (MQ=255)   ctctGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcggg                                                                    <  2:717323‑M1/76‑5 (MQ=255)     ctGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggca                                                                  >  2:36065‑M1/1‑70 (MQ=255)          aaCCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgc                                                             >  2:603461‑M1/1‑65 (MQ=255)          aaCCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgc                                                             >  1:601675‑M1/1‑65 (MQ=255)              gAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacg                                                         <  1:8349‑M1/76‑16 (MQ=255)              gAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacg                                                         <  2:409996‑M1/76‑16 (MQ=255)              gAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacg                                                         <  1:921914‑M1/76‑16 (MQ=255)              gAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacg                                                         >  1:1023710‑M1/1‑61 (MQ=255)                 tCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaac                                                      >  2:1065353‑M1/1‑58 (MQ=255)                 tCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaac                                                      <  2:1136312‑M1/76‑19 (MQ=255)                  ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaact                                                     >  1:458946‑M1/1‑57 (MQ=255)                  ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaact                                                     >  1:39855‑M1/1‑57 (MQ=255)                           ccGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcct                                            <  2:484281‑M1/76‑29 (MQ=255)                            cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGTCCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctg                                           <  2:628719‑M1/76‑30 (MQ=255)                            cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAATCAGCCGcgggcaaacgcgacgaactgcgcggcctg                                           <  2:1136748‑M1/76‑30 (MQ=255)                            cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctg                                           <  2:550702‑M1/76‑30 (MQ=255)                             gcgGCATCGTTCCAGGAGACCGCTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctga                                          <  2:499273‑M1/76‑31 (MQ=255)                             gcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctga                                          <  2:1269810‑M1/76‑31 (MQ=255)                              cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaa                                         <  1:14101‑M1/76‑32 (MQ=255)                               ggCATCGTTCCAGGAGCCCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaa                                        <  1:444528‑M1/76‑33 (MQ=255)                               ggCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaa                                        <  2:315120‑M1/76‑33 (MQ=255)                                  aTCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagag                                     <  1:1065353‑M1/76‑36 (MQ=255)                                   tCGTTCCAGGAGACCACTCGCGTGGTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagaga                                    <  1:761319‑M1/76‑37 (MQ=255)                                    cGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaagcgaaacgaactgggcggcctgatagagga                                   >  2:456473‑M1/1‑39 (MQ=255)                                      ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacg                                 >  2:1175314‑M1/1‑37 (MQ=255)                                       tCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacacgccgaacaccccggcccgacagagaacgt                                >  2:352122‑M1/1‑36 (MQ=255)                                          aGGAGACCCCTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttat                             <  2:1018912‑M1/76‑44 (MQ=255)                                                ccACTCGCGTGCTGACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtggg                       >  1:1137991‑M1/1‑27 (MQ=255)                                                             gACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtgggtcgtctgatcccg          >  1:68654‑M1/1‑14 (MQ=255)                                                             gACCGAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtgggtcgtctgatcccg          >  2:753937‑M1/1‑14 (MQ=255)                                                                 gAAGCAGCCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtgggacgtctgatcccggcag      >  2:999065‑M1/1‑10 (MQ=255)                                                                     cagcCGcgggcaaacgcgacgaactgcgcggcctgaaagagaacgttatcgtgggtcgtctgatcccggcaggtac  >  2:1084604‑M1/1‑6 (MQ=255)                                                                          |                                                                       GTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGG  >  REL606/4168896‑4169039

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