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

GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184409‑4184680                                                                                                                                       |                                                                                                                                          ggTTAAACGCCGTATCCAGCCTGATGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCAAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGCAGACAGCGCACTGAACCAGCTGGGCg                                                                                                                                   >  2:342123/1‑143 (MQ=255)               tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATg                                                                                                                     <  2:403321/143‑1 (MQ=255)                 cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGACATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGaa                                                                                                                   >  2:358628/1‑143 (MQ=255)                         tGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTc                                                                                                           >  1:554358/1‑143 (MQ=255)                          gACAAGATGGCAGGTCGGCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                          <  2:51587/143‑1 (MQ=255)                                   gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAGCCAGCTGGGCGTACCGTCTCGCATGAACATGGGTCAGATCCTCGa                                                                                                 >  1:552850/1‑143 (MQ=255)                                   gCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGa                                                                                                 >  1:624199/1‑143 (MQ=255)                                     aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGaaa                                                                                               <  1:176253/143‑1 (MQ=255)                                      ggTCGTCACGGTAACAAGGGTGTAAGTTCTAAGACCAACCCGAGCGAAGATATGCCTTACGATGAAAACGGTACTCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGCCTCGCAAGCACATCGGTCAGATCCTCGAAAc                                                                                              <  2:402605/143‑1 (MQ=255)                                             aCGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTg                                                                                       >  1:567745/1‑143 (MQ=255)                                                  aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGAACCAGCTGGGCGTACCCTCGAGTACAACCATCGGGCCGAGCATCGAAGCCCGGCTGGGtgt                                                                                  >  2:422115/1‑141 (MQ=255)                                                  aaCAAGGGCGTAATTTCTAAGATCAACCCGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCCCGTATGACCATCGGCCAGATCCCCGACACCCACCTGGGCAt                                                                                  >  2:97078/1‑143 (MQ=255)                                                   aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGCTATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATg                                                                                 >  1:425958/1‑143 (MQ=255)                                                   aCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTgggggtg                                                                                 >  1:356702/1‑139 (MQ=255)                                                        tgtgtAATTTCTAAGATCAACCCGATCGAAGATATCCCTTACGATGAAAACGGTACGCCGGTAGACTTCGTACTGATCCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGc                                                                            <  2:86181/142‑1 (MQ=255)                                                           gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                         <  1:488518/143‑1 (MQ=255)                                                            tAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAACAGCTGGGCGTACCGTCACGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGGATGGCTGCGaaa                                                                        >  2:623844/1‑143 (MQ=255)                                                                ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                                                                    >  1:483135/1‑143 (MQ=255)                                                                  ccaaGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCTGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATc                                                                  <  1:14903/141‑1 (MQ=255)                                                                   tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGACGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                 >  2:245446/1‑143 (MQ=255)                                                                      gcggaaCCGGATCGAAGAGATGCCTTACGATGAGCGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGTGCGTACCGTCTGGTATGAACATCGGTCAGATCGTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCg                                                              <  2:501224/139‑1 (MQ=255)                                                                         cAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAGAGGTCTCGGCGACa                                                           >  1:137804/1‑143 (MQ=255)                                                                          aaCCCGATCGAAGATATGCCTGACGCTGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                          <  2:617835/143‑1 (MQ=255)                                                                           aCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGGAAGGTGTCgggggcggg                                                         >  1:132778/1‑136 (MQ=255)                                                                             ccGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCACGTATGAACATCGGTCAGATCCTCGAAGCCGACCTGGGTATCGCTGCGAAAGGTATCGGCGACAAGac                                                       >  2:167004/1‑142 (MQ=255)                                                                                     ggATATGCCTTACGATGAAAACGGTACGCCGGTAGACACCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCa                                               <  1:131177/142‑1 (MQ=255)                                                                                                aCGATGAAAACGGTCCGCCGAGAGACATCGCACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGCGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                    <  1:320474/143‑1 (MQ=255)                                                                                                 cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                   >  1:426374/1‑143 (MQ=255)                                                                                                  gacgAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                  <  1:130988/140‑1 (MQ=255)                                                                                                  gATGAAAACGGTAGGCCGGGAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACGTCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                  <  1:484047/143‑1 (MQ=255)                                                                                                     gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGCTCAACGCCATGCTGAAACAgcagca                               >  1:81666/1‑143 (MQ=255)                                                                                                     gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCAACAAGATCAACGCCATGCTGAAACAgaagca                               >  2:276890/1‑143 (MQ=255)                                                                                                     gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCAGGGGATGGCTGCGAAAGGTATCGGCGACAAGAAGAACGAGATGAGCAAACAgcagca                               >  2:133680/1‑143 (MQ=255)                                                                                                       aaaCGGTGCGACGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                             <  1:50403/143‑1 (MQ=255)                                                                                                        aaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCCCGCAACCCACCTGGGTATGGCCGCGAAAGGTCTCCGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                            >  2:583285/1‑143 (MQ=255)                                                                                                                                  gAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTAACGGCGACAAGATCAACGCCAAGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATcc  >  1:557943/1‑143 (MQ=255)                                                                                                                                       |                                                                                                                                          GGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCC  >  NC_000913/4184409‑4184680

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