Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,161,283	G→T	R143L (CGT→CTT)	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,161,283	0	G	T	100.0%	231.6 / NA	59	R143L (CGT→CTT)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (27/32);  total (27/32)

AACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTT  >  REL606/4161211‑4161354                                                                         |                                                                        aaCAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCttgt                                                                      >  1:794711/1‑76 (MQ=255)  aCAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGtt                                                                     <  1:1620625/76‑1 (MQ=255)   cAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTa                                                                    >  1:305203/1‑76 (MQ=255)   cAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTa                                                                    >  2:1500487/1‑76 (MQ=255)   cAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTa                                                                    >  2:1149028/1‑76 (MQ=255)    aagaagTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTAt                                                                   >  2:133828/1‑76 (MQ=255)    aagaagTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTAt                                                                   >  1:1272/1‑76 (MQ=255)       aagTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGt                                                                >  1:659491/1‑76 (MQ=255)        agTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGATTGTTATCGtt                                                               >  1:232401/1‑76 (MQ=255)           cTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTcc                                                            >  1:1639349/1‑76 (MQ=255)             aCATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCa                                                          <  1:527779/76‑1 (MQ=255)              cATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAg                                                         <  2:80984/76‑1 (MQ=255)              cATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAg                                                         >  2:1715412/1‑76 (MQ=255)               aTGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGc                                                        >  1:951717/1‑76 (MQ=255)                tGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCt                                                       <  1:996025/76‑1 (MQ=255)                 gggCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTg                                                      >  1:1114296/1‑76 (MQ=255)                  ggCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGc                                                     <  1:837265/76‑1 (MQ=255)                    cGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCAc                                                   <  1:299219/76‑1 (MQ=255)                      aaaTTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCg                                                 >  1:1105471/1‑76 (MQ=255)                         ttCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAg                                              <  1:130530/76‑1 (MQ=255)                         ttCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGATTGTTATCGTTTCCCAGCTGCACCGTAg                                              <  2:1472607/76‑1 (MQ=255)                           ccGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTc                                            >  2:1351263/1‑76 (MQ=255)                            cGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTcc                                           <  1:982978/76‑1 (MQ=255)                             gCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCg                                          <  2:795617/76‑1 (MQ=255)                              cTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCgg                                         >  1:879659/1‑76 (MQ=255)                                cATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGc                                       <  2:779135/76‑1 (MQ=255)                                     cagacaACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtctt                                  <  2:1671620/76‑1 (MQ=255)                                     cagacaACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtctt                                  >  1:1555870/1‑76 (MQ=255)                                     cagacaACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtctt                                  <  2:1103365/76‑1 (MQ=255)                                       gacaACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGtcttct                                <  1:1215171/76‑1 (MQ=255)                                         caACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCttt                              <  1:416394/76‑1 (MQ=255)                                         caACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCttt                              <  2:1324020/76‑1 (MQ=255)                                         caACGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCttt                              <  2:811489/76‑1 (MQ=255)                                          aaCGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTg                             >  2:402764/1‑76 (MQ=255)                                          aaCGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTg                             >  2:1604086/1‑76 (MQ=255)                                          aaCGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTg                             >  2:212182/1‑76 (MQ=255)                                           aCGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGa                            <  2:63575/76‑1 (MQ=255)                                            cGGTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGAc                           <  1:478360/76‑1 (MQ=255)                                              gTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTc                         <  1:619859/76‑1 (MQ=255)                                              gTACCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTc                         <  2:715034/76‑1 (MQ=255)                                                aCCTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCg                       >  1:1461554/1‑76 (MQ=255)                                                  cTTTGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGAc                     <  1:1038800/76‑1 (MQ=255)                                                   tttGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACa                    >  1:823646/1‑76 (MQ=255)                                                    ntGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACaa                   >  1:1112846/2‑76 (MQ=255)                                                     tGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACaaa                  >  2:743973/1‑76 (MQ=255)                                                     tGTTATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACaaa                  <  2:557379/76‑1 (MQ=255)                                                       ttATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAgg                <  1:1518765/76‑1 (MQ=255)                                                        tATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTTTTCTTTGACTCCGACAAAGGt               <  2:967091/76‑1 (MQ=255)                                                        tATCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTTTTCTTTGACTCCGACAAAGGt               <  2:967081/76‑1 (MQ=255)                                                         aTCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTa              <  2:528165/76‑1 (MQ=255)                                                         aTCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTa              >  1:188041/1‑76 (MQ=255)                                                         aTCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTa              <  1:368175/76‑1 (MQ=255)                                                          tCAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaa             >  2:1164307/1‑76 (MQ=255)                                                           cAACGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaaa            >  1:347257/1‑76 (MQ=255)                                                            aaCGGTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTaaaa           >  2:563787/1‑76 (MQ=255)                                                                gTACTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCa       >  2:751729/1‑76 (MQ=255)                                                                   cTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCActc    <  2:1493726/76‑1 (MQ=255)                                                                   cTGAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCActc    >  1:473909/1‑76 (MQ=255)                                                                     gAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCtt  <  2:290992/76‑1 (MQ=255)                                                                     gAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCtt  <  1:1132492/76‑1 (MQ=255)                                                                     gAGCTTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCtt  <  1:1740084/76‑1 (MQ=255)                                                                         |                                                                        AACAAGAAGTCTACATGGGCGAAATTCCGCTCATGACAGACAACGGTACCTTTGTTATCAACGGTACTGAGCGTGTTATCGTTTCCCAGCTGCACCGTAGTCCGGGCGTCTTCTTTGACTCCGACAAAGGTAAAACCCACTCTT  >  REL606/4161211‑4161354

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