Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,569	A→C	I572L (ATC→CTC)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,162,569	0	A	C	100.0%	111.3 / NA	27	I572L (ATC→CTC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/14);  total (13/14)

TCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  REL606/4162496‑4162640                                                                          |                                                                        tCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTc                                                                       >  2:522245/1‑76 (MQ=255)    agaCGTACACCCGACTCACTACGGTCGCGTATGTACAGTCGAAACCCNTGCAGGTCCGAACATCGGTCTGCTCAAc                                                                    >  1:102363/1‑76 (MQ=255)           cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGt                                                             >  2:364421/1‑76 (MQ=255)                    cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACg                                                    <  2:1117723/76‑1 (MQ=255)                    cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGCCCGTGTACg                                                    >  2:123953/1‑76 (MQ=255)                       tACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGcac                                                 >  1:1140366/1‑76 (MQ=255)                       tACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGcac                                                 >  1:932350/1‑76 (MQ=255)                             cgcgTATGTCCAATCGAAACCCCTGCAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                           <  2:1140365/76‑1 (MQ=255)                                 tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                       >  1:483512/1‑76 (MQ=255)                                   tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                     <  1:456973/76‑1 (MQ=255)                                      ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                  >  2:190482/1‑76 (MQ=255)                                       cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                 >  2:859344/1‑76 (MQ=255)                                       cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                 <  1:341404/76‑1 (MQ=255)                                           cGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                             <  2:404324/76‑1 (MQ=255)                                             aaaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                           >  2:307052/1‑76 (MQ=255)                                               aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                         <  1:776116/76‑1 (MQ=255)                                               aCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                         <  2:770491/76‑1 (MQ=255)                                                 cccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                       <  2:384324/76‑1 (MQ=255)                                                         gTCCGAACATCGGTNTGCTCAACTCTCTGTCCGTNTACGCGCANACTAACGAATACGGCTTCCTTGAGACTCCGTa               <  2:964905/76‑1 (MQ=255)                                                         gTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa               <  1:753054/76‑1 (MQ=255)                                                          tCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt              >  2:1041145/1‑76 (MQ=255)                                                             gAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt           <  2:818576/76‑1 (MQ=255)                                                                 aTCGGTCTGCTCAACTCTCTGTCCGTGGACGCACAGACTAACGCNTANGGCTTCCTTGAGACNCCGTATCGTAAAg       >  2:1054120/1‑76 (MQ=255)                                                                 aTCGCTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg       >  2:1148659/1‑76 (MQ=255)                                                                   cGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg     <  1:533154/76‑1 (MQ=255)                                                                      tCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  <  1:190482/76‑1 (MQ=255)                                                                      tCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc  <  1:548071/76‑1 (MQ=255)                                                                          |                                                                        TCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACC  >  REL606/4162496‑4162640

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