Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,927,010 T→C K2E (AAA→GAA)  kdgA ← bifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,927,0100TC94.4% 48.4 / ‑3.5 18K2E (AAA→GAA) kdgAbifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (8/9);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

CCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGTTTT  >  NZ_CP009273/1926926‑1927095
                                                                                    |                                                                                     
ccGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                                  >  2:298773/1‑90 (MQ=255)
  gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                                >  1:332897/1‑90 (MQ=255)
  gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                                >  2:257553/1‑90 (MQ=255)
     tgtTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGtt                                                                             <  1:499624/90‑1 (MQ=255)
            gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTtctctc                                                                      >  2:120219/1‑90 (MQ=255)
               tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                                   <  2:554423/90‑1 (MQ=255)
               tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                                   <  1:474270/90‑1 (MQ=255)
                                ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                                  <  1:353418/90‑1 (MQ=255)
                                ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                                  <  2:298243/90‑1 (MQ=255)
                                      aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                            >  2:383697/1‑90 (MQ=255)
                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                          <  1:427571/90‑1 (MQ=255)
                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                          <  2:450493/90‑1 (MQ=255)
                                        cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                          <  2:48139/90‑1 (MQ=255)
                                           gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                                       >  2:283857/1‑90 (MQ=255)
                                                           gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                       >  1:304383/1‑90 (MQ=255)
                                                                     ttGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCAc             >  1:406012/1‑90 (MQ=255)
                                                                          ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc        <  1:383697/90‑1 (MQ=255)
                                                                          ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc        <  2:412932/90‑1 (MQ=255)
                                                                              aGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGtt    >  1:31783/1‑90 (MQ=255)
                                                                                ttttcCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGtttt  <  1:427832/90‑1 (MQ=255)
                                                                                ttttcCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGtttt  >  2:202722/1‑90 (MQ=255)
                                                                                    |                                                                                     
CCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACAGTTTT  >  NZ_CP009273/1926926‑1927095

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