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,5420AG96.8% 82.9 / NA 31K287E (AAG→GAG) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (6/24);  total (6/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGAT  >  NZ_CP009273/1673458‑1673616
                                                                                    |                                                                          
atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       <  1:1071814/90‑1 (MQ=255)
atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       <  2:1512221/90‑1 (MQ=255)
atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       <  2:1045479/90‑1 (MQ=255)
atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       <  1:585036/90‑1 (MQ=255)
  ggTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAgg                                                                     <  1:394500/90‑1 (MQ=255)
      ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                 <  1:846391/90‑1 (MQ=255)
      ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                 <  2:830181/90‑1 (MQ=255)
      ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                 <  2:1104588/90‑1 (MQ=255)
       cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                <  1:5984/90‑1 (MQ=255)
        gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                               <  2:920124/90‑1 (MQ=255)
        gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                               <  2:1319893/90‑1 (MQ=255)
        gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                >  2:257635/1‑89 (MQ=255)
        gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                <  1:257635/89‑1 (MQ=255)
             aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                 <  1:189110/83‑1 (MQ=255)
             aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                 >  2:189110/1‑83 (MQ=255)
             aCAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTAtt                                                          >  1:1057792/1‑90 (MQ=255)
                          ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       <  1:1364798/64‑1 (MQ=255)
                          ctctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                       >  2:1364798/1‑64 (MQ=255)
                              ggCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCAccagccag                                         <  1:762506/90‑1 (MQ=255)
                                      cgacAACACGCCCAAATTCTCCTATAAGGTGGCGAAACGGGCGGGGGAGAAGGGTGATAATTATTCCCACATCCCCCGCCAGTGGCAGGa                                 <  1:965054/87‑1 (MQ=255)
                                          aaCATGCCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAAAAAGGTGCTAAGTATTCCAACACCACCCGCCAGTGGCAGGAACGt                             <  1:310168/90‑1 (MQ=255)
                                          aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGGGGAAAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGt                             <  1:450912/90‑1 (MQ=255)
                                          aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGt                             <  2:557193/90‑1 (MQ=255)
                                          aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGGGCCGCGGAGAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGt                             <  1:275432/90‑1 (MQ=255)
                                          aaCATCCCCGAACTCTCCTATCGGGTTGCCAAACGTGCCGCGGAGAAAGATGTTAATTATTACCAAATCACCCGCCAGTGGCAGGAACGt                             <  2:718398/90‑1 (MQ=255)
                                              tGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGta                         <  2:1057792/90‑1 (MQ=255)
                                                aCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaat                       >  1:1549037/1‑90 (MQ=255)
                                                            tATCAGTTTGCGAAACGTGCCGGGGAAAAGGGTGTTAATTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAATCTGACCGTc           <  1:370637/90‑1 (MQ=255)
                                                                aGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAgcg       >  1:86462/1‑90 (MQ=255)
                                                                 gtttGCAAAACGTGCCGGGGAGAAAGGTGATAATTATTACCCCATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc      <  2:728064/88‑1 (MQ=255)
                                                                     gagAAAGGTGCGGAGAAGAAGGGTTTAAAATTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGAt  <  1:158780/88‑1 (MQ=255)
                                                                                    |                                                                          
ATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGAT  >  NZ_CP009273/1673458‑1673616

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