Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	REL606	4,163,094	Δ3 bp	coding (2239‑2241/4029 nt)	rpoB →	DNA‑directed RNA polymerase subunit beta

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	REL606	4163094	4163096	3	54 [0]	[0] 55	rpoB	DNA‑directed RNA polymerase subunit beta

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	REL606	= 4163093	0 (0.000)	54 (1.010)	42/150	0.2	100%	coding (2238/4029 nt)	rpoB	DNA‑directed RNA polymerase subunit beta
	?	REL606	4163097 =	0 (0.000)					coding (2242/4029 nt)	rpoB	DNA‑directed RNA polymerase subunit beta

CGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTC  >  REL606/4163024‑4163172                                                                          |                                                                            tgtcgttcagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATc                                                                           <  2:950730‑M2/3‑1 (MQ=255)  gtcgttcagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCg                                                                          >  2:1030267‑M2/73‑76 (MQ=255)    cgttcagtacgtggatgcttcccgtatcgctatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGAc                                                                        <  2:292156‑M2/6‑1 (MQ=255)      ttcagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACAt                                                                      >  1:89233‑M2/69‑76 (MQ=255)       tcagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATc                                                                     >  1:899270‑M2/68‑76 (MQ=255)        cagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCt                                                                    <  1:429222‑M2/10‑1 (MQ=255)        cagtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagaaATCGACATCt                                                                    >  2:190190‑M2/67‑76 (MQ=255)         agtacgtggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTa                                                                   >  1:404625‑M2/66‑76 (MQ=255)               tggatgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCt                                                             >  2:1580900‑M2/60‑76 (MQ=255)                  atgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGAc                                                          >  1:620815‑M2/57‑76 (MQ=255)                  atgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGAc                                                          <  1:1060440‑M2/20‑1 (MQ=255)                   tgcttcccgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGAcc                                                         <  1:1331593‑M2/21‑1 (MQ=255)                          cgtatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATaca                                                  <  2:1548577‑M2/28‑1 (MQ=255)                            tatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACAcc                                                >  2:616868‑M2/47‑76 (MQ=255)                            tatcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACAcc                                                >  2:792543‑M2/47‑76 (MQ=255)                             atcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACAccc                                               >  1:342797‑M2/46‑76 (MQ=255)                             atcgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACAccc                                               >  2:1337698‑M2/46‑76 (MQ=255)                               cgttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGt                                             >  1:67174‑M2/44‑76 (MQ=255)                                gttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGtt                                            >  1:1233402‑M2/43‑76 (MQ=255)                                gttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGtt                                            <  2:899088‑M2/34‑1 (MQ=255)                                 ttatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTc                                           >  2:695701‑M2/42‑76 (MQ=255)                                  tatcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCt                                          >  1:1291712‑M2/41‑76 (MQ=255)                                   atcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTa                                         >  2:1049343‑M2/40‑76 (MQ=255)                                   atcaaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTa                                         <  1:12292‑M2/37‑1 (MQ=255)                                      aaagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAAcc                                      >  2:1074510‑M2/37‑76 (MQ=255)                                       aagttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCa                                     >  2:1203931‑M2/36‑76 (MQ=255)                                        agttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAg                                    <  1:696159‑M2/42‑1 (MQ=255)                                         gttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGGCCAAATACACCCGTTCTAACCAGa                                   <  2:699627‑M2/43‑1 (MQ=255)                                         gttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGa                                   <  1:179742‑M2/43‑1 (MQ=255)                                         gttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGa                                   <  2:699616‑M2/43‑1 (MQ=255)                                         gttaacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGa                                   <  2:1650460‑M2/43‑1 (MQ=255)                                           taacgaagacgagatgtatccgggtgacgcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAac                                 >  2:1512761‑M2/32‑76 (MQ=255)                                            aacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAaca                                <  2:1206142‑M2/46‑1 (MQ=255)                                            aacgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAaca                                >  1:1349772‑M2/31‑76 (MQ=255)                                             acgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAacac                               <  2:1380216‑M2/47‑1 (MQ=255)                                             acgaagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAacac                               <  2:375990‑M2/47‑1 (MQ=255)                                                aagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTg                            <  2:801193‑M2/50‑1 (MQ=255)                                                aagacgagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTg                            <  1:1545712‑M2/50‑1 (MQ=255)                                                     gagatgtatccgggtgaagcgATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCa                       <  1:557088‑M2/55‑1 (MQ=255)                                                     gagatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCa                       <  1:151737‑M2/55‑1 (MQ=255)                                                      agatgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCaa                      <  1:1497081‑M2/56‑1 (MQ=255)                                                        atgtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAAcc                    >  1:43480‑M2/19‑76 (MQ=255)                                                          gtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAg                  <  1:430011‑M2/60‑1 (MQ=255)                                                          gtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAg                  <  2:1299743‑M2/60‑1 (MQ=255)                                                          gtatccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAg                  >  1:495274‑M2/17‑76 (MQ=255)                                                              ccgggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGc              <  2:1129927‑M2/64‑1 (MQ=255)                                                                gggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCg            >  2:1607841‑M2/11‑76 (MQ=255)                                                                gggtgaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCg            >  2:763664‑M2/11‑76 (MQ=255)                                                                    gaagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCgtgtg        <  2:686506‑M2/70‑1 (MQ=255)                                                                     aagcaATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCgtgtgt       >  2:1349999‑M2/6‑76 (MQ=255)                                                                        caATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtc    >  1:1269408‑M2/3‑76 (MQ=255)                                                                        caATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtc    >  1:1233467‑M2/3‑76 (MQ=255)                                                                         aATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtct   >  1:1621195‑M2/2‑76 (MQ=255)                                                                         aATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtct   >  1:814791‑M2/2‑76 (MQ=255)                                                                          aTCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtcta  <  1:1361423/76‑2 (MQ=255)                                                                          |                                                                            CGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTC  >  REL606/4163024‑4163172

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