Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,803	C→T	P520L (CCG→CTG)	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,803	0	C	T	100.0%	160.1 / NA	44	P520L (CCG→CTG)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/19);  total (25/19)

AGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCG  >  NC_000913/4182659‑4182945                                                                                                                                                 |                                                                                                                                               agCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACctgctg                                                                                                                                            <  2:542518/149‑1 (MQ=255)     tGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTAAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTgag                                                                                                                                      <  2:297375/151‑1 (MQ=255)         gTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTa                                                                                                                                  >  2:147661/1‑151 (MQ=255)             aaGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGca                                                                                                                              <  2:456370/151‑1 (MQ=255)                          tctGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTAtc                                                                                                                 <  2:434705/151‑1 (MQ=255)                          tctGGGCGACCTGGATGCCCTGATGCCACAGGATATGATCAACGCCACGCCGATTTCCGCAGCAGTGAAAGAGTGCTTCGGTTCCAGCCACCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCCCACACGTCGTAtc                                                                                                                 <  1:426117/151‑1 (MQ=255)                                cGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGc                                                                                                            <  2:19882/150‑1 (MQ=255)                                    cTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGCTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCg                                                                                                       >  2:304754/1‑151 (MQ=255)                                      ggATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCgg                                                                                                      >  2:267694/1‑150 (MQ=255)                                        aTACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGccc                                                                                                   <  1:583290/151‑1 (MQ=255)                                           cccTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAgg                                                                                                <  2:346951/151‑1 (MQ=255)                                                aTGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt                                                                                            >  2:657039/1‑150 (MQ=255)                                                        ggATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                                                                                   >  2:259960/1‑151 (MQ=255)                                                        ggATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                                                                                   >  2:31496/1‑151 (MQ=255)                                                        ggATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                                                                                   >  1:283166/1‑151 (MQ=255)                                                        ggATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTACGCACTCGGCCCAGGCGGTCTGACCCGt                                                                                   >  2:236036/1‑151 (MQ=255)                                                         gATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTg                                                                                  >  1:763677/1‑151 (MQ=255)                                                            atGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAAc                                                                               <  1:350871/151‑1 (MQ=255)                                                              gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGt                                                                             >  1:301376/1‑151 (MQ=255)                                                              gATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGCGACCGt                                                                             >  1:98628/1‑151 (MQ=255)                                                                   aCGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAgg                                                                        >  2:736385/1‑151 (MQ=255)                                                                           ccGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa                                                                 <  2:735116/150‑1 (MQ=255)                                                                                        cagTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTaca                                                   >  2:93531/1‑151 (MQ=255)                                                                                          gTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                                  >  1:221410/1‑150 (MQ=255)                                                                                           tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                                                <  2:321978/151‑1 (MQ=255)                                                                                           tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                                 <  1:298633/150‑1 (MQ=255)                                                                                           tGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                                 <  2:677210/150‑1 (MQ=255)                                                                                                     cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTAc                                      <  2:42895/151‑1 (MQ=255)                                                                                                     cttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTAc                                      <  2:232168/151‑1 (MQ=255)                                                                                                        cGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGt                                   <  1:726962/151‑1 (MQ=255)                                                                                                          gTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg                                 >  2:721293/1‑151 (MQ=255)                                                                                                             ccagccagCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                              >  1:287474/1‑151 (MQ=255)                                                                                                             ccagccagCTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGt                              >  2:373244/1‑151 (MQ=255)                                                                                                                 ccagcTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAt                            <  2:758590/149‑1 (MQ=255)                                                                                                                     cTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc                        <  1:657032/149‑1 (MQ=255)                                                                                                                     cTGTCTCAGTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc                        <  2:408160/149‑1 (MQ=255)                                                                                                                             gTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAc               >  1:392158/1‑150 (MQ=255)                                                                                                                             gTTTATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAc               >  2:865465/1‑150 (MQ=255)                                                                                                                              tttATGGACCAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAAccc             <  2:782595/151‑1 (MQ=255)                                                                                                                                      ccAGAACAACCTGCTGTCTGAGATTACGCACAGACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt     >  2:348316/1‑151 (MQ=255)                                                                                                                                      ccAGAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGt     >  2:745428/1‑151 (MQ=255)                                                                                                                                         gAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg  >  2:439502/1‑151 (MQ=255)                                                                                                                                         gAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg  >  1:337668/1‑151 (MQ=255)                                                                                                                                         gAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg  >  1:726037/1‑151 (MQ=255)                                                                                                                                         gAACAACCTGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg  >  1:673480/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                               AGCGTGCGGTGAAAGAGCGTCTGTCTCTGGGCGATCTGGATACCCTGATGCCACAGGATATGATCAACGCCAAGCCGATTTCCGCAGCAGTGAAAGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCG  >  NC_000913/4182659‑4182945

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