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,5420AG91.7% 52.5 / ‑0.2 24K287E (AAG→GAG) ydgHDUF1471 family protein YdgH
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/14);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

CGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATA  >  NZ_CP009273/1673453‑1673627
                                                                                         |                                                                                     
cGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCga                                                                                       <  1:512104/90‑1 (MQ=255)
     atgGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGaaa                                                                                  <  1:813013/90‑1 (MQ=255)
           ccGTTCGATAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCt                                                                            <  1:376979/90‑1 (MQ=255)
           ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCt                                                                            <  2:437584/90‑1 (MQ=255)
            cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:636089/90‑1 (MQ=255)
            cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  2:280873/90‑1 (MQ=255)
            cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           >  2:288197/1‑90 (MQ=255)
            cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTa                                                                           <  1:268824/90‑1 (MQ=255)
             gTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTaa                                                                          >  2:175466/1‑90 (MQ=255)
               tCGCCAGTATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGt                                                                        <  2:781686/90‑1 (MQ=255)
               tCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGt                                                                        <  2:276985/90‑1 (MQ=255)
                              tctcTGGAAACTCCGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCcccc                                                         <  2:188804/90‑1 (MQ=255)
                                 ctGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGcc                                                      >  1:469030/1‑90 (MQ=255)
                                             gCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAAc                                          <  1:288197/90‑1 (MQ=255)
                                               aaCATGACCGAAGTCTCCTATCAGGTTGCGAACCGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGt                                        <  1:175466/90‑1 (MQ=255)
                                               aaCATGACCGAAGTCCCCTATCAGGTTGCGAAACGTCCCGCGGAAAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGt                                        <  1:596917/90‑1 (MQ=255)
                                                     aCCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGtaat                                  >  1:518116/1‑90 (MQ=255)
                                                            tctcCTATAAGGTTGCGAAGCGTCCCGGGGAGAAAGGTGATAAATTTTACCACACCGCCCGCCAGTGGCAGGAACGTGGTAATAACCTGa                           <  1:446833/90‑1 (MQ=255)
                                                                      ggTTGCAAACCGTGCGGGGGAGAAGGGTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                 <  2:13053/90‑1 (MQ=255)
                                                                       gTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCa                >  1:566161/1‑90 (MQ=255)
                                                                         tGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGa              >  1:204948/1‑90 (MQ=255)
                                                                         tGCGAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGGGGCAGGAAAGGGGGAATAACGTTTCCGGCCGCGCAGa              >  2:546883/1‑90 (MQ=255)
                                                                            gAAACGTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCt           >  2:246447/1‑90 (MQ=255)
                                                                                 gTGCCGCGGAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGGAATAACCTGACCGTCAGCGCAGATCTGtata      >  2:201364/1‑90 (MQ=255)
                                                                                     cgcgGAAAAGGGTGTTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaata  <  2:51633/90‑1 (MQ=255)
                                                                                         |                                                                                     
CGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATA  >  NZ_CP009273/1673453‑1673627

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