Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,506,028 G→A A40T (GCC→ACC)  ydcS → ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,506,0280GA92.6% 99.6 / NA 27A40T (GCC→ACC) ydcSABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/11);  minor base C (0/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 2.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGT  >  NZ_CP009273/1505893‑1506153
                                                                                                                                       |                                                                                                                             
cGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCt                                                                                                                            <  1:492278/139‑1 (MQ=255)
  aCAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTgg                                                                                                                          <  2:315750/139‑1 (MQ=255)
                       cAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGa                                                                                                     >  2:467552/1‑139 (MQ=255)
                                       agcagcCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACaa                                                                                      >  1:620284/1‑138 (MQ=255)
                                           gcCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 <  1:420951/139‑1 (MQ=255)
                                            ccTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGAc                                                                                >  1:743501/1‑139 (MQ=255)
                                                     gcTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACg                                                                       <  2:159313/139‑1 (MQ=255)
                                                         aGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGt                                                                   >  2:620948/1‑139 (MQ=255)
                                                          gCATGACAATAATGACCGCTCACGCCGCCGAATCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGtt                                                                  >  2:439010/1‑139 (MQ=255)
                                                                cAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGaaa                                                             >  2:423072/1‑138 (MQ=255)
                                                                   taatGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 >  2:263722/1‑115 (MQ=255)
                                                                   taatGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 <  1:263722/115‑1 (MQ=255)
                                                                       gACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 >  1:82007/1‑111 (MQ=255)
                                                                       gACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 <  2:82007/111‑1 (MQ=255)
                                                                          cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 >  1:547437/1‑108 (MQ=255)
                                                                          cGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGa                                                                                 <  2:547437/108‑1 (MQ=255)
                                                                           gCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACa                                                                                                   >  2:833224/1‑89 (MQ=255)
                                                                           gCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACa                                                                                                   <  1:833224/89‑1 (MQ=255)
                                                                              cACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcg                                              >  1:546203/1‑139 (MQ=255)
                                                                               aCGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTgcg                                              <  2:145027/138‑1 (MQ=255)
                                                                                 gccgccGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCgg                                           <  1:467552/139‑1 (MQ=255)
                                                                                        aaCCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAAtgt                                    >  2:607647/1‑139 (MQ=255)
                                                                                                              aCCGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGAt              >  2:264663/1‑139 (MQ=255)
                                                                                                               ccGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACa                                                                                                   <  2:165876/53‑1 (MQ=255)
                                                                                                               ccGGAAGGGCGACTGGATATTATCACCTGGCCGGGATACATCGAACGCGGACa                                                                                                   >  1:165876/1‑53 (MQ=255)
                                                                                                                       gCGACTGGATATTATCCCCTGGCCGGGATACATCGAACGCGGACAAACTGAAAAACAATACGACGGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTc     <  2:422935/139‑1 (MQ=255)
                                                                                                                          agccagtattctCCCCTGGCCGGGATACAGCGAACGCGGACAAACTGATAAACAATCCGCCTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGTTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGt  <  2:770348/133‑1 (MQ=255)
                                                                                                                                       |                                                                                                                             
CGACAGGAGTAAGACCTTATGAGCAAGACATTTGCCCGCAGCAGCCTGTGTGCGCTCAGCATGACAATAATGACCGCTCACGCCGCCGAACCGCCTACCAATTTAGATAAACCGGAAGGGCGACTGGATATTATCGCCTGGCCGGGATACATCGAACGCGGACAAACTGATAAACAATACGACTGGGTAACGCAGTTCGAAAAAGAGACAGGCTGCGCGGTGAATGTGAAAACCGCCGCGACTTCCGATGAAATGGTCAGT  >  NZ_CP009273/1505893‑1506153

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: