Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,505,430 A→G K391E (AAA→GAA)  ydcR → PLP‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,505,4300AG89.5% 40.5 / ‑0.2 19K391E (AAA→GAA) ydcRPLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (8/9);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

TTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCAC  >  NZ_CP009273/1505350‑1505514
                                                                                |                                                                                    
ttCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtcatc                                                                             <  2:58259/90‑1 (MQ=255)
    tcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATaa                                                                         >  2:464826/1‑90 (MQ=255)
       ccTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGa                                                                      <  2:601163/90‑1 (MQ=255)
               gCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGtt                                                              <  1:704609/90‑1 (MQ=255)
                       gCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtctc                                                      <  1:775860/90‑1 (MQ=255)
                       gCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtctc                                                      <  2:429305/90‑1 (MQ=255)
                                      cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGTGCTcccc                                       >  1:13222/1‑90 (MQ=255)
                                        tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGa                                     >  1:666145/1‑90 (MQ=255)
                                         gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtca                                                                               >  1:22122/1‑47 (MQ=255)
                                         gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtca                                                                               <  2:22122/47‑1 (MQ=255)
                                         gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTcatcat                                                                           <  2:362989/51‑1 (MQ=255)
                                         gCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTcatcat                                                                           >  1:362989/1‑51 (MQ=255)
                                                   caCTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATg                          >  2:778484/1‑90 (MQ=255)
                                                        ctgcTTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCCCCCCGAGCCGTTAGATGCCGGc                     <  2:98458/90‑1 (MQ=255)
                                                        ctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGTTTCCTTTCTTTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGc                     <  1:485325/90‑1 (MQ=255)
                                                                 ctgcctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGc            >  2:61353/1‑90 (MQ=255)
                                                                  tgcctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGcc           >  2:485325/1‑90 (MQ=255)
                                                                    cctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTg         >  1:32245/1‑90 (MQ=255)
                                                                           aaGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCAc  >  2:828490/1‑90 (MQ=255)
                                                                                |                                                                                    
TTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCAC  >  NZ_CP009273/1505350‑1505514

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