Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 259,276 A→G G35G (GGT→GGC ykfI ← type IV toxin‑antitoxin system toxin YkfI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273259,2760AG94.7% 115.6 / ‑2.4 38G35G (GGT→GGCykfItype IV toxin‑antitoxin system toxin YkfI
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (16/20);  total (17/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TCACGGCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGACCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCACCGCCCGCTGAGT  >  NZ_CP009273/259191‑259359
                                                                                     |                                                                                   
tCACGGCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCAt                                                                                 <  2:365657/90‑1 (MQ=255)
    ggCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTg                                                                             <  2:907304/90‑1 (MQ=255)
    ggCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTg                                                                             <  1:792679/90‑1 (MQ=255)
     gCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGc                                                                            <  1:814051/90‑1 (MQ=255)
                ggTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCcagcag                                                                 <  2:135557/90‑1 (MQ=255)
                  tGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGAc                                                               <  2:319089/90‑1 (MQ=255)
                   gATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACg                                                              >  2:274666/1‑90 (MQ=255)
                   gATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACg                                                              <  1:4783/90‑1 (MQ=255)
                     tACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACgtg                                                            <  1:328998/90‑1 (MQ=255)
                         ggCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATTGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAg                                                        >  1:916172/1‑90 (MQ=255)
                              cGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACa                                                   >  1:20196/1‑90 (MQ=255)
                                   tgtTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGc                                              >  2:556499/1‑90 (MQ=255)
                                   tgtTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGc                                              >  2:540001/1‑90 (MQ=255)
                                        cTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGAc                                         >  1:911015/1‑90 (MQ=255)
                                              cacaGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCAc                                   <  2:972442/90‑1 (MQ=255)
                                              cacaGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCAc                                   <  2:453298/90‑1 (MQ=255)
                                                caGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACgg                                 >  1:25815/1‑90 (MQ=255)
                                                caGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAATGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCCCgg                                 >  1:499500/1‑90 (MQ=21)
                                                caGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGACCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACgg                                 >  1:785692/1‑90 (MQ=255)
                                                 aGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGAc                                                               >  1:869702/1‑59 (MQ=255)
                                                 aGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGAc                                                               <  2:869702/59‑1 (MQ=255)
                                                         tGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGAcaggc                        <  1:180572/90‑1 (MQ=255)
                                                         tGCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGAcaggc                        <  1:84759/90‑1 (MQ=255)
                                                          gCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACg                                  >  1:898331/1‑79 (MQ=255)
                                                          gCAGAATGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACg                                  <  2:898331/79‑1 (MQ=255)
                                                            aGAATGGCGTGTCGTTTATTGTCAGACCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGAcaggcagg                     <  2:785692/90‑1 (MQ=255)
                                                               aTGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCt                  >  1:291105/1‑90 (MQ=255)
                                                                tGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACAtt                                                 >  2:656154/1‑58 (MQ=255)
                                                                tGGCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACAtt                                                 <  1:656154/58‑1 (MQ=255)
                                                                  gCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCa               >  2:128085/1‑90 (MQ=255)
                                                                  gCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCa               >  1:908165/1‑90 (MQ=255)
                                                                  gCGTGTCGTTTATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCa               >  1:744780/1‑90 (MQ=255)
                                                                       tCGTTTATTGTCAGGCCATAGAGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCAccgcc          >  2:541335/1‑90 (MQ=255)
                                                                          tttATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCAccgcccgc       <  1:359210/90‑1 (MQ=255)
                                                                          tttATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCAccgcccgc       <  2:17259/90‑1 (MQ=255)
                                                                           ttATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCACCGCCCGCt      <  2:932925/90‑1 (MQ=255)
                                                                           ttATTGTCAGGCCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCACCGCCCGCt      <  2:952596/90‑1 (MQ=255)
                                                                               tGTCAGGCCATAGAGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCACCGCCCGCTGAGt  <  1:541335/90‑1 (MQ=255)
                                                                                     |                                                                                   
TCACGGCATCGGCTAGGGTGATACCGGCATCGATGTGTTCCTTAATCACAGCCTCATTGCAGAATGGCGTGTCGTTTATTGTCAGACCATAGTGCTGTTCCAGCAGACGTGTCAGTAACATTTGCCAGACAGCCACGGGTGACAGGCAGGGCTTCACCGCCCGCTGAGT  >  NZ_CP009273/259191‑259359

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