Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,162,570	T→A	I572N (ATC→AAC)	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,570	0	T	A	100.0%	277.9 / NA	69	I572N (ATC→AAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (35/34);  total (35/34)

GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGA  >  REL606/4162498‑4162642                                                                         |                                                                         gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACaa                                                                       >  2:1370599/1‑76 (MQ=255) gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACaa                                                                       >  2:98094/1‑76 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTTAACAActct                                                                   <  1:1680539/76‑1 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActct                                                                   <  2:1325558/76‑1 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActct                                                                   >  2:459008/1‑76 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActct                                                                   <  2:1648349/76‑1 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActct                                                                   >  2:1797889/1‑76 (MQ=255)     cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActca                                                                   <  2:500955/76‑2 (MQ=255)      gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAActctc                                                                  >  1:854044/1‑76 (MQ=255)        acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTg                                                                >  2:1726071/1‑76 (MQ=255)               aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCgtgt                                                         <  1:331524/76‑1 (MQ=255)                cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTa                                                        <  2:1466795/76‑1 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTAca                                                      <  1:1704563/76‑2 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                      >  1:1269628/1‑76 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                      >  2:574050/1‑76 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                      >  1:1593438/1‑76 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                      >  1:1381048/1‑76 (MQ=255)                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACg                                                      <  2:749754/76‑1 (MQ=255)                    cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGcc                                                    >  2:1185474/1‑75 (MQ=255)                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGa                                                >  1:793107/1‑76 (MQ=255)                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTa                                             <  1:1282816/76‑1 (MQ=255)                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTa                                             <  1:1408996/76‑1 (MQ=255)                            gcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTaa                                            <  1:992175/76‑1 (MQ=255)                             cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAAc                                           >  1:1774550/1‑76 (MQ=255)                             cgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAAc                                           >  2:866104/1‑76 (MQ=255)                              gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACg                                          <  1:743253/76‑1 (MQ=255)                              gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACg                                          <  1:589078/76‑1 (MQ=255)                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGa                                         >  1:1017267/1‑76 (MQ=255)                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGa                                         >  1:283119/1‑76 (MQ=255)                                 tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                       <  1:53039/76‑1 (MQ=255)                                 tGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                                       >  2:1107940/1‑76 (MQ=255)                                  gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                      <  1:991375/76‑1 (MQ=255)                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                     <  1:440449/76‑1 (MQ=255)                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                     <  1:1224365/76‑1 (MQ=255)                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                     <  1:1914059/76‑1 (MQ=255)                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                     <  1:954076/76‑1 (MQ=255)                                    ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                    >  1:1172368/1‑76 (MQ=255)                                      aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                  >  1:1415986/1‑76 (MQ=255)                                      aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                  <  1:116186/76‑1 (MQ=255)                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                <  1:505195/76‑1 (MQ=255)                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                <  1:862250/76‑1 (MQ=255)                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                <  1:305034/76‑1 (MQ=255)                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                <  2:459407/76‑1 (MQ=255)                                         cGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                               >  1:549714/1‑76 (MQ=255)                                         cGAAACCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc                               >  1:1093762/1‑76 (MQ=255)                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGACTACGGCTTCCt                             >  2:1504406/1‑76 (MQ=255)                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                             >  1:1852809/1‑76 (MQ=255)                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                             >  2:1504390/1‑76 (MQ=255)                                              ccccTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTga                          >  2:636891/1‑76 (MQ=255)                                                 cTGAAGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc                       >  1:105321/1‑76 (MQ=255)                                                    aaGGTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTTAGACTcc                    >  1:1051633/1‑76 (MQ=255)                                                      ggTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt                  <  1:925579/76‑1 (MQ=255)                                                      ggTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt                  <  2:340375/76‑1 (MQ=255)                                                      ggTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGt                  >  2:587729/1‑76 (MQ=255)                                                       gTCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTa                 <  2:79274/76‑1 (MQ=255)                                                        tCCGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt                >  2:1100255/1‑76 (MQ=255)                                                         ccGAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATc               >  1:340475/1‑76 (MQ=255)                                                           gAACATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt             >  2:1289917/1‑76 (MQ=255)                                                              cATCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa          <  1:1379641/76‑1 (MQ=255)                                                               aTCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg         >  2:1683263/1‑76 (MQ=255)                                                                tCGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGt        <  1:256962/76‑1 (MQ=255)                                                                 cGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg       <  1:1251055/76‑1 (MQ=255)                                                                 cGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg       >  1:591143/1‑76 (MQ=255)                                                                 cGGTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg       >  1:69079/1‑76 (MQ=255)                                                                  ggTCTGAACAACTCTCTGTCCGTGTACGCTCAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa      <  2:330607/76‑1 (MQ=255)                                                                  ggTCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa      <  2:35964/76‑1 (MQ=255)                                                                    tCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc    <  1:597792/76‑1 (MQ=255)                                                                    tCTGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc    <  1:134018/76‑1 (MQ=255)                                                                      tGAACAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGa  >  1:458849/1‑76 (MQ=255)                                                                         |                                                                         GAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGA  >  REL606/4162498‑4162642

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