Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,002,903	T→C	S261P (TCA→CCA)	ygeY →	YgeY family selenium metabolism‑linked hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,002,903	0	T	C	100.0%	212.9 / NA	61	S261P (TCA→CCA)	ygeY	YgeY family selenium metabolism‑linked hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/30);  total (31/30)

GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCT  >  NZ_CP009273/3002816‑3002986                                                                                        |                                                                                    gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa                                                                                   >  1:236872/1‑90 (MQ=255) gAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCa                                                                                   >  2:114525/1‑90 (MQ=255)  aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa                                                                                  >  1:11038/1‑90 (MQ=255)  aaTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCaa                                                                                  >  2:264687/1‑90 (MQ=255)   aTTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAAt                                                                                 >  1:322254/1‑90 (MQ=255)        tCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGAcc                                                                            >  2:1233459/1‑90 (MQ=255)          ggCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgt                                                                          <  2:627716/90‑1 (MQ=255)          ggCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgt                                                                          <  1:1233459/90‑1 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      >  2:1075146/1‑90 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      <  2:977381/90‑1 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      >  2:546641/1‑90 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      <  2:28491/90‑1 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      <  1:608965/90‑1 (MQ=255)              aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                      <  1:268116/90‑1 (MQ=255)                   aCCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGAcc                                                                 >  2:1088525/1‑90 (MQ=255)                     ccTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg                                                               >  1:81018/1‑90 (MQ=255)                     ccTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTg                                                               >  2:570613/1‑90 (MQ=255)                        cACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTgggg                                                            >  1:138695/1‑90 (MQ=255)                          ccGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCg                                                          <  2:566762/90‑1 (MQ=255)                          ccGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCg                                                          <  1:725852/90‑1 (MQ=255)                           cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                         <  2:847802/90‑1 (MQ=255)                           cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                         <  2:216243/90‑1 (MQ=255)                           cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                         <  2:1007833/90‑1 (MQ=255)                           cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                         <  2:82121/90‑1 (MQ=255)                           cGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                         <  2:127957/90‑1 (MQ=255)                            gTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaa                                                        <  2:337839/90‑1 (MQ=255)                             tttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGaaa                                                       <  2:160043/90‑1 (MQ=255)                              ttCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAAc                                                      >  1:144959/1‑90 (MQ=255)                                       cttcttCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAgg                                             >  2:124701/1‑90 (MQ=255)                                          cttcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgc                                          >  1:618985/1‑90 (MQ=255)                                             cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt                                       <  2:40514/90‑1 (MQ=255)                                                           gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcg                                                                           <  2:1227647/40‑1 (MQ=39)                                                           gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcg                                                                           >  1:1227647/1‑40 (MQ=39)                                                           gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                         >  2:110066/1‑90 (MQ=255)                                                           gTTGCGCTGTAGCAGACAGCTGCACAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                         >  1:470701/1‑90 (MQ=255)                                                             tGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGcc                       <  1:680939/90‑1 (MQ=255)                                                                gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                    <  1:360242/90‑1 (MQ=255)                                                                gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                    >  1:1095265/1‑90 (MQ=255)                                                                  tGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcc                  <  2:363916/90‑1 (MQ=255)                                                                   gTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcct                 >  2:1116857/1‑90 (MQ=255)                                                                     aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  1:1088525/90‑1 (MQ=255)                                                                     aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  1:945330/90‑1 (MQ=255)                                                                     aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  2:144959/90‑1 (MQ=255)                                                                     aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  1:48895/90‑1 (MQ=255)                                                                     aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  1:483580/90‑1 (MQ=255)                                                                      gCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCa              <  1:197562/90‑1 (MQ=255)                                                                       cAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAg             <  1:1170523/90‑1 (MQ=255)                                                                           cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa         >  2:1143082/1‑90 (MQ=255)                                                                           cAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACa         >  1:392691/1‑90 (MQ=255)                                                                             gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa       >  2:1045589/1‑90 (MQ=255)                                                                             gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa       >  1:877659/1‑90 (MQ=255)                                                                             gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa       >  2:643659/1‑90 (MQ=255)                                                                             gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGa       >  2:738923/1‑90 (MQ=255)                                                                               tGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  2:638748/80‑1 (MQ=255)                                                                               tGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               >  1:638748/1‑80 (MQ=255)                                                                               tGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGaaa     >  1:435553/1‑90 (MQ=255)                                                                                gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               >  1:831153/1‑79 (MQ=255)                                                                                gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc               <  2:831153/79‑1 (MQ=255)                                                                                gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg    >  1:370788/1‑90 (MQ=255)                                                                                gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg    >  2:571306/1‑90 (MQ=255)                                                                                  gcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCt  <  2:522089/90‑1 (MQ=255)                                                                                        |                                                                                    GAATTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAGCT  >  NZ_CP009273/3002816‑3002986

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