Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,890	A→G	D549G (GAC→GGC)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,890	0	A	G	100.0%	169.1 / NA	51	D549G (GAC→GGC)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/25);  total (26/25)

GCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  NC_000913/4182742‑4183030                                                                                                                                                     |                                                                                                                                             gcagcaGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCg                                                                                                                                            >  2:590920/1‑151 (MQ=255)  cagcagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGt                                                                                                                                           <  2:765187/151‑1 (MQ=255)         gAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCg                                                                                                                                    >  2:660301/1‑151 (MQ=255)                ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCAc                                                                                                                              >  1:294516/1‑150 (MQ=255)                 tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCGCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACt                                                                                                                             >  2:275930/1‑150 (MQ=255)                 tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACt                                                                                                                             <  2:774564/150‑1 (MQ=255)                 tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACt                                                                                                                             <  2:511408/150‑1 (MQ=255)                    tcGTTTCCAGCCAGTTGTCTCAGTTTATGGACCAGAACAACCCGTTGTCTGAGATTACGCACAACCGTCGTATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTAc                                                                                                                           <  1:827388/149‑1 (MQ=255)                     cGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGt                                                                                                                        >  2:216376/1‑151 (MQ=255)                        ttCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTcgc                                                                                                                     >  1:393764/1‑151 (MQ=255)                         tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGGCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTcgcg                                                                                                                    >  1:349856/1‑151 (MQ=255)                           cagccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGt                                                                                                                   >  1:809554/1‑150 (MQ=255)                             gccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATg                                                                                                                <  2:45609/151‑1 (MQ=255)                                   tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAAt                                                                                                          >  2:607953/1‑151 (MQ=255)                                        cAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa                                                                                                      <  2:12249/150‑1 (MQ=255)                                            ttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAcccc                                                                                                 >  1:401837/1‑151 (MQ=255)                                              aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                                <  2:515360/150‑1 (MQ=255)                                               tGGACCAGAACAACCCGCTGTCTGAGATTACGCCCAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGTCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                                                                                <  1:572334/149‑1 (MQ=255)                                                ggACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                                                                                             <  1:221713/151‑1 (MQ=255)                                                        acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGa                                                                                      >  1:328182/1‑150 (MQ=255)                                                          aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACa                                                                                   <  2:333691/151‑1 (MQ=255)                                                            cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc                                                                                 >  2:685155/1‑151 (MQ=255)                                                               gCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGt                                                                              >  1:429400/1‑151 (MQ=255)                                                                      gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATc                                                                        <  2:699629/150‑1 (MQ=255)                                                                                acaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCGACTCTCTGTc                                                             >  2:36524/1‑151 (MQ=255)                                                                                 caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc                                                            >  2:794199/1‑151 (MQ=255)                                                                                    aCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtg                                                         <  2:809610/151‑1 (MQ=255)                                                                                      gtcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCATGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                                        >  1:68300/1‑150 (MQ=255)                                                                                        cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                                     <  1:61506/151‑1 (MQ=255)                                                                                        cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                                     <  2:412842/151‑1 (MQ=255)                                                                                         gtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                    >  1:803100/1‑151 (MQ=255)                                                                                          tATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                                    >  2:220929/1‑150 (MQ=255)                                                                                           aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                  >  1:475532/1‑151 (MQ=255)                                                                                               ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                                                  <  2:39718/147‑1 (MQ=255)                                                                                                cGCACTCGGCCCAGGCGGTCTGACCCGTGGACGTGCAGGCTTCGAAGGTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAAATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGa                                               >  1:832523/1‑149 (MQ=255)                                                                                                cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                             <  2:840419/151‑1 (MQ=255)                                                                                                cGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                             <  1:386177/151‑1 (MQ=255)                                                                                                 gCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACt                                             >  1:392491/1‑150 (MQ=255)                                                                                                      cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACg                                         <  1:270003/149‑1 (MQ=255)                                                                                                      cGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        <  2:827008/150‑1 (MQ=255)                                                                                                         cccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                    >  2:609611/1‑151 (MQ=255)                                                                                                             ggCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCt                                <  1:383484/151‑1 (MQ=255)                                                                                                                ggTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                             <  2:824924/151‑1 (MQ=255)                                                                                                                    tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          >  2:650422/1‑150 (MQ=255)                                                                                                                     gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          <  2:142354/149‑1 (MQ=255)                                                                                                                        ccGTGAACGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt                     >  2:823012/1‑151 (MQ=255)                                                                                                                              aCGTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTAt               >  1:352042/1‑151 (MQ=255)                                                                                                                                gTGCAGGCTTCGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg             >  2:49121/1‑151 (MQ=255)                                                                                                                                          cGAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc   <  1:208297/151‑1 (MQ=255)                                                                                                                                           gAAGTTCGAGGCGTACGCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGAcc   <  1:279697/150‑1 (MQ=255)                                                                                                                                           gAAGTTCGAGGCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  <  2:176654/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             GCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  NC_000913/4182742‑4183030

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