Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,673,542 A→G K287E (AAG→GAG)  ydgH → DUF1471 family protein YdgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,673,5420AG90.9% 21.5 / ‑5.7 11K287E (AAG→GAG) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (4/6);  total (4/7)
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.82e-01

GATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATA  >  NZ_CP009273/1673457‑1673623
                                                                                     |                                                                                  
gatgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCggagaa                                                                                >  2:348536/1‑90 (MQ=255)
 atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                               <  1:45502/90‑1 (MQ=255)
  tgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAg                                                                              <  1:363029/90‑1 (MQ=255)
                                           aaCATGACCGAAGTCTCCTATCGGGTTGCAAAACGGGCCGGGGAAAAAGGTGC‑TAATTTTTACCCCATCCCCCGCCAGTGGCAGGAACGt                                    <  1:348536/90‑1 (MQ=255)
                                           aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAAAAAGGTGC‑TAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGt                                    <  2:147399/90‑1 (MQ=255)
                                            aCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTg                                   >  1:183042/1‑90 (MQ=255)
                                              aTGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTggg                                 >  1:46188/1‑89 (MQ=255)
                                                                  ggTTGCGAAACGGGCCGGGGAAAAGGGTGC‑TAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc             <  2:183042/90‑1 (MQ=255)
                                                                  ggTTGCAAACGGTGCGGGGAA‑AAAGGTTCTTAAGTTTTCCCCCTTCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc             <  1:337953/90‑1 (MQ=255)
                                                                           aCGTGCCGCGGAGAAAGGTGC‑TAAGTATTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGta    <  2:46188/90‑1 (MQ=255)
                                                                             gTGCCGCGGAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGtata  >  2:149950/1‑90 (MQ=255)
                                                                             gTGCCGCGGAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGATCTGtatc  >  2:337953/1‑89 (MQ=255)
                                                                                     |                                                                                  
GATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGC‑TAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATA  >  NZ_CP009273/1673457‑1673623

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