Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,754,231	G→T	R127S (CGC→AGC)	selA ←	L‑seryl‑tRNA(Sec) selenium transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,754,231	0	G	T	100.0%	282.1 / NA	73	R127S (CGC→AGC)	selA	L‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (45/28);  total (45/28)

CATTGTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754150‑3754316                                                                                   |                                                                                      cATTGTTATTGATGATACACGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:941939/90‑1 (MQ=255) cATTGTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:1123475/90‑1 (MQ=255) cATTGTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:280772/90‑1 (MQ=255) cATTGTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:148061/90‑1 (MQ=255)     gTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCgtcgtc                                                                           >  1:1315543/1‑90 (MQ=255)      ttATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:142635/90‑1 (MQ=255)       tATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         <  2:248557/90‑1 (MQ=255)       tATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         >  1:1471203/1‑90 (MQ=255)       tATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCa                                                                         <  2:1315543/90‑1 (MQ=255)        aTTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAg                                                                        <  2:239303/90‑1 (MQ=255)          tGACGATACAGGCATCTTCCGCCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                >  1:1190096/1‑81 (MQ=255)          tGACGATACAGGCATCTTCCGCCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                <  2:1190096/81‑1 (MQ=255)          tGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  2:481499/1‑90 (MQ=255)          tGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  1:375233/1‑90 (MQ=255)             cGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCAt                                                                   >  1:656415/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 <  2:1909/90‑1 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:23172/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:265921/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:20801/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:1172423/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:382625/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:1070530/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:955131/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:1073419/1‑90 (MQ=255)               aTACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:748574/1‑90 (MQ=255)                   aGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             <  2:1021977/90‑1 (MQ=255)                       aTCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGt                                                         >  2:1430752/1‑90 (MQ=255)                         cTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  2:299715/1‑90 (MQ=255)                         cTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  2:593362/1‑90 (MQ=255)                         cTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  1:855813/1‑90 (MQ=255)                         cTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  2:94334/1‑90 (MQ=255)                          ttCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      >  1:585392/1‑90 (MQ=255)                          ttCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  2:60204/90‑1 (MQ=255)                          ttCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  2:78165/90‑1 (MQ=255)                          ttCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  1:740571/90‑1 (MQ=255)                              g‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             >  2:404673/1‑79 (MQ=255)                              g‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             <  1:404673/79‑1 (MQ=255)                                 ccccGTAATATGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:447349/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:784754/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:499968/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:520366/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:567571/90‑1 (MQ=255)                                   ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                                              >  2:935862/1‑90 (MQ=255)                                   ccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAAc                                              >  1:1159763/1‑90 (MQ=255)                                    cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACg                                             >  1:993403/1‑90 (MQ=255)                                      tAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                           <  2:599908/90‑1 (MQ=255)                                          aCGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                                       >  2:1137577/1‑90 (MQ=255)                                           cGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                      <  1:1020888/90‑1 (MQ=255)                                            ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCt                                     >  1:461656/1‑90 (MQ=255)                                            ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCt                                     >  2:326732/1‑90 (MQ=255)                                            ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCt                                     >  1:1238158/1‑90 (MQ=255)                                             gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  2:399765/1‑90 (MQ=255)                                             gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  1:540569/1‑90 (MQ=255)                                              cacaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                   >  1:1282329/1‑90 (MQ=255)                                              cacaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                   >  2:1109128/1‑90 (MQ=255)                                              cacaGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgc                                   <  2:1210561/90‑1 (MQ=255)                                                 agcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCa                                >  2:1255230/1‑90 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          >  1:1364691/1‑47 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:1364691/47‑1 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATa                                                                  >  1:1340515/1‑55 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATa                                                                  <  2:1340515/55‑1 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               >  1:495927/1‑90 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               >  2:1266699/1‑90 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               >  1:1585572/1‑90 (MQ=255)                                                     gcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACgg                            >  2:1584745/1‑90 (MQ=255)                                                       tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                           >  1:1093017/1‑73 (MQ=255)                                                       tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                           <  2:1093017/73‑1 (MQ=255)                                                       tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                          >  2:1529723/1‑90 (MQ=255)                                                       tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGCCGTCCAGATCATCCTCGGGGGTCACTGGCGAACGCATAGCCTGCGCAACCGCt                          >  1:215137/1‑90 (MQ=255)                                                        gcgcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCtt                         <  1:1255230/90‑1 (MQ=255)                                                                       cgcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCg          <  2:111946/90‑1 (MQ=255)                                                                             gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTa    >  1:992193/1‑90 (MQ=255)                                                                               ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  2:885205/1‑90 (MQ=255)                                                                                   |                                                                                      CATTGTTATTGACGATACAGGCATCTTCCG‑CCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754150‑3754316

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