Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,220,044 A→G S236P (TCG→CCG)  yohF ← SDR family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,220,0440AG94.3% 150.6 / 1.6 53S236P (TCG→CCG) yohFSDR family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (3/0);  new base G (28/22);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

ACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219957‑2220129
                                                                                       |                                                                                     
aCCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCt                                                                                     >  1:245712/1‑90 (MQ=255)
   cAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGcc                                                                                  >  2:178047/1‑90 (MQ=255)
   cAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGcc                                                                                  >  1:1555113/1‑90 (MQ=255)
    aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:1011033/90‑1 (MQ=255)
    aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:874616/90‑1 (MQ=255)
    aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:594766/90‑1 (MQ=255)
    aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  2:582279/90‑1 (MQ=255)
    aaTGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGccc                                                                                 <  1:1428269/90‑1 (MQ=255)
      tGAAGTATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCTATCCACTATCAACGGCTGCCCgg                                                                               <  1:1000453/90‑1 (MQ=255)
      tGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCgg                                                                               <  2:408410/90‑1 (MQ=255)
                      cTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCNCTGTGGTGTAATTTGCGcc                                                               >  1:107152/1‑90 (MQ=255)
                       tCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGccc                                                              >  1:197025/1‑90 (MQ=255)
                         aCTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTc                                                            <  2:751830/90‑1 (MQ=255)
                              gcTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAAc                                                       <  1:47087/90‑1 (MQ=255)
                               cTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACa                                                      >  2:1491931/1‑90 (MQ=255)
                               cTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACa                                                      >  2:1169205/1‑90 (MQ=255)
                                 aTTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACaaa                                                    <  1:990038/90‑1 (MQ=255)
                                      gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  2:985653/1‑90 (MQ=255)
                                      gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               <  1:772072/90‑1 (MQ=255)
                                      gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               <  2:1506772/90‑1 (MQ=255)
                                      gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  2:171904/1‑90 (MQ=255)
                                      gggTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               <  2:245712/90‑1 (MQ=255)
                                        gTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCacac                                             <  2:747424/90‑1 (MQ=255)
                                           gAACTGTGGATTCGCCAACGTAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:1082876/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcgc                                                                <  2:478362/68‑1 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTgcgc                                                                >  1:478362/1‑68 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:1487354/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:748617/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:707899/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:118832/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:272920/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:1080044/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:462877/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  2:899457/1‑90 (MQ=255)
                                           gAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCAcaccac                                          >  1:1024202/1‑90 (MQ=255)
                                           gAACAGTGGATTCGCCAACATAAAGCCGGCATCGACTATCATCGGCTGCCCGGTTTTGGAATTAGCGTCCTCCGAAGAAATCCAcatcac                                          >  1:38014/1‑90 (MQ=255)
                                                gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                     >  2:603511/1‑90 (MQ=255)
                                                gtgGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGc                                     >  2:1302617/1‑90 (MQ=255)
                                                     ttCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCa                                >  1:247593/1‑90 (MQ=255)
                                                      tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  2:225784/90‑1 (MQ=255)
                                                      tCGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCaa                               <  1:9515/90‑1 (MQ=255)
                                                       cGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               >  2:573204/1‑73 (MQ=255)
                                                       cGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCcac                                               <  1:573204/73‑1 (MQ=255)
                                                       cGCCAACATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAAt                              >  1:271371/1‑90 (MQ=255)
                                                            aCATAAAGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCAt                         >  1:1361559/1‑90 (MQ=255)
                                                                 aaGCCGCCATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGtt                    <  2:153345/90‑1 (MQ=255)
                                                                   gccgccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTgc                  <  2:748617/90‑1 (MQ=255)
                                                                      gccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGcc               >  1:209739/1‑90 (MQ=255)
                                                                       ccATCCACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCa              >  1:23851/1‑90 (MQ=255)
                                                                           ccACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAAc          <  2:37025/90‑1 (MQ=255)
                                                                           ccACTATCAACGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAAc          <  2:124216/90‑1 (MQ=255)
                                                                                   aaCGGCTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  1:219574/1‑90 (MQ=255)
                                                                                   aaCGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAg  >  1:139985/1‑90 (MQ=255)
                                                                                       |                                                                                     
ACCCAATGAAATATCCACGCACCTCACTGCGCTATTCTGGGTTGAACTGTGGATTCGCCAACATAAAGCCGCCATCCACTATCAACGACTGCCCGGTGGTGTAATTTGCGCCCTCCGAACAAAGCCACACCACCAGGCTGGCAATCTCATGCGTTGCGCCAAAACGCCGCAAG  >  NZ_CP009273/2219957‑2220129

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