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

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%	158.3 / NA	55	P1100Q (CCG→CAG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/30);  total (25/30)

TAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184412‑4184682                                                                                                                                    |                                                                                                                                            tcaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCAGGGTAACAAGGGTGTCATTTCTAAGAGCATCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAc                                                                                                                                  <  2:658887/141‑1 (MQ=255)  aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAcc                                                                                                                                 >  1:113670/1‑143 (MQ=255)  aaaCGCCGTATCCAGCCTGGAGACAAGATGGCAGGTAGTCACGGTAACAAGGATGTAAACAACAGGTTCTAGCCGAGCGAAGATATGCCTTCCGATGAAAACGCTACGCCGGTAGACATCGACCTGCACCAGCTGGGCggccc                                                                                                                                 >  2:394645/1‑139 (MQ=255)    aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACACGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGCTAGACATCGGACTGAACCAGCTGGGCGCAccct                                                                                                                               >  2:261088/1‑141 (MQ=255)       ccGTCTCGAGTCTGGTGACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGAACTGtctc                                                                                                                            <  2:661050/143‑1 (MQ=255)       ccGTATGCAGCCGTGTGACAAGATGGCAGGTCGTCGCGGTAGCGTGGGTGTAATTTCTAAGATCAACCCGATAGACGATATGCCTTACGATGAAGACGGTACGCCGGTAGTTTTCGTGCTGAACCAGCTGGGCGTACCGtctc                                                                                                                            <  2:304817/143‑1 (MQ=255)       ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtctc                                                                                                                            >  1:576114/1‑143 (MQ=255)         gTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGt                                                                                                                          <  2:669387/143‑1 (MQ=255)          tATCCAGCCTGGTGACATGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTa                                                                                                                         <  2:517070/143‑1 (MQ=255)           aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGCAACCGCGCCGGCGGTATTCATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                        <  1:188558/143‑1 (MQ=255)           aTCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                        <  1:77424/143‑1 (MQ=255)                   tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCACATCg                                                                                                                >  1:882163/1‑143 (MQ=255)                   tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCg                                                                                                                >  1:922051/1‑143 (MQ=255)                       gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAGCATCGGGCa                                                                                                            >  1:593235/1‑143 (MQ=255)                        aCAAGATGGCAGGTCGTCACGGTAACAAGGATGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAg                                                                                                           >  1:554335/1‑143 (MQ=255)                         cAAGATGGCAGGCCGTCACGGCAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                          <  1:753494/143‑1 (MQ=255)                          gaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGCAAATCGGTCAGAt                                                                                                         >  2:867667/2‑143 (MQ=255)                          aaGATGGCAGGTCGTCACGGTAACAAGGGTGTACACTTAGTTTTTAATCCGATCGAAGATATGCCCTACGGGTATATCGGTTTGCCGGGAGACATCGTACGGAACCAGCTGGGCGGGGCGTCGCGTGTTAAGAGCGGTCGGAt                                                                                                         >  2:599151/1‑143 (MQ=255)                          aaGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGAt                                                                                                         >  2:415000/1‑143 (MQ=255)                            gctGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCATACCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGATCCAGCTGGGCGTACTGTCTCGTATGAACATCGGTCAGATcc                                                                                                       <  2:778702/141‑1 (MQ=255)                              tGGCGGGTCGTCACGGTAACAAGGGTGTAAGTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTc                                                                                                     <  1:823735/143‑1 (MQ=255)                                gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                   >  1:620688/1‑143 (MQ=255)                                gCAGGTCGCCACCGTAACAAACGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                   <  1:273235/143‑1 (MQ=255)                                 cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                                                                  >  2:231309/1‑143 (MQ=255)                                  cggTCGTCACGGTAACAAGGGTGTAATGTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                 <  2:143463/142‑1 (MQ=255)                                  aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATGAACCCGATCGAAGAGATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                                 <  2:768005/143‑1 (MQ=255)                                    gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGAGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAcc                                                                                               <  2:866897/143‑1 (MQ=255)                                     tcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCCGCGCCGCGAAAccc                                                                                              >  2:127631/1‑143 (MQ=255)                                       gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTGCGCCGGTAGACATCGTACTGGACCAGCTGGCCGTACCGTCTCGTATAAACATCGGGCAGATCGTCGAAAgccgc                                                                                            >  1:236493/1‑141 (MQ=255)                                          aCGGTAACAAGTGTGTAATTTCTAAGCTCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGCACCGCCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                         <  2:624708/143‑1 (MQ=255)                                          aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAACTg                                                                                         >  1:294501/1‑143 (MQ=255)                                                  aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACCTCGGTCAGATCCTCGAAACCCACCCGGGTATGGc                                                                                 >  2:46227/1‑143 (MQ=255)                                                     ggTGTACTTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                              <  2:481798/143‑1 (MQ=255)                                                         tattttCTAAGTTCAACCCGATCGCAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                          <  1:40168/140‑1 (MQ=255)                                                         tagcttCTAAGATCAACCCGATCGAAGATATGCCCGCCGATGAAAACGGCACGCCGGTCGCCATCGTACTGAACCAGCCGGGCGCACCGTCGCGTATGAACATCGGTCAGCTCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                          <  2:515638/139‑1 (MQ=255)                                                         tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                          <  1:86434/143‑1 (MQ=255)                                                         tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTCGGTATGGCTGCGaaa                                                                          >  2:781746/1‑143 (MQ=255)                                                               gcacgATCAATCCGATCGAAGATATGCCGTACGATGAAGACGGTACGCGGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGAACATCGGTCAGAGCCTCGAATCCCACCTTGGTATGGCTGCGAAAGGTATc                                                                    <  2:371545/139‑1 (MQ=255)                                                               cTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGGACTGAACCAGCTGGGCGTACCGTCACGTATGAACATCGGACAGATCCTCGAAACCCACCGGGGGATGGCTGCGAAAgggagc                                                                    >  2:191265/1‑139 (MQ=255)                                                                tAAGATCAAGCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTTGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                   <  2:392981/143‑1 (MQ=255)                                                                         cccGATCGACGAGATCGCTGGTGATGCAAACGGGACGCCGGTTGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATAGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                          <  2:469735/143‑1 (MQ=255)                                                                         cccGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                          <  1:127331/143‑1 (MQ=255)                                                                          ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAGACCCACATGGGTATGGCGGAGAAAGGTATCGGCGACCAGAt                                                         >  2:705739/1‑143 (MQ=255)                                                                          ccGAGCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAt                                                         <  2:761063/143‑1 (MQ=255)                                                                                gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATACTCGAAACCCACCTGGGTATGGCTGCGACACGGATCGACGACAAGATCAACGc                                                   >  2:871459/1‑143 (MQ=255)                                                                                                aTGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagcag                                   >  1:543954/1‑143 (MQ=255)                                                                                                 tGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                  <  1:645847/143‑1 (MQ=255)                                                                                                      aCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                             <  1:599151/143‑1 (MQ=255)                                                                                                             gCCGTTAGTCACCGTACTGAACCATCTCTGCTTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                      <  2:328231/143‑1 (MQ=255)                                                                                                                     aCATCGTACTGAACCAGCTGGGCGTACCGACTCGTATGAACACCGGTCAGAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc              <  1:521575/143‑1 (MQ=255)                                                                                                                      cATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGcgcg             >  1:575311/1‑143 (MQ=255)                                                                                                                       aTCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGa            >  2:975440/1‑143 (MQ=255)                                                                                                                         cGTACTGAGCCAGCTGGGCGTACAGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAAt          <  2:781108/143‑1 (MQ=255)                                                                                                                                 aCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAg  >  1:568247/1‑143 (MQ=255)                                                                                                                                 aCCAGCTGGGCGTACCGTCTCGTATGAACAGCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAg  <  1:249717/143‑1 (MQ=255)                                                                                                                                    |                                                                                                                                            TAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCAG  >  NC_000913/4184412‑4184682

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