Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,477,150	C→T	K90K (AAG→AAA)	fkpA ←	peptidyl‑prolyl cis‑trans isomerase FkpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,477,150	0	C	T	100.0%	108.7 / NA	31	K90K (AAG→AAA)	fkpA	peptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/17);  total (14/17)

CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196                                               |                                               cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                                              <  1:835769/49‑1 (MQ=255) cGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                                             >  1:1853987/1‑50 (MQ=255)  gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            >  1:249853/1‑50 (MQ=255)  gAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                                            <  1:2102448/50‑1 (MQ=255)    aTGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCa                                          <  1:855755/50‑1 (MQ=255)       cTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGc                                        <  1:2148656/49‑1 (MQ=255)        ttGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                                       >  1:919136/1‑49 (MQ=255)         tGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                                     >  1:1252544/1‑50 (MQ=255)           tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                                    <  1:488705/49‑1 (MQ=255)           tAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                                   >  1:1244507/1‑50 (MQ=255)             gagTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                                 >  1:1951493/1‑50 (MQ=255)               gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               <  1:2231672/50‑1 (MQ=255)               gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               <  1:1654111/50‑1 (MQ=255)               gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                               <  1:31216/50‑1 (MQ=255)                    ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg                            <  1:2487385/48‑1 (MQ=255)                      cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAac                         <  1:171876/49‑1 (MQ=255)                       gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:130306/1‑50 (MQ=255)                        aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAacac                       >  1:2478245/1‑49 (MQ=255)                           tctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   >  1:1374080/1‑50 (MQ=255)                             ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:1941167/47‑1 (MQ=255)                             ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAg                    <  1:1543804/47‑1 (MQ=255)                             ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGc                   <  1:2562146/48‑1 (MQ=255)                               ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCg                  <  1:253302/47‑1 (MQ=255)                               ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  1:1789763/1‑49 (MQ=255)                               ggTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGAt                >  1:1124940/1‑49 (MQ=255)                                gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCa              <  1:893184/50‑1 (MQ=255)                                  cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGc            >  1:241342/1‑50 (MQ=255)                                      gagTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGAt        <  1:2547680/50‑1 (MQ=255)                                        gTTTGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCt      <  1:2572534/50‑1 (MQ=255)                                         tttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt     <  1:493856/50‑1 (MQ=255)                                          ttGCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCtt     >  1:21440/1‑49 (MQ=255)                                            gCTTTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAt  >  1:1827062/1‑50 (MQ=255)                                               |                                               CGAATGCTTGTAGAGTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGATCAGCTGATCTTTAT  >  NC_000913/3477104‑3477196

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑