Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,171,725	G→C	D185H (GAC→CAC)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,171,725	0	G	C	100.0%	349.2 / NA	92	D185H (GAC→CAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (42/50);  total (42/50)

TCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171647‑4171809                                                                               |                                                                                     tCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTc                                                                           <  2:1691892/90‑1 (MQ=255)   cGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGa                                                                         <  2:928715/90‑1 (MQ=255)      cAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATcc                                                                      <  2:1407711/90‑1 (MQ=255)       aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCg                                                                     <  2:1908782/90‑1 (MQ=255)       aaaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCg                                                                     <  2:1222807/90‑1 (MQ=255)        aaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGa                                                                    <  1:1709760/90‑1 (MQ=255)         aGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGaa                                                                   <  1:1366747/90‑1 (MQ=255)         aGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGaa                                                                   <  1:2130107/90‑1 (MQ=255)         aGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGaa                                                                   <  2:1052916/90‑1 (MQ=255)             aaaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  1:567758/90‑1 (MQ=255)             aaaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  2:2048068/90‑1 (MQ=255)              aaaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACa                                                              >  2:569294/1‑90 (MQ=255)               aaCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACaa                                                             <  1:1110438/90‑1 (MQ=255)                 cccACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  2:651952/90‑1 (MQ=255)                   cACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTg                                                         >  1:1471923/1‑90 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAATGACAACCTGt                                                        >  2:1495731/1‑90 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               >  1:639056/1‑83 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  2:639056/83‑1 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:1932384/1‑90 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        <  2:1277957/90‑1 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:1553281/1‑90 (MQ=255)                    aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        <  1:2141847/90‑1 (MQ=255)                        ttCGGGTAAAGTTCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:2133991/90‑1 (MQ=255)                        ttCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  2:1553281/90‑1 (MQ=255)                        ttCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTACTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  2:1494210/90‑1 (MQ=255)                              tAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTAt                                              >  2:1395103/1‑90 (MQ=255)                              tAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTAt                                              >  2:1306838/1‑90 (MQ=255)                               aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATc                                             >  2:994737/1‑90 (MQ=255)                               aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATc                                             >  1:360303/1‑90 (MQ=255)                                aaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCg                                            <  1:1269937/90‑1 (MQ=255)                                 aGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGa                                           >  1:33691/1‑90 (MQ=255)                                 aGTGCTGTATAACGAGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    >  2:1757022/1‑81 (MQ=255)                                 aGTGCTGTATAACGAGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:1757022/81‑1 (MQ=255)                                     cTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:2141618/90‑1 (MQ=255)                                     cTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:311397/90‑1 (MQ=255)                                      tGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTc                                      >  1:347223/1‑90 (MQ=255)                                       gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACaa                                                             <  1:1176067/66‑1 (MQ=255)                                       gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACaa                                                             >  2:1176067/1‑66 (MQ=255)                                           aaCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgta                                                   >  2:2017123/1‑72 (MQ=255)                                           aaCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgta                                                   <  1:2017123/72‑1 (MQ=255)                                            aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCt                                                          <  2:1689104/64‑1 (MQ=255)                                            aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCt                                                          >  1:1689104/1‑64 (MQ=255)                                            aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  1:535000/1‑90 (MQ=255)                                              gcgcgTATAATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCg                                                                          <  1:640035/46‑1 (MQ=255)                                              gcgcgTATAATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCg                                                                          >  2:640035/1‑46 (MQ=255)                                               cgggTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAc                             <  1:130173/87‑1 (MQ=255)                                               cgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCTGTAAAc                             <  1:1855578/90‑1 (MQ=255)                                               cgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAc                             >  1:308546/1‑90 (MQ=255)                                                 cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActg                           >  2:1021028/1‑90 (MQ=255)                                                 cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAc                             <  2:742466/88‑1 (MQ=255)                                                 cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAc                             >  1:742466/1‑88 (MQ=255)                                                    atcatcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActgcct                        <  2:308546/90‑1 (MQ=255)                                                     tcatcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActgcctg                       <  1:794122/90‑1 (MQ=255)                                                     tcatcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActgcctg                       <  1:1328906/90‑1 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               >  2:2124315/1‑48 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  1:2124315/48‑1 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       >  1:800824/1‑72 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:800824/72‑1 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     >  1:1883220/1‑90 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  2:347223/90‑1 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  2:414349/90‑1 (MQ=255)                                                       atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  2:1971219/90‑1 (MQ=255)                                                           cTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACca                 >  2:1390829/1‑90 (MQ=255)                                                            ttACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    >  2:277947/1‑54 (MQ=255)                                                            ttACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:277947/54‑1 (MQ=255)                                                             tACCGTGGTTCCTTGCTGCACTTCGAATTCAATCCGAAGGACAACCTGGTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatc               >  1:804091/1‑90 (MQ=255)                                                             tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    >  2:2072837/1‑53 (MQ=255)                                                             tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:2072837/53‑1 (MQ=255)                                                             tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatc               >  2:585474/1‑90 (MQ=255)                                                             tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatc               >  2:485977/1‑90 (MQ=255)                                                             tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatc               >  2:1612469/1‑90 (MQ=255)                                                              aCCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcatca              <  1:835784/90‑1 (MQ=255)                                                                cGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:2009116/63‑1 (MQ=255)                                                                cGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       >  1:2009116/1‑63 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCTATCCGAAGGACAACCTGTTc                                                      >  2:1463033/1‑45 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCTATCCGAAGGACAACCTGTTc                                                      <  1:1463033/45‑1 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCTATCCGAAGGACAACCTGTTc                                                      >  2:1669522/1‑45 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCTATCCGAAGGACAACCTGTTc                                                      <  1:1669522/45‑1 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        <  2:319399/43‑1 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        <  2:222974/43‑1 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:222974/1‑43 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:319399/1‑43 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTc                                                      <  1:1822289/45‑1 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTc                                                      >  2:1822289/1‑45 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTg         >  2:1203311/1‑90 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTg         >  1:836681/1‑90 (MQ=255)                                                                   ggTTCCTGGCTGCACTTCGAATTCGATCCAAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTg         >  2:280415/1‑90 (MQ=255)                                                                     ttCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcg       <  2:1932384/90‑1 (MQ=255)                                                                      tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  2:360303/90‑1 (MQ=255)                                                                      tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      >  2:986202/1‑90 (MQ=255)                                                                      tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      >  2:794979/1‑90 (MQ=255)                                                                      tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  1:1295503/90‑1 (MQ=255)                                                                        ctggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTaaa                              <  2:607913/64‑1 (MQ=255)                                                                        ctggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTaaa                              >  1:607913/1‑64 (MQ=255)                                                                        ctggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgcgc    <  2:836681/90‑1 (MQ=255)                                                                         tggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGcc   >  2:428211/1‑90 (MQ=255)                                                                          ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  1:1252947/90‑1 (MQ=255)                                                                          ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  2:1502959/90‑1 (MQ=255)                                                                               |                                                                                     TCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171647‑4171809

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