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,4300AG88.2% 42.8 / ‑6.2 17K391E (AAA→GAA) ydcRPLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (8/7);  total (8/9)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGC  >  NZ_CP009273/1505358‑1505509
                                                                        |                                                                               
cTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGAc                                                                >  1:307351/1‑90 (MQ=255)
  gCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAg                                                              >  2:974052/1‑90 (MQ=255)
  gCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAg                                                              <  2:104316/90‑1 (MQ=255)
        ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTa                                                        <  2:988497/90‑1 (MQ=255)
              tGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtct                                                  <  1:681419/90‑1 (MQ=255)
               gCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtctc                                                 <  2:401551/90‑1 (MQ=255)
                      gTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTgg                                          >  2:361625/1‑90 (MQ=255)
                              cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTcccc                                  >  1:432929/1‑90 (MQ=255)
                                tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGa                                <  1:279202/90‑1 (MQ=255)
                                     ggcaggcaCTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGt                           >  1:135781/1‑90 (MQ=255)
                                       caggcaCTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTa                         <  2:102134/90‑1 (MQ=255)
                                               gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                 >  2:216632/1‑90 (MQ=255)
                                                     ttATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAAtt           >  2:288776/1‑90 (MQ=255)
                                                            cctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTg    >  2:653626/1‑90 (MQ=255)
                                                              tgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc  <  2:135781/90‑1 (MQ=255)
                                                              tgcGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc  <  2:933984/90‑1 (MQ=255)
                                                              cggggAGGTGAAAATTCATCATAATGACAGGGTTTTCTTTTTCGGGTGGGAGCCCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc  <  1:653626/87‑1 (MQ=255)
                                                                        |                                                                               
CTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGC  >  NZ_CP009273/1505358‑1505509

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