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

TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184405‑4184681                                                                                                                                           |                                                                                                                                            tGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTCTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACCGAACCAGCTg                                                                                                                                         >  2:637847/1‑143 (MQ=255)    cGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGc                                                                                                                                      <  1:261113/143‑1 (MQ=255)          aaCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCg                                                                                                                                >  2:686040/1‑143 (MQ=255)           aCGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGt                                                                                                                               >  2:25815/1‑143 (MQ=255)            cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGtc                                                                                                                              <  2:663024/143‑1 (MQ=255)                  actcAGCCTGGTGAGAAGTTGGGAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCCTCCCGATCCAAGATATGCCTTACGATGAAAGCGGTACGCCGGTAGACACCGTACTGAACCAGCTGGGCGTACCGTCTCGTAt                                                                                                                        <  2:444835/140‑1 (MQ=255)                    ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGACCCAGCTGGGCGTACCGTCTCGTATGa                                                                                                                      >  1:2576/1‑143 (MQ=255)                     cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAATATATGCCTTACGATGACACCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGa‑a                                                                                                                    <  2:656443/143‑1 (MQ=255)                                     tGGTAGGTCGTCACGGTAACAAGGGTGTAGTTTCTAAGATCACCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTc                                                                                                    <  2:521990/143‑1 (MQ=255)                                     tGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTc                                                                                                    <  2:376258/143‑1 (MQ=255)                                              gtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCAc                                                                                           <  1:686040/143‑1 (MQ=255)                                               tcACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCAcc                                                                                          >  1:729803/1‑143 (MQ=255)                                                cACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCt                                                                                         <  1:203539/143‑1 (MQ=255)                                                      aaCAAGGGTGTAATTTCTAAGATCAACCCGATCGACGCTCTGCATTACCCTGCAAACGGCACGCCGGTAGACATGGTATTGGGCCAGCTGGGCGTCCCGT‑AAATAAGAGACAGCGGTCAGAACCTCGAAACTCACCTGGGTAt                                                                                   >  1:499568/1‑143 (MQ=255)                                                         aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                                >  1:202677/1‑143 (MQ=255)                                                         aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGc                                                                                >  1:634126/1‑143 (MQ=255)                                                         aaGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCGTACGATGAAAACGGTACGCCGGTAGACAGCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAGACCCGCCCGGGGGTGGc                                                                                >  1:520268/1‑143 (MQ=255)                                                               gtAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCGCGTATGA‑ACATCGGTCAGATCCACGAAACCCACCGGGGAATGgcggcggc                                                                          >  2:160481/1‑141 (MQ=255)                                                               atAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGaa                                                                          <  1:569033/142‑1 (MQ=255)                                                                           ctCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGa                                                              <  1:383451/142‑1 (MQ=255)                                                                              aaCCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACAGCGGTCAGATCCCCGAAACCCCCCTGGGTATGGCTGCGAAAGGTATCGGCGACaa                                                           <  1:770397/143‑1 (MQ=255)                                                                                   gATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCGGCTCAAGGTATCGGCGACAAGATCa                                                      >  2:677562/1‑143 (MQ=255)                                                                                           atatGCCTTCACATGAAAACGGTACGCCGCTAGACATCGTACTGAACCAGCCTGGCGTACCGTCTAGAATGA‑ACACCCGACAGTGCCTCGAAACCCACCTGGGTATGGCTGCGTAAGGTATCGGCGACAAGATCAACGCCATg                                              <  1:114384/143‑1 (MQ=255)                                                                                             atGCCTTACGATAAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCt                                            >  1:42717/1‑143 (MQ=255)                                                                                                 cTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGaaa                                        <  1:729476/143‑1 (MQ=255)                                                                                                     cGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAAcagc                                    <  1:25815/143‑1 (MQ=255)                                                                                                        tGAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACagcagc                                 >  1:337811/1‑143 (MQ=255)                                                                                                         gAAAACGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATGAACGCCATGCTGAAACAgcagca                                >  1:538517/1‑143 (MQ=255)                                                                                                          aaaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaa                               <  2:337811/143‑1 (MQ=255)                                                                                                           aaaCGGTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaag                              <  1:214108/143‑1 (MQ=255)                                                                                                             aCGGTACGCCGGGAGACATCGCACGGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGCACCCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCaagaa                            <  1:301262/143‑1 (MQ=255)                                                                                                                gTACGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCAGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTc                         <  1:416623/143‑1 (MQ=255)                                                                                                                 tACGCCGGTAGAAATAGTACTGAACCAGCCGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcg                        <  2:729803/143‑1 (MQ=255)                                                                                                                   cGCCGGTAGACATCGTCCTGACCCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgcg                      <  1:463588/143‑1 (MQ=255)                                                                                                                   cGCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTcgcg                      >  1:34174/1‑143 (MQ=255)                                                                                                                    gCCGGTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTAGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGa                     >  1:632837/1‑143 (MQ=255)                                                                                                                       ggTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAAc                  >  1:252983/1‑143 (MQ=255)                                                                                                                        gTAGACATCGTACTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACACACCTGGGTATGGCTGCGAAAGGTATCGGCGCCAAGCTCCACGCCCTGCTGAAACAGCAGCAAGAAGTCGCGAAAcc                 >  1:159242/1‑142 (MQ=255)                                                                                                                            aCATCGTACTGAACCAGCTGGGCGTCCCGTCTCGTATGA‑TCATCGGTCAGCTCCTCGAACCCCACCCGGGTATGGCTGCGACAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTgcgc             <  2:702277/143‑1 (MQ=255)                                                                                                                                    cTGAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCAt     >  1:662241/1‑143 (MQ=255)                                                                                                                                     tGACCCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc    <  1:564170/143‑1 (MQ=255)                                                                                                                                     tGAACCAGCTGGGCGTACCGTCTCGTATTT‑TCATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATc    >  2:487427/1‑143 (MQ=255)                                                                                                                                      gAACCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCGCGAAACCCACCTGGGTATGGCGGCGAAAGGTATCGGCGACAAGAACAAAGCCATGCTGAACCAGCAGCAAGAAGGCGCGAAACTGCGCGAAATCAccc   >  2:473456/1‑143 (MQ=255)                                                                                                                                       aaCCAGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGCATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCCAAACTGCGCGACATCATACa  >  2:798160/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                            TGGCGGTTAAACGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGTCTCGTATGA‑ACATCGGTCAGATCCTCGAAACCCACCTGGGTATGGCTGCGAAAGGTATCGGCGACAAGATCAACGCCATGCTGAAACAGCAGCAAGAAGTCGCGAAACTGCGCGAATTCATCCA  >  NC_000913/4184405‑4184681

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