Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,499,566 A→G S258P (TCC→CCC)  ydcO ← BenE family transporter YdcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,499,5660AG96.4% 92.8 / ‑2.7 28coding (772/1176 nt)ydcOBenE family transporter YdcO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/14);  minor base . (1/0);  total (14/14)
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.89e-01

TTATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGAAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGAT  >  NZ_CP009273/1499487‑1499647
                                                                               |                                                                                 
ttATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGt                                                                         >  2:321598/1‑90 (MQ=255)
 tATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTg                                                                        <  2:674270/90‑1 (MQ=255)
 tATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTg                                                                        <  2:444086/90‑1 (MQ=255)
 tATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTg                                                                        <  1:104244/90‑1 (MQ=255)
   tCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGcc                                                                      <  2:663688/90‑1 (MQ=255)
       gATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCa                                                                  >  2:284689/1‑90 (MQ=255)
       gATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCa                                                                  >  2:270424/1‑90 (MQ=255)
       gATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCa                                                                  >  1:648694/1‑90 (MQ=255)
       gATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCa                                                                  >  1:488980/1‑90 (MQ=255)
          gcgcTTCCGGGTTTTGGCAAAAAGCCGCGGTGATTGCCGCAATCCGGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                               <  1:12696/90‑1 (MQ=255)
          gcgcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                               >  2:105040/1‑90 (MQ=255)
          gcgcTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATc                                                               <  1:105040/90‑1 (MQ=255)
                ccGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCAACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTaa                                                         >  2:12696/1‑90 (MQ=255)
                                   cgcgGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAAc                                      <  2:198572/90‑1 (MQ=255)
                                        tGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAg                                 <  1:321598/90‑1 (MQ=255)
                                                     aCCGACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgc                    <  1:169137/90‑1 (MQ=255)
                                                        gACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgctgc                 >  1:113504/1‑90 (MQ=255)
                                                        gACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgctgc                 >  1:279460/1‑90 (MQ=255)
                                                        gACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgc                    <  2:180958/87‑1 (MQ=255)
                                                        gACGGAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAgc                    >  1:180958/1‑87 (MQ=255)
                                                            gAATAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTc             <  2:488980/90‑1 (MQ=255)
                                                             aaTAAACGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCa            >  1:340724/1‑90 (MQ=255)
                                                                 aaCGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATtgc        <  2:279460/90‑1 (MQ=255)
                                                                 aaCGCCGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATtgc        <  2:648694/90‑1 (MQ=255)
                                                                      cGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATg                                             <  2:665276/48‑1 (MQ=255)
                                                                      cGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATg                                             >  1:665276/1‑48 (MQ=255)
                                                                      cGAAAGGGGGAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGa   >  1:747555/1‑90 (MQ=255)
                                                                      cGAAAGGGGGAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGAt  >  1:224521/1‑90 (MQ=255)
                                                                               |                                                                                 
TTATCCGGATGCGCTTCCGGGCTTTGGCAAATAGCCGCGGTGATTGCCGCAATACCGACGGAATAAACGCCGAAAGGGGAAAAAACCAGTGCCAGCAATCCAGTAAATACAATTAATGGCGAAACAGGAGCCGAATATCCAGCTGCTTTCATTGCTGCGAT  >  NZ_CP009273/1499487‑1499647

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