Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,011,010 A→G E531G (GAA→GGA)  ygfK → putative selenate reductase subunit YgfK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,011,0100AG96.7% 187.1 / ‑1.0 61E531G (GAA→GGA) ygfKputative selenate reductase subunit YgfK
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (35/24);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGC  >  NZ_CP009273/3010923‑3011091
                                                                                       |                                                                                 
tttGCGACCCCCAGTGCAAATAAAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGgaa                                                                                 <  2:446829/90‑1 (MQ=255)
        cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                         >  2:596417/1‑90 (MQ=255)
        cACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTgg                                                                         >  1:1097561/1‑90 (MQ=255)
          ccAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                       <  2:373366/90‑1 (MQ=255)
          ccAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                       >  1:593896/1‑90 (MQ=255)
          ccAGTGCAAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGa                                                                       <  2:223173/90‑1 (MQ=255)
           cAGTGCAAATACACCTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGAt                                                                      <  2:122779/90‑1 (MQ=255)
            agccccAATACAACTGTCCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATg                                                                     <  1:202872/86‑1 (MQ=255)
                  aaTACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAAtata                                                               >  2:1027052/1‑90 (MQ=255)
                    tACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAg                                                             >  2:834416/1‑90 (MQ=255)
                       aaCTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCaa                                                          >  1:366558/1‑90 (MQ=255)
                          tGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc                                                       <  2:366558/90‑1 (MQ=255)
                           gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc                                                       <  2:282819/89‑1 (MQ=255)
                           gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGc                                                       >  1:282819/1‑89 (MQ=255)
                           gTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCt                                                      >  2:769650/1‑90 (MQ=255)
                             aCCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTgg                                                    >  1:1355953/1‑90 (MQ=255)
                                cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 <  2:765179/90‑1 (MQ=255)
                                cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 >  2:895436/1‑90 (MQ=255)
                                cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 >  2:1354056/1‑90 (MQ=255)
                                cGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAATCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 >  1:377987/1‑90 (MQ=255)
                                 gCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcaca                                                <  2:568329/90‑1 (MQ=255)
                                     ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGc                                                            <  1:310073/74‑1 (MQ=255)
                                     ggATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGc                                                            >  2:310073/1‑74 (MQ=255)
                                         tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                        <  2:377574/90‑1 (MQ=255)
                                         tACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                        <  1:895436/90‑1 (MQ=255)
                                          aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                        <  1:843610/89‑1 (MQ=255)
                                          aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGcc                                        >  2:843610/1‑89 (MQ=255)
                                          aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                       >  1:1420228/1‑90 (MQ=255)
                                          aCGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCg                                       >  2:1275111/1‑90 (MQ=255)
                                              cAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTc                                   <  2:1210600/90‑1 (MQ=255)
                                               aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 >  1:1440294/1‑75 (MQ=255)
                                               aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGcac                                                 <  2:1440294/75‑1 (MQ=255)
                                               aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                  >  1:405660/1‑90 (MQ=255)
                                               aGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCt                                  >  2:1088307/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                      <  2:348983/53‑1 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                      <  1:129448/53‑1 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                      >  1:348983/1‑53 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGAt                                                                      >  2:129448/1‑53 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  2:320554/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  1:1186552/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  2:713560/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  1:268264/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  2:1027028/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTg                                 >  1:944630/1‑90 (MQ=255)
                                                gTGCGCTGAATATCCGCGAACTGAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTgg                                >  2:632071/1‑90 (MQ=255)
                                                  gcgcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTggg                               >  1:370873/1‑89 (MQ=255)
                                                  gcgcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGt                               <  1:989530/90‑1 (MQ=255)
                                                    gcTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTc                             <  2:1033409/90‑1 (MQ=255)
                                                     cTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                           <  2:235797/59‑1 (MQ=255)
                                                     cTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCa                                                           >  1:235797/1‑59 (MQ=255)
                                                         atatCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGcc                        >  2:314056/1‑90 (MQ=255)
                                                           atCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCAt                      <  1:769650/90‑1 (MQ=255)
                                                            tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                     >  2:278564/1‑90 (MQ=255)
                                                            tCCGCGAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATc                     >  2:1260958/1‑90 (MQ=255)
                                                               gcgAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCgg                  >  1:130506/1‑90 (MQ=255)
                                                                 gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt                >  1:1155891/1‑90 (MQ=255)
                                                                 gAACTGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGtt                >  1:1118305/1‑90 (MQ=255)
                                                                     tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg            <  2:293355/90‑1 (MQ=255)
                                                                     tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg            <  1:320554/90‑1 (MQ=255)
                                                                     tGAAAAAAGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCg            <  1:1185317/90‑1 (MQ=255)
                                                                          aaaGTCGCGCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGAtt       <  1:1275111/90‑1 (MQ=255)
                                                                               cgcgCTGGGAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGc  <  1:871305/90‑1 (MQ=255)
                                                                                       |                                                                                 
TTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATCCGGTTGCCGTGATTGGTGC  >  NZ_CP009273/3010923‑3011091

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: