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_CP0092733,754,2310GT100.0% 76.2 / NA 21R127S (CGC→AGC) selAL‑seryl‑tRNA(Sec) selenium transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (11/10);  total (11/10)

CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754150‑3754316
                                                                                 |                                                                                     
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:468560/90‑1 (MQ=255)
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  1:335697/90‑1 (MQ=255)
cATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:211451/90‑1 (MQ=255)
     ttATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTcc                                                                          <  2:364851/90‑1 (MQ=255)
         tGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGAt                                                                      >  2:190836/1‑90 (MQ=255)
            cGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCAt                                                                   >  2:38501/1‑90 (MQ=255)
              aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGCTGTCCGCTTCCGGCGTCGTCCAGATCAtcc                                                                 >  2:407546/1‑88 (MQ=255)
              aTACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATAc                                                                 >  1:581363/1‑90 (MQ=255)
                         ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      >  2:538457/1‑90 (MQ=255)
                         ttCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCAc                                                      <  1:364958/90‑1 (MQ=255)
                               ccccGTAATACGGCACAGCAGCTGCGCCCGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGa                                                <  1:157608/90‑1 (MQ=255)
                                  cGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACg                                             >  1:213221/1‑90 (MQ=255)
                                           gCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTg                                    >  1:523238/1‑90 (MQ=255)
                                                gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               >  1:408174/1‑42 (MQ=255)
                                                gcagcTGCGCCAGCGCCCGATCGCGATGTCCGCTTCCGGcgt                                                                               <  2:408174/42‑1 (MQ=255)
                                                        gcCAGCGCCCGATCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTcc                       >  2:154913/1‑90 (MQ=255)
                                                                    tCGCGATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTcc           <  2:43195/90‑1 (MQ=255)
                                                                     cgcgATGTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCg          <  2:523238/90‑1 (MQ=255)
                                                                           gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                           >  1:513953/1‑51 (MQ=255)
                                                                           gTCCGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCa                                           <  2:513953/51‑1 (MQ=255)
                                                                             ccGCTTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCGCCGCGGCTTCCGCCTGTaaa  >  2:440451/1‑90 (MQ=255)
                                                                                 |                                                                                     
CATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATGTCCGCGTCCGGCGTCGTCCAGATCATACTCGAGGGTCACTGGCGAACGCATAGCCTGCGCAACGGCTTCCACCGCGGCTTCCGCCTGTAAA  >  NZ_CP009273/3754150‑3754316

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