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

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