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,9860AG97.3% 119.3 / ‑3.8 37F285S (TTC→TCC) ybdKYbdK family carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (13/23);  total (13/24)
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.81e-01

AATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCT  >  NZ_CP009273/601899‑602067
                                                                                       |                                                                                 
aaTCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGgaa                                                                                 <  1:754916/90‑1 (MQ=255)
 aTCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAAc                                                                                <  2:888379/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                            <  2:765711/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                            <  2:531019/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                            <  1:315141/90‑1 (MQ=255)
     cAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGtt                                                                            <  2:172106/90‑1 (MQ=255)
                 cGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATAcag                                                                <  1:147156/90‑1 (MQ=255)
                  gTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATAcagc                                                               <  1:539373/90‑1 (MQ=255)
                      gCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGt                                                           <  2:469571/90‑1 (MQ=255)
                            gacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCtt                                                     <  1:682384/90‑1 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                        >  2:120319/1‑90 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                        >  1:99909/1‑90 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                        >  1:58181/1‑90 (MQ=255)
                                         gtgCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGttt                                        >  1:2482/1‑90 (MQ=255)
                                           gCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTaa                                      >  1:424097/1‑90 (MQ=255)
                                           gCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTaa                                      >  1:414036/1‑90 (MQ=255)
                                            cGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTaaa                                     >  1:77383/1‑90 (MQ=255)
                                                   gTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGt                              >  2:278683/1‑90 (MQ=255)
                                                     gatgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTc                            <  2:296713/90‑1 (MQ=255)
                                                      atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTcc                           <  2:99909/90‑1 (MQ=255)
                                                      atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTcc                           <  2:2482/90‑1 (MQ=255)
                                                        gaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGt                         <  2:77383/90‑1 (MQ=255)
                                                        gaCGCCTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGt                         >  2:486918/1‑90 (MQ=255)
                                                             cTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTcc                           <  2:405691/83‑1 (MQ=255)
                                                             cTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTcc                           >  1:405691/1‑83 (MQ=255)
                                                             cTTCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTa                    >  1:152382/1‑90 (MQ=255)
                                                               tCAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAAc                  <  2:999424/90‑1 (MQ=255)
                                                                cAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAAcc                 <  1:120319/90‑1 (MQ=255)
                                                                cAAGCCCATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAAcc                 <  2:797564/90‑1 (MQ=255)
                                                                      cATAGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTggg           >  2:316110/1‑90 (MQ=255)
                                                                         aGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg        <  2:734647/90‑1 (MQ=255)
                                                                         aGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg        <  1:828469/90‑1 (MQ=255)
                                                                         aGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg        <  1:479136/90‑1 (MQ=255)
                                                                         aGCGACAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCgg        <  1:486918/90‑1 (MQ=255)
                                                                              cAGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGc  <  2:815734/90‑2 (MQ=255)
                                                                               aGGCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCt  <  1:157545/90‑1 (MQ=255)
                                                                                ggCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGcc   <  2:82410/88‑1 (MQ=255)
                                                                                ggCCTGGGAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGcc   >  1:82410/1‑88 (MQ=255)
                                                                                       |                                                                                 
AATCGCAAGGTATCTTCCGTTAGCGGTCGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCT  >  NZ_CP009273/601899‑602067

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