Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,163,896	T→C	L1014P (CTG→CCG)	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,163,896	0	T	C	100.0%	137.8 / NA	34	L1014P (CTG→CCG)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/13);  total (21/13)

GACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCA  >  REL606/4163823‑4163968                                                                          |                                                                         gACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAgccgg                                                                        >  2:290599/1‑76 (MQ=255)    aaaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAgccggctg                                                                     <  2:163286/76‑1 (MQ=255)     aaCTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGa                                                                    >  1:52148/1‑76 (MQ=255)         gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAg                                                                >  1:1043415/1‑76 (MQ=255)         gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAg                                                                >  1:1204195/1‑76 (MQ=255)         gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAg                                                                >  2:1516042/1‑76 (MQ=255)         gCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAg                                                                >  1:284307/1‑76 (MQ=255)          ccGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGt                                                               <  2:1573017/76‑1 (MQ=255)            gcgcgATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTAt                                                             <  2:442147/76‑1 (MQ=255)            gcgcgATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTAt                                                             <  1:41236/76‑1 (MQ=255)             cgcgATCGCTGGCTGGAGCTGGGCCTGACTGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATg                                                            <  2:1126289/76‑1 (MQ=255)             cgcgATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATg                                                            <  1:1271452/76‑1 (MQ=255)             cgcgATCGCGGGCTGGAGCTGGGCCTGGCCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATg                                                            <  2:776064/76‑1 (MQ=255)              gcgATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGa                                                           <  1:619269/76‑1 (MQ=255)                       ggctggAGCTGGGCCTGACCGACAAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGaa                                                  <  1:231591/76‑1 (MQ=255)                            gAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACg                                             >  2:319804/1‑76 (MQ=255)                                 gggCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTc                                        >  2:1024274/1‑76 (MQ=255)                                   gCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCga                                      >  1:604508/1‑76 (MQ=255)                                        aCCGGCGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCgggatgt                                 <  2:170913/76‑4 (MQ=255)                                         ccGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAgaagaa                                >  2:1297029/1‑76 (MQ=255)                                         ccGACGAAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAgaagaa                                >  1:1037818/1‑76 (MQ=255)                                          cgacgaAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGaaa                               <  2:723787/76‑1 (MQ=255)                                            acgaAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACt                             <  2:1024403/76‑1 (MQ=255)                                            acgaAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACt                             >  1:749590/1‑76 (MQ=255)                                              gaAGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCg                           >  1:1322845/1‑76 (MQ=255)                                               aaGAGAAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGa                          >  1:1203828/1‑76 (MQ=255)                                                  agaAACAAAATCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGc                       >  2:1042435/1‑76 (MQ=255)                                                      aCAAAATCAGCGGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGaaa                   <  1:176040/76‑1 (MQ=255)                                                        aaaaTCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcg                 >  1:865770/1‑76 (MQ=255)                                                         aaaTCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcgc                >  1:1338313/1‑76 (MQ=255)                                                         aaaTCAGCTGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAAcgc                >  1:822475/1‑76 (MQ=255)                                                                 tGGAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAAt        >  2:176142/1‑76 (MQ=255)                                                                   gAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCa      >  2:1027803/1‑76 (MQ=255)                                                                   gAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCa      >  2:46661/1‑76 (MQ=255)                                                                   gAACAGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCa      <  1:646045/76‑1 (MQ=255)                                                                       aGCCGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCa  >  1:855573/1‑76 (MQ=255)                                                                          |                                                                         GACAAACTGCCGCGCGATCGCTGGCTGGAGCTGGGCCTGACCGACGAAGAGAAACAAAATCAGCTGGAACAGCTGGCTGAGCAGTATGACGAACTGAAACACGAGTTCGAGAAGAAACTCGAAGCGAAACGCCGCAAAATCACCCA  >  REL606/4163823‑4163968

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