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%	197.5 / NA	50	R127S (CGC→AGC)	selA	L‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/28);  total (22/28)

GCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754149‑3754316                                                                                   |                                                                                      gcATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcg                                                                                <  2:192089/90‑1 (MQ=255)  cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:98642/90‑1 (MQ=255)  cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:757590/90‑1 (MQ=255)  cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:446399/90‑1 (MQ=255)  cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:739857/90‑1 (MQ=255)  cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:741981/90‑1 (MQ=255)       ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:139218/90‑1 (MQ=255)       ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:576040/90‑1 (MQ=255)       ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:7280/90‑1 (MQ=255)                aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  2:10688/1‑90 (MQ=255)                aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:860898/1‑90 (MQ=255)                aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:675666/1‑90 (MQ=255)                   cAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCg                                                              >  2:255973/1‑90 (MQ=255)                    aGGCATCTACCGCCCCCGTAATACGGCACAGCCGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             <  1:956124/90‑1 (MQ=255)                     ggCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAg                                                            >  1:38560/1‑90 (MQ=255)                          cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  1:733629/1‑66 (MQ=255)                          cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:733629/66‑1 (MQ=255)                          cTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCa                                                       >  2:449370/1‑90 (MQ=255)                               gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  >  1:611503/1‑90 (MQ=255)                               gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  >  1:794528/1‑90 (MQ=255)                               gCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGc                                                  >  2:397650/1‑90 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:761779/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:881642/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:675666/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:860898/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:441971/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:255973/90‑1 (MQ=255)                                 ccccGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  2:780527/90‑1 (MQ=255)                                    cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACg                                             >  2:836117/1‑90 (MQ=255)                                     gTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGc                                            >  2:931050/1‑90 (MQ=255)                                          aCGGCACAGAAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGc                                       >  2:340381/1‑90 (MQ=255)                                           cGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                      <  2:338658/90‑1 (MQ=255)                                           cGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                      <  1:62067/90‑1 (MQ=255)                                           cGGCACAGAAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGcc                                      <  1:340381/90‑1 (MQ=255)                                            ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  1:720808/1‑48 (MQ=255)                                            ggCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:720808/48‑1 (MQ=255)                                             gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:762088/78‑1 (MQ=255)                                             gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                >  2:762088/1‑78 (MQ=255)                                                cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             <  2:108870/62‑1 (MQ=255)                                                cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGa                                                             >  1:108870/1‑62 (MQ=255)                                                cagcagCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTgcgc                                 >  1:598803/1‑90 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  2:354413/90‑1 (MQ=255)                                                  gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCaa                               <  2:611503/90‑1 (MQ=255)                                                       tGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCt                          >  2:273580/1‑90 (MQ=255)                                                                  ccGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGc               >  2:65065/1‑90 (MQ=255)                                                                    gATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             <  1:163341/90‑1 (MQ=255)                                                                        gcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGc         >  2:664307/1‑90 (MQ=255)                                                                             gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             <  1:966942/81‑1 (MQ=255)                                                                             gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCtt             >  2:966942/1‑81 (MQ=255)                                                                               ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTaaa  >  2:690205/1‑90 (MQ=255)                                                                                   |                                                                                      GCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754149‑3754316

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