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

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184405‑4184679                                                                                                                                           |                                                                                                                                         tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGGTGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                      >  1:365114/1‑143 (MQ=255)  ggCGGTTAAACGCCGTATCCAGCCTGGGGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                     <  1:452218/143‑1 (MQ=255)    cGGTTAACCGCCGTATCCAGCCTGGTGTCAAGATGGCAGGTCGTCACGGTAACAAGAGTGTAATTTCTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGc                                                                                                                                   <  2:271971/143‑1 (MQ=255)           aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGCTAACAAGGGTGTAATTACTAAGATCAACCCGACCGAAGATATGCCTTACGCTGACAACGGCACGCCGGCAGACACCGTACTGCACCAGATGGACATAACGt                                                                                                                            >  2:299170/1‑143 (MQ=255)                  acccAGCCTGCTGACAACATTGCAGGTCGTCACGGTAACAAGGGTGTAATTCCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                     <  1:23874/141‑1 (MQ=255)                        ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGTTGAGAACTTTGCGTGGGTAGACACCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAt                                                                                                               >  1:289122/1‑143 (MQ=255)                          tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCAGTAGAGATCGTACTGCACCAGCTGGGCGGACCGCGACGTATGAACGTCg                                                                                                             >  2:270580/1‑143 (MQ=255)                            gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGt                                                                                                           >  1:243066/1‑143 (MQ=255)                             tGACAAGCTGGCAGGTCGTCACGGTAAGAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGCGTACTGTCTCGTATGAACATCGGTc                                                                                                          <  2:29480/143‑1 (MQ=255)                             gggtaggATCGCGGGTCGTGACGGGAACAAGGGTGTCAGTTCTAACATCAACCCTATCGAAGATATGCCGTACGATGAAAACGTGACGCGGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                          <  2:364766/137‑1 (MQ=255)                                cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGa                                                                                                       >  1:674114/1‑143 (MQ=255)                                   gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACAACGGTCAGATcc                                                                                                    >  2:400251/1‑143 (MQ=255)                                        cAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaa                                                                                               <  1:490018/143‑1 (MQ=255)                                         aGGTCGTCACGGTAACAAGGGTGTAATCACTAAGATCATCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                              <  1:324161/143‑1 (MQ=255)                                          ggTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAAc                                                                                             <  1:39743/143‑1 (MQ=255)                                                 aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCAcccg                                                                                      >  2:648632/1‑141 (MQ=255)                                                  cGGTAACAAGGGTGTAATTTCTAAGATCAACACGATAGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgg                                                                                     <  1:26130/143‑1 (MQ=255)                                                   ggTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTggg                                                                                    <  1:61573/143‑1 (MQ=255)                                                        caaagttgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATgg                                                                               <  2:29464/137‑1 (MQ=255)                                                         aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGCTCCTCGAAACCCACCTGGGtctcgc                                                                              >  2:66809/1‑138 (MQ=255)                                                              tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTCTGGCTGCgc                                                                         >  2:670054/1‑142 (MQ=255)                                                                     tCTAAGATCAACACGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAt                                                                  >  1:706227/1‑143 (MQ=255)                                                                        aaGATCAACCCGATCGAAGAGATGTCTTAGGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCgg                                                               <  2:652954/143‑1 (MQ=255)                                                                            tCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGGAACGGCGAc                                                           >  2:603061/1‑143 (MQ=255)                                                                                   gATCGAAGATATGCCTTGCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                                                    <  1:405235/143‑1 (MQ=255)                                                                                   gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                                                    <  1:405232/143‑1 (MQ=255)                                                                                   gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCa                                                    <  1:405227/143‑1 (MQ=255)                                                                                       gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTTGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGc                                                >  1:96798/1‑143 (MQ=255)                                                                                       gAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGACGTACCGTCTCGTATGAACCTCGGTCAGACCCTCGACACCCCCCTGGGTATCTCTGCGAAAGGAAAGAGCAAAAAGATCAAcac                                                >  1:192284/1‑141 (MQ=255)                                                                                                   gaCGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAca                                    <  1:644920/142‑1 (MQ=255)                                                                                                      gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCCGCGACAAGATCAACGCCATGCTGAAAcagcc                                 >  1:536063/1‑142 (MQ=255)                                                                                                        tgacgaCGGCACGCTGGTAGACACCGTTCTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGACCCTCCACTCCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                               <  2:297474/138‑1 (MQ=255)                                                                                                                          aGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgc             <  2:622210/143‑1 (MQ=255)                                                                                                                           gACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcg            >  1:97388/1‑143 (MQ=255)                                                                                                                                     tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCACCGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATGCATc  >  2:741745/1‑143 (MQ=255)                                                                                                                                     tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAATCCGCAGCACGGAGTCGCGACACTGCGCGacgtcgtc  >  1:207611/1‑136 (MQ=255)                                                                                                                                           |                                                                                                                                         TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184405‑4184679

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