Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 601,986 A→G F285S (TTC→TCC)  ybdK ← YbdK family carboxylate‑amine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273601,9860AG87.5% 54.5 / 2.3 24F285S (TTC→TCC) ybdKYbdK family carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (13/8);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTG  >  NZ_CP009273/601899‑602068
                                                                                       |                                                                                  
aaTCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGgaa                                                                                  <  2:831075/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                             <  2:333181/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                             <  1:273775/90‑1 (MQ=255)
           aTCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAAttt                                                                       >  1:217606/1‑90 (MQ=255)
                  gTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATAcagc                                                                >  2:659093/1‑90 (MQ=255)
                              cgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCtttt                                                    <  2:712291/90‑1 (MQ=255)
                                    ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGa                                              <  2:438979/90‑1 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                         >  2:791666/1‑90 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                         >  1:511052/1‑90 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                         >  2:343847/1‑90 (MQ=255)
                                           gCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTaa                                       >  1:729384/1‑90 (MQ=255)
                                                 cGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGc                                 >  2:681561/1‑90 (MQ=255)
                                                     gatgaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTc                             >  1:506569/1‑90 (MQ=255)
                                                        gaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGt                          <  2:297884/90‑1 (MQ=255)
                                                                cAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAAcc                  <  1:266494/90‑1 (MQ=255)
                                                                       aTAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGc           >  1:734306/1‑90 (MQ=255)
                                                                         aGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg         <  1:214319/90‑1 (MQ=255)
                                                                         aGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg         <  2:606051/90‑1 (MQ=255)
                                                                          gCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGt        <  1:623627/90‑1 (MQ=255)
                                                                            gACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAg      >  1:86966/1‑90 (MQ=255)
                                                                            gACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAg      >  1:277878/1‑90 (MQ=255)
                                                                              cAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGcc    <  2:367643/90‑1 (MQ=255)
                                                                              cAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGcc    >  2:826053/1‑90 (MQ=255)
                                                                                ggCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTg  >  2:426051/1‑90 (MQ=255)
                                                                                ggCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTg  >  2:28322/1‑90 (MQ=255)
                                                                                       |                                                                                  
AATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTG  >  NZ_CP009273/601899‑602068

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