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_CP0092733,002,9030TC100.0% 87.2 / NA 27S261P (TCA→CCA) ygeYYgeY family selenium metabolism‑linked hydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/14);  total (13/14)

TTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAG  >  NZ_CP009273/3002819‑3002984
                                                                                    |                                                                                 
tccgtCGGCAAAGGCACCCTCCCCGTTTCTGAAATCTTGTTCACATCCCAAAGCCGTTGCGCTGTAGCAGACAGCTTCGCAGTTCCAATc                                                                              <  2:293006/86‑1 (MQ=255)
    tCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGAcc                                                                          >  1:485617/1‑90 (MQ=255)
       gCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACcgtc                                                                       <  1:50660/90‑1 (MQ=255)
          aaGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    <  1:425538/90‑1 (MQ=255)
                    cACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTgggg                                                          >  1:1437/1‑90 (MQ=255)
                     aCCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGc                                                         <  1:399108/90‑1 (MQ=255)
                                     tcttcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcg                                         >  1:307184/1‑90 (MQ=255)
                                      cttcACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGgcgc                                        >  1:205784/1‑90 (MQ=255)
                                         cacaTCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCACAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCt                                     <  2:441606/90‑1 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                       >  1:184873/1‑68 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGa                                                       <  2:184873/68‑1 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAgg                                           >  2:344713/1‑80 (MQ=255)
                                             tCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAgg                                           <  1:344713/80‑1 (MQ=255)
                                              ccccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACg                                <  2:205784/90‑1 (MQ=255)
                                                ccAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGaa                              >  2:224295/1‑90 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    >  2:208228/1‑45 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCgtcgtc                                                                    <  1:208228/45‑1 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                       >  1:395948/1‑90 (MQ=255)
                                                       gTTGCGCTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCcgcg                       >  1:9879/1‑90 (MQ=255)
                                                            gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                  <  1:135974/90‑1 (MQ=255)
                                                            gcTGTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctg                  <  1:244538/90‑1 (MQ=255)
                                                               gTAGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCcctgcct               >  1:430921/1‑90 (MQ=255)
                                                                 aGCAGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             <  1:161827/90‑1 (MQ=255)
                                                                    aGACAGCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGt          <  1:224295/90‑1 (MQ=255)
                                                                         gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             <  2:225973/82‑1 (MQ=255)
                                                                         gCTGCGCAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCctgcctgc             >  1:225973/1‑82 (MQ=255)
                                                                            gcgcAGTTCCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAg  >  2:268852/1‑90 (MQ=255)
                                                                                    |                                                                                 
TTCCTCGGCAAAGGCACCCTCACCGTTTCTGAAATCTTCTTCACATCCCCAAGCCGTTGCGCTGTAGCAGACAGCTGCGCAGTTTCAATCGACCGTCGTCTGACCTGGGGCGAAACCTGGGAAGGCGCGCTGGACGAAATCCGCGCCCTGCCTGCAGTACAGAAAG  >  NZ_CP009273/3002819‑3002984

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