Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,184,543	C→A	P1100Q (CCG→CAG)	rpoB →	RNA polymerase, beta subunit

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

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184405‑4184679                                                                                                                                           |                                                                                                                                          tGGCGGTTAAACGCCGCCTCCAGCCTGGTGACAAGATGGCAGGTCGTCACTGCAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTg                                                                                                                                       <  1:676658/143‑1 (MQ=255)  gggggTGAACCGCCGGATCCAGCCTGGTGACGAGATGGCACGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGAGGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                      <  2:696851/140‑1 (MQ=255)  ggCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTgg                                                                                                                                      <  1:812067/143‑1 (MQ=255)     ggTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCg                                                                                                                                   <  1:732892/143‑1 (MQ=255)         aaaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTAcc                                                                                                                               >  1:84983/1‑143 (MQ=255)         aaaCGCCGTATCCAGCCTGGGGACAAGATAGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGCTCAACCCGATCGAAGATGTGCCTTACGATGAAAACGGTACGCCGGTCGACATCGTACTGAACCAGCGGGGCGTAcc                                                                                                                               >  1:445835/1‑143 (MQ=255)             gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGGGTACCGtct                                                                                                                           >  2:259230/1‑143 (MQ=255)                 tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGTCGTACCGTCTCGTa                                                                                                                       <  1:474239/143‑1 (MQ=255)                        ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTCTGAACAt                                                                                                                >  1:325484/1‑143 (MQ=255)                         cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATc                                                                                                               >  1:860942/1‑143 (MQ=255)                              gACATGATGGCAGGGCGTCCCGGTAACGAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCa                                                                                                          <  2:853038/143‑1 (MQ=255)                                 aggATGGGAGGTCGTCACAGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGACGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GAt                                                                                                      <  2:616499/141‑1 (MQ=255)                                    aTGGCAGGTCGCCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCt                                                                                                   <  2:149970/143‑1 (MQ=255)                                     tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTc                                                                                                  <  2:859228/143‑1 (MQ=255)                                      ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTCAGATCAACCCGATCGAAGATAAGCGTTACGACGAAAACGGTACGCCGGTAGACATCGTAATGAACCAGCTGGACGTCGCGTCTCGCATGAACCTCGGTCA‑ggaaacgc                                                                                                 >  2:832353/1‑136 (MQ=255)                                         aGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCaaaa                                                                                              >  2:589500/1‑143 (MQ=255)                                         aGGTCGTCACGGGCACACGGGTGTAATTTCTAACCTCAACTCGATCGCAGTTTTGCCTTACGATGAAAACGGTTCGCCGGTTGCCATCGTACTGAACCAGCTGGGCGTACCCTCTCGTATGAACATCGGTCA‑GATCCTCGaaa                                                                                              <  1:339020/143‑1 (MQ=255)                                                cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCt                                                                                       >  1:799413/1‑143 (MQ=255)                                                     tAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTa                                                                                  <  1:538201/143‑1 (MQ=255)                                                      aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTATAGTTACTACATCGGTCA‑GATCCTCGAAACCCACCTGGGTac                                                                                 >  2:196975/1‑142 (MQ=255)                                                              tgtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATACTCAATACGCACCCGGGCAAGACTGCGa                                                                         >  2:621317/1‑143 (MQ=255)                                                                taactccTAAGAGCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGaaa                                                                       <  2:703049/137‑1 (MQ=255)                                                                    ttCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTa                                                                   >  1:787678/1‑143 (MQ=255)                                                                    ttCTAAGATCAACCCGATCGAAGATATACCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGGCA‑GATCCTCGAAACCCACCCGGGGATGGCTGCGAAAGGtg                                                                   >  1:384945/1‑142 (MQ=255)                                                                       tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCg                                                                >  1:853687/1‑143 (MQ=255)                                                                       tAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACATGGGTATGGCTGCGAAAGGCAtgg                                                                >  1:626668/1‑141 (MQ=255)                                                                         aGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGc                                                              >  1:492707/1‑143 (MQ=255)                                                                              aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGCCaa                                                         >  1:403360/1‑143 (MQ=255)                                                                                cccGATCGCAGAGATGCCGTGCGGTGATGACGGTACGCCGGTAGGCAACTTACGGAACCAGCTGGGGGGACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGa                                                       <  2:423536/143‑1 (MQ=255)                                                                                  cGATCGAAGATACGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATc                                                     <  2:555609/143‑1 (MQ=255)                                                                                     tCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCAACTGGGTATGGATGCGAAAGGTATCGGCGCCAAGATCAAc                                                  >  1:140888/1‑143 (MQ=255)                                                                                        aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTCTGGCTGCGCAAGGTGTCGGCGACAAGATCAACGcc                                               >  1:609452/1‑143 (MQ=255)                                                                                        aaGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGcc                                               <  1:700525/143‑1 (MQ=255)                                                                                            tatGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGc                                           <  2:157303/143‑1 (MQ=255)                                                                                                   caCGCTGCAACCGGTACGCCGCTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCGCGAACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAca                                    <  2:138398/142‑1 (MQ=255)                                                                                                    aCGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTTGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGACCAAGCCCATGCCGAAAcag                                   >  2:786062/1‑143 (MQ=255)                                                                                                    aCGAGGAAAGCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcag                                   <  2:801992/143‑1 (MQ=255)                                                                                                     cGATGACAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GAGCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                  <  2:635309/143‑1 (MQ=255)                                                                                                      gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAACcagca                                 >  1:782831/1‑143 (MQ=255)                                                                                                      gATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagca                                 >  2:693874/1‑143 (MQ=255)                                                                                                           aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                            <  1:889177/143‑1 (MQ=255)                                                                                                            aaCTGTACGCCAGTAGAAATCGGACTGAACCAGCTGGGCGTGCCGGAGCGTATGAACATCGGTCA‑GATCCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaaga                           <  1:381339/143‑1 (MQ=255)                                                                                                                 tACGCCGGTAGACATCGCACTGAACCAGCGGGGCGGACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                      <  1:7637/143‑1 (MQ=255)                                                                                                                    gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACACCGCTAATGAT‑AACTAATCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                   <  1:441686/143‑1 (MQ=255)                                                                                                                     ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaa                  >  1:549813/1‑143 (MQ=255)                                                                                                                     ccGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCGGAAACAGCAGCAAGAGGTCGCGaa                  >  1:144645/1‑143 (MQ=255)                                                                                                                      cGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGaaa                 >  1:564924/1‑143 (MQ=255)                                                                                                                               tCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGaa        <  1:406207/143‑1 (MQ=255)                                                                                                                                  tACTGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTc     <  1:196975/143‑1 (MQ=255)                                                                                                                                   aCTGAACCAGCTGGGCGTCCCGCTGCGTATGAACATCGGTCA‑GATCCACGCAACCCACCTGGGGATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCa    <  2:33232/143‑1 (MQ=255)                                                                                                                                    cTGAACCAGCTGGGCGTACCGTCTAGTATGAACATCGGTGA‑GATCCTCGAAACGCACCTGGGTATGGCTGGGAAAGGTATCGGCGACAAGATGAACGCGATGCTGAAACAGCAGCAAGAAGGCGCGAAAATGCGCGAGTTCAt   >  1:578571/1‑143 (MQ=255)                                                                                                                                     tGAACCAGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAGGAAGTCGCGAAACTGCGCGAATTCATc  >  1:483947/1‑143 (MQ=255)                                                                                                                                     tGAACCAGCTGGGCGTACCGTCTCGTATGAAAATAGCAGC‑TACTCGCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGAACAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc  <  2:420396/143‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGAACATCGGTCA‑GATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATC  >  NC_000913/4184405‑4184679

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