Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,543	C→A	P1100Q (CCG→CAG)	rpoB →	RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,184,543	0	C	A	100.0%	177.9 / NA	56	P1100Q (CCG→CAG)	rpoB	RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (36/20);  total (36/20)

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184405‑4184678                                                                                                                                           |                                                                                                                                        tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                     <  1:580580/143‑1 (MQ=255)    cGGTTAAGCGCCGTAGCGCGTCTGGTGACAAGATGGCATGTCGTCACGGTAACAAGGGTGTAATTTCTAACATCATCCCTTTCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                  <  2:722943/143‑1 (MQ=255)       ttAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTa                                                                                                                               <  1:690443/143‑1 (MQ=255)            cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGAAGACATCGTACTGAACCAGCCGGGCGACCCGtc                                                                                                                          >  2:396152/1‑143 (MQ=255)            cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGCATGTAATTTCTAAGATCAACCCGATAGAAGATATGCCATACGATGAAAACGGAACGGCGCTAAAAATCGGACTAAAACAGCTAGGCgcaacgac                                                                                                                          >  2:127400/1‑136 (MQ=255)             gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtcg                                                                                                                         >  1:307197/1‑142 (MQ=255)             gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGCTAGACATCGTACTGAACCAGCTGGGCTTATCgcaa                                                                                                                         >  2:454895/1‑140 (MQ=255)              cggTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCCACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                        <  2:443105/141‑1 (MQ=255)                 tcgtcAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                     <  2:666735/139‑1 (MQ=255)                    ccAGCCTGGTCAACGGATGGCAGGGCGTCAGGGTAACAAGGGTCTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGa                                                                                                                  <  2:415486/143‑1 (MQ=255)                     cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGaa                                                                                                                 >  1:419320/1‑143 (MQ=255)                        ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                              >  1:277347/1‑143 (MQ=255)                          cggTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                            <  2:769577/142‑1 (MQ=255)                                 aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGGCCCGGCTCGTCTGAACGTCGGTCAGAt                                                                                                     >  2:716464/1‑143 (MQ=255)                                   gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                            >  2:494489/1‑134 (MQ=255)                                   gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                            <  1:494489/134‑1 (MQ=255)                                     tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                 <  1:316611/143‑1 (MQ=255)                                      ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCg                                                                                                >  1:538490/1‑143 (MQ=255)                                       gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGCAGATATAACTTCAGATGAAGATGGTACCCCGGTAGACATCGTACTGAACCAGCTCGGCGTACAGTCTCGCATGCACGTCGGTCCGGTCCTCGa                                                                                               >  1:700632/1‑143 (MQ=255)                                       gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACACGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGAATAAAAAGCGGTCAGAACCTCGa                                                                                               >  2:187000/1‑143 (MQ=255)                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGGATGAACATCGGTCAGACCCTCGCAAcccccc                                                                                       >  2:137891/1‑143 (MQ=255)                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGAATCGTAAGAACAACGGCCAGATCGACGGAAccgcca                                                                                       >  2:176118/1‑139 (MQ=255)                                                cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGACGTACCGTCTCGTATGCAAATCGGGCAGATCCTCGAGACCCAcca                                                                                      >  2:746888/1‑142 (MQ=255)                                                cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTGACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCt                                                                                      <  1:478928/143‑1 (MQ=255)                                                  cGGTAACAAGGGTGTAATTCCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGTACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                    <  2:275689/143‑1 (MQ=255)                                                        cAAGGGTGTAATTTCTAAGCTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                              <  1:693359/143‑1 (MQ=255)                                                        cAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                              <  2:220519/143‑1 (MQ=255)                                                          aGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCt                                                                            <  1:509899/143‑1 (MQ=255)                                                             gtgtAATTTCTAAGTTGAACCGGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCg                                                                         <  2:59904/143‑1 (MQ=255)                                                               gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATGGGTCAGATGCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                       >  2:200449/1‑143 (MQ=255)                                                                  aTTTCTAAGATCAACCCGATCGAAGATATGCCTTNCGATGANAACGGTNCGCCGGTAGACATCGTACNGAACCAGCTGGGCGTACCGTCCCGTCTGANCATCGCTCAGATCCTCGACACCCACCTGGCTATGGCCGCCACAgg                                                                    >  2:405408/1‑143 (MQ=255)                                                                     tCTAAGATCGACCCGATCGAAGATATGCCTTACGATGAGAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                                 <  2:695417/143‑1 (MQ=255)                                                                          gATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCg                                                            >  2:64657/1‑143 (MQ=255)                                                                              cgcccGTTCGAAGACATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                        <  1:205353/141‑1 (MQ=255)                                                                                 ccGATTGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACCAGAt                                                     >  1:698756/1‑143 (MQ=255)                                                                                 ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaggga                                                     >  2:158695/1‑139 (MQ=255)                                                                                   gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCCCCTGGGTGTGGCTGCGAGAGGTATCGGGGACAAGAGCa                                                   >  2:611613/1‑143 (MQ=255)                                                                                       gAAGATATGCGTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                               <  1:646745/143‑1 (MQ=255)                                                                                       gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                               >  1:644103/1‑143 (MQ=255)                                                                                        aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                              >  1:258764/1‑143 (MQ=255)                                                                                                  ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                    >  1:474957/1‑143 (MQ=255)                                                                                                  ttACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGCCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAc                                    >  1:377000/1‑143 (MQ=255)                                                                                                      gATGAAAACGGTACGCCAGTAGACATCGTACAGAACCAGCTGGGCGTACCGGCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                <  1:120678/143‑1 (MQ=255)                                                                                                        tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTTCGAAAGGTATCGGCGACAAGATCAACCCCCTGCTGAAACagcagc                              >  2:211379/1‑143 (MQ=255)                                                                                                        tGAAAACGGTACGCCGGTAGACATCGTACTGAAAAAGCTGGGCGTACCGTGTCGTATGAACATAGGTCAGATCCTCGAAAGCAAACTGGCTATGGCTGCGAAAGGTAACGGAGAGAAGATCAACGAAATGCTGAAAcagcaga                              >  2:218409/1‑142 (MQ=255)                                                                                                         gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAgcagca                             >  2:269692/1‑143 (MQ=255)                                                                                                              cGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCacgaag                        >  1:130601/1‑143 (MQ=255)                                                                                                                gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCACGAAACCGACCTGGGTATGGCTGCGAAAGGTATCGGGGACCAGATCAACGCCGTGCTGAAGCAGCAGCaagcacac                      >  2:512295/1‑138 (MQ=255)                                                                                                                    gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACACGCTCACCGCCACGCTGCACCAGCCGCAAGCCGTcgcgc                  >  1:252502/1‑142 (MQ=255)                                                                                                                      cGGTAGACATCGTACTGAACTAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                >  2:623338/1‑143 (MQ=255)                                                                                                                      cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCACGAGGTCGCGaaa                >  1:40363/1‑143 (MQ=255)                                                                                                                        gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCGAGAAGTCGCGAAAcc              >  2:675541/1‑142 (MQ=255)                                                                                                                            aCAGCGTACTGGACCAGGTGGGCGTACCGTCTCGTATGAACATCGGTCAGATGCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc          <  2:252502/143‑1 (MQ=255)                                                                                                                               tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa       >  1:552249/1‑143 (MQ=255)                                                                                                                                   aCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCCAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAACCTGCGCGGATTcc   >  1:295192/1‑142 (MQ=255)                                                                                                                                    cTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGACGTCGCGCAACTGCGCGACTTCAt  >  1:38882/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                        TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAT  >  NC_000913/4184405‑4184678

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