Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,183,814	T→A	V857E (GTG→GAG)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,183,814	0	T	A	100.0%	201.7 / NA	52	V857E (GTG→GAG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (30/22);  total (30/22)

GGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTT  >  NC_000913/4183669‑4183961                                                                                                                                                  |                                                                                                                                                    ggTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGt                                                                                                                                                 <  2:254964/150‑1 (MQ=255)        aCTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTgc                                                                                                                                         >  2:266058/1‑151 (MQ=255)          ttCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTgcgc                                                                                                                                       >  1:760763/1‑151 (MQ=255)             gcaGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGATTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGATGGGTCCGGAAGAGATCACCGCTGACATCACGAACGAGGGTGAAGCTGCGccct                                                                                                                                    >  2:521214/3‑151 (MQ=255)                gACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTcc                                                                                                                                  >  2:419384/1‑150 (MQ=255)                   tCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAAc                                                                                                                              >  1:850643/1‑151 (MQ=255)                    ccATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACt                                                                                                                             <  2:550703/151‑1 (MQ=255)                     cATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTACAAAATg                                                                                                                            >  2:305878/1‑151 (MQ=255)                            gTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAAt                                                                                                                     >  1:186362/1‑151 (MQ=255)                            gTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAAt                                                                                                                     >  1:26495/1‑151 (MQ=255)                             tATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATc                                                                                                                    >  2:763382/1‑151 (MQ=255)                                ccGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCgg                                                                                                                 >  1:260002/1‑151 (MQ=255)                                   aGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTAt                                                                                                              >  2:199372/1‑151 (MQ=255)                                           gttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTa                                                                                                        <  1:334342/149‑1 (MQ=255)                                           gttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTAc                                                                                                       <  1:514117/150‑1 (MQ=255)                                            ttCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACAt                                                                                                     >  2:501157/1‑151 (MQ=255)                                                   aGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCg                                                                                              >  1:428867/1‑151 (MQ=255)                                                        gTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGt                                                                                         >  2:224133/1‑151 (MQ=255)                                                              ccaccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCgg                                                                                   <  1:92670/151‑1 (MQ=255)                                                               caccaTCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGt                                                                                  <  2:450625/151‑1 (MQ=255)                                                                 ccaTCCACATTCAGGAACTGGCGTGTGTGTTCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTgg                                                                                >  2:516675/1‑151 (MQ=255)                                                                     ccACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGAc                                                                            >  1:214748/1‑151 (MQ=255)                                                                      cacaTTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACa                                                                           <  2:192204/151‑1 (MQ=255)                                                                        caTTCAGGATCTTGCGTGTGGGTACCGTGACACCAAGCTGGCTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGGGAAGCTGCGCTCTCCAAACTGGATGAATCTGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGGGAAAtt                                                                         <  2:152000/151‑1 (MQ=255)                                                                            cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTTCGGAAGTGACCGGTGGCGACATTCt                                                                       <  2:409058/149‑1 (MQ=255)                                                                            cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg                                                                     >  1:633970/1‑151 (MQ=255)                                                                            cAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCGCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtgg                                                                     <  2:751960/151‑1 (MQ=255)                                                                               gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttg                                                                  <  1:23888/151‑1 (MQ=255)                                                                               gAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATGCtggttg                                                                  <  2:449456/151‑1 (MQ=255)                                                                                 aCTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggttggt                                                                <  2:63061/151‑1 (MQ=255)                                                                                  cTGGCGTGTGTGTCCCGTGAAAACAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGACTCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCtggctggt                                                                >  2:212618/1‑150 (MQ=255)                                                                                    ggCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggta                                                               >  2:730221/1‑149 (MQ=255)                                                                                         gtgtgtCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTggtaaggtaa                                                         <  2:457290/150‑1 (MQ=255)                                                                                               cccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGa                                                   >  2:109907/1‑150 (MQ=255)                                                                                                ccGTGACACTAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaaa                                                 <  1:100999/151‑1 (MQ=255)                                                                                                ccGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGaaa                                                 <  1:514391/151‑1 (MQ=255)                                                                                                  gTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAgg                                               >  2:773338/1‑151 (MQ=255)                                                                                                    gACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTg                                             <  1:564977/151‑1 (MQ=255)                                                                                                        ccAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCAGAAAGGTGAAAc                                         >  2:464054/1‑151 (MQ=255)                                                                                                         cAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACt                                        <  1:485646/151‑1 (MQ=255)                                                                                                               gggTCCGGAAGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTg                                  >  1:321007/1‑151 (MQ=255)                                                                                                                     ggAAGAGCTCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACCCCGAAAGGCGAAACTCAGCTGAccc                              >  2:164358/1‑149 (MQ=255)                                                                                                                      gAAGAGATCACCGCTGACATACCGAACGAGGGAGAAGCTGCGCTGTCCAAACTGGATGAATCCGGTATCGTGTACATTGGAGAGGAAGTGACCGGGGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGGGAAAATCAGCTGa                                 >  2:152578/1‑145 (MQ=255)                                                                                                                      gAAGAGACCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTACGGTAACGCCGAAAGGTGAAACTCAGCTGAc                                >  2:625539/1‑146 (MQ=255)                                                                                                                       aaGAGATCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGCTGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGAcccc                             >  2:276425/1‑148 (MQ=255)                                                                                                                         gagaTCACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCagaa                         <  2:609437/150‑1 (MQ=255)                                                                                                                              cACCGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGaaaaa                   <  1:426293/151‑1 (MQ=255)                                                                                                                                ccGCTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGCTGAATCCGGTATCGTTTAAATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAAct                 >  2:710940/1‑151 (MQ=255)                                                                                                                                   cTGACATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAActgct              >  1:830112/1‑151 (MQ=255)                                                                                                                                      aCATCCCGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCtgcg           <  1:305077/151‑1 (MQ=255)                                                                                                                                           ccGAACGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCtgcgtgcg       >  2:712776/1‑150 (MQ=255)                                                                                                                                               aCGAGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCtt  <  2:259941/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GGTTACAACTTCGAAGACTCCATCCTCGTATCCGAGCGTGTTGTTCAGGAAGACCGTTTCACCACCATCCACATTCAGGAACTGGCGTGTGTGTCCCGTGACACCAAGCTGGGTCCGGAAGAGATCACCGCTGACATCCCGAACGTGGGTGAAGCTGCGCTCTCCAAACTGGATGAATCCGGTATCGTTTACATTGGTGCGGAAGTGACCGGTGGCGACATTCTGGTTGGTAAGGTAACGCCGAAAGGTGAAACTCAGCTGACCCCAGAAGAAAAACTGCTGCGTGCGATCTT  >  NC_000913/4183669‑4183961

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