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

CCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCG  >  NZ_CP009273/1505345‑1505510
                                                                                     |                                                                                
ccATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAAt                                                                              <  2:99300/90‑1 (MQ=255)
 cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAAtt                                                                             <  1:556265/90‑1 (MQ=255)
 cATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAAtt                                                                             >  1:6633/1‑90 (MQ=255)
   tCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtca                                                                           <  1:511776/90‑1 (MQ=255)
   tCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtca                                                                           >  2:361648/1‑90 (MQ=255)
      tcgtcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATtcatca                                                                        >  1:798477/1‑90 (MQ=255)
        gtcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATa                                                                      <  2:6633/90‑1 (MQ=255)
         tcgCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATaa                                                                     >  2:54916/1‑90 (MQ=255)
             cTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGAc                                                                 >  2:734858/1‑90 (MQ=255)
             cTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGAc                                                                 >  2:560325/1‑90 (MQ=255)
                     ccAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTa                                                         <  2:43069/90‑1 (MQ=255)
                                      aacaacGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGc                                        <  2:26414/90‑1 (MQ=255)
                                           cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtctc                                                  <  1:651495/75‑1 (MQ=255)
                                           cGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTtctc                                                  >  2:651495/1‑75 (MQ=255)
                                            gTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCg                                  >  1:917741/1‑90 (MQ=255)
                                             tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTATCATAATGACa                                                                >  1:905027/1‑58 (MQ=255)
                                             tGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGGAAATTATCATAATGACa                                                                <  2:905027/58‑1 (MQ=255)
                                                            gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                  >  1:523715/1‑90 (MQ=255)
                                                            gctgcGTTATCTGCCTGCGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCgg                  >  2:910953/1‑90 (MQ=255)
                                                                        gcctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCt      >  2:655429/1‑90 (MQ=255)
                                                                         cctgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTg     >  2:30580/1‑90 (MQ=255)
                                                                           tgcGGAGGTGAAAATTCATAATAATGAAGGTGTTTTCTTTCTCGGGTGGGAGCCCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc   <  2:156132/90‑1 (MQ=255)
                                                                           tgcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggc   <  1:30580/90‑1 (MQ=255)
                                                                            gcGGAAGTGGAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTggcg  >  2:35332/1‑90 (MQ=255)
                                                                                     |                                                                                
CCATCTTCGTCGCCTGCGTCGCCAGCTTGCGGAACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCG  >  NZ_CP009273/1505345‑1505510

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