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%	353.5 / NA	84	K90K (AAG→AAA)	fkpA	peptidyl‑prolyl cis‑trans isomerase FkpA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/55);  total (29/55)

GTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGA  >  NC_000913/3477118‑3477181                                 |                                gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                               >  1:995703/1‑35 (MQ=255) gTCTGTTCGATCTCTTGGTCGGAGAGTTTGCtttt                               >  1:6544119/1‑35 (MQ=255)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              >  1:1182816/1‑35 (MQ=38)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              >  1:7485780/1‑35 (MQ=255)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              <  1:4677562/35‑1 (MQ=255)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              <  1:1311925/35‑1 (MQ=255)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              <  1:7001535/35‑1 (MQ=255)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              >  1:132120/1‑35 (MQ=38)  tCTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTa                              <  1:1452527/35‑1 (MQ=38)   cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                             <  1:3382186/35‑1 (MQ=255)   cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                             <  1:3501863/35‑1 (MQ=255)   cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                             <  1:3764867/35‑1 (MQ=38)   cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                             <  1:6355850/35‑1 (MQ=255)   cTGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTAt                             <  1:715717/35‑1 (MQ=38)    tGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                            >  1:4827323/1‑35 (MQ=39)    tGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                            <  1:3956191/35‑1 (MQ=38)    tGTTCGATCTCTTGGTCGGAGAGTTTGCTTTTATc                            <  1:1209846/35‑1 (MQ=255)      ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                          <  1:4693476/35‑1 (MQ=255)      ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                          <  1:5643762/35‑1 (MQ=39)      ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                          >  1:40334/1‑35 (MQ=39)      ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                          <  1:7738451/35‑1 (MQ=38)      ttCGATCTCTTGGTCGGAGAGTTTGCTTTTATCAg                          >  1:7844114/1‑35 (MQ=39)        cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                        <  1:3488942/35‑1 (MQ=255)        cGATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCa                        <  1:6104987/35‑1 (MQ=255)         gATCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaa                       >  1:770584/1‑35 (MQ=39)          aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                      >  1:321569/1‑35 (MQ=255)          aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                      >  1:599823/1‑35 (MQ=39)          aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                      >  1:1088498/1‑35 (MQ=39)          aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                      <  1:4158264/35‑1 (MQ=255)          aTCTCTTGGTCGGAGAGTTTGCTTTTATCAGCaaa                      >  1:747002/1‑35 (MQ=39)           tctcTTGGTCGGAGAGTTTGCTTTTATCAGCAAAt                     >  1:3451069/1‑35 (MQ=39)            ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                    <  1:1236195/35‑1 (MQ=38)            ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                    <  1:1927423/35‑1 (MQ=255)            ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                    <  1:5157494/35‑1 (MQ=38)            ctctTGGTCGGAGAGTTTGCTTTTATCAGCAAATg                    <  1:1862549/35‑1 (MQ=38)            ctctTGGTCGAAGAGTTTGCTTTTATCAGCAAATg                    <  1:1862502/35‑1 (MQ=37)              ctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                  >  1:3665429/1‑35 (MQ=38)              ctTGGTCGGAGAGTTTGCTTTTATCAGCAAATGCa                  <  1:7634089/35‑1 (MQ=38)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 >  1:2545555/1‑35 (MQ=255)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:7262496/35‑1 (MQ=255)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:7023739/35‑1 (MQ=255)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:6378553/35‑1 (MQ=255)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:638771/35‑1 (MQ=38)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:3558110/35‑1 (MQ=38)               ttGGTCGGAGAGTTTGCTTTTATCAGCAAATGCAt                 <  1:360018/35‑1 (MQ=255)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:4364779/35‑1 (MQ=255)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:6448279/35‑1 (MQ=38)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:2207067/35‑1 (MQ=255)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:1639945/35‑1 (MQ=255)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:5387257/35‑1 (MQ=38)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:4965367/35‑1 (MQ=38)                tGGTCGGAGAGTTTGCTTTTATCAGCAAATGCATc                <  1:1640019/35‑1 (MQ=38)                  gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt              >  1:1298410/1‑35 (MQ=255)                  gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt              <  1:2051009/35‑1 (MQ=38)                  gTCGGAGAGTTTGCTTTTATCAGCAAATGCATCCt              >  1:6715632/1‑35 (MQ=255)                   tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg             <  1:5866036/35‑1 (MQ=38)                   tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg             >  1:4117558/1‑35 (MQ=255)                   tCGGAGAGTTTGCTTTTATCAGCAAATGCATCCTg             <  1:6592224/35‑1 (MQ=38)                    cGGAGAGTTTGCTTTTATCAGCAAATGCATCCTGa            <  1:3579281/35‑1 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:878372/35‑1 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:786731/35‑1 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           >  1:1314295/1‑35 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           >  1:6802972/1‑35 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:1028836/35‑1 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:1147912/35‑1 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:5376447/35‑1 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           >  1:3143612/1‑35 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           >  1:3673459/1‑35 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:4760608/35‑1 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           >  1:1382884/1‑35 (MQ=38)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:5376482/35‑1 (MQ=255)                     ggAGAGTTTGCTTTTATCAGCAAATGCATCCTGaa           <  1:1604256/35‑1 (MQ=255)                      gagaGTTTGCTTTTATCAGCAAATGCATCCTGAac          >  1:5820341/1‑35 (MQ=38)                       agagTTTGCTTTTATCAGCAAATGCATCCTGAaca         <  1:7124698/35‑1 (MQ=38)                       agagTTTGCTTTTATCAGCAAATGCATCCTGAaca         >  1:3849367/1‑35 (MQ=255)                        gagTTTGCTTTTATCAGCAAATGCATCCTGAacac        <  1:491001/35‑1 (MQ=38)                        gagTTTGCTTTTATCAGCAAATGCATCCTGAacac        <  1:427835/35‑1 (MQ=38)                         agTTTGCTTTTATCAGCAAATGCATCCTGAACAcc       <  1:5297527/35‑1 (MQ=38)                         agTTTGCTTTTATCAGCAAATGCATCCTGAACAcc       <  1:1803712/35‑1 (MQ=38)                          gTTTGCTTTTATCAGCAAATGCATCCTGAACACCa      >  1:5576842/1‑35 (MQ=38)                           tttGCTTTTATCAGCAAATGCATCCTGAACACCAg     <  1:1006406/35‑1 (MQ=255)                           tttGCTTTTATCAGCAAATGCATCCTGAACACCAg     >  1:6632872/1‑35 (MQ=38)                            ttGCTTTTATCAGCAAATGCATCCTGAACACCAGc    >  1:797317/1‑35 (MQ=38)                             tGCTTTTATCAGCAAATGCATCCTGAACACCAGCg   >  1:708700/1‑35 (MQ=38)                             tGCTTTTATCAGCAAATGCATCCTGAACACCAGCg   >  1:6228422/1‑35 (MQ=38)                              gCTTTTATCAGCAAATGCATCCTGAACACCAGCGa  <  1:6736038/35‑1 (MQ=255)                                 |                                GTCTGTTCGATCTCTTGGTCGGAGAGTTTGCTCTTATCAGCAAATGCATCCTGAACACCAGCGA  >  NC_000913/3477118‑3477181

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