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,0100TC98.4% 200.8 / ‑3.7 62K2E (AAA→GAA) kdgAbifunctional 4‑hydroxy‑2‑oxoglutarate aldolase/2‑dehydro‑3‑deoxy‑phosphogluconate aldolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/41);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

GCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACA  >  NZ_CP009273/1926923‑1927090
                                                                                       |                                                                                
gCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTgttttccagttttcca                                                                                >  1:735882/1‑90 (MQ=255)
   ccGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTccatca                                                                             <  2:1564133/90‑1 (MQ=255)
    cgcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCag                                                                            <  2:1007300/90‑1 (MQ=255)
     gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                           >  2:1292124/1‑90 (MQ=255)
     gcgTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCaga                                                                           >  2:1290288/1‑90 (MQ=255)
               gTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTtctctc                                                                 <  2:1050914/90‑1 (MQ=255)
                  tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                              <  2:1429019/90‑1 (MQ=255)
                  tttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGcc                                                              <  1:868030/90‑1 (MQ=255)
                   ttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCt                                                             >  2:280797/1‑90 (MQ=255)
                   ttttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGACTGATTCTGCACTTGTTTTCCAGTTTCCCATCAGAGTTTTCTCTCGCCt                                                             <  1:984623/90‑1 (MQ=255)
                    tttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTg                                                            <  2:1565537/90‑1 (MQ=255)
                    tttACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTg                                                            <  1:59798/90‑1 (MQ=255)
                         cACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTAc                                                       <  2:192629/90‑1 (MQ=255)
                         cACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTAc                                                       >  1:992453/1‑90 (MQ=255)
                            gATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACaaa                                                    >  2:661259/1‑90 (MQ=255)
                                   ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                             <  1:1555274/90‑1 (MQ=255)
                                   ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                             <  1:1543877/90‑1 (MQ=255)
                                   ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                             <  2:503229/90‑1 (MQ=255)
                                   ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                             <  1:1305487/90‑1 (MQ=255)
                                   ggTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcg                                             <  2:992453/90‑1 (MQ=255)
                                     tACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcgtc                                           >  2:340349/1‑90 (MQ=255)
                                     tACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTgtcgtc                                           >  1:83771/1‑90 (MQ=255)
                                         aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                       >  2:175308/1‑90 (MQ=255)
                                         aCCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                       <  2:57955/90‑1 (MQ=255)
                                           cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                     <  1:369646/90‑1 (MQ=255)
                                           cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                     <  1:1406178/90‑1 (MQ=255)
                                           cGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaa                                     <  2:1172378/90‑1 (MQ=255)
                                            gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                                    >  1:1027006/1‑90 (MQ=255)
                                            gggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaaaaa                                    >  1:1083280/1‑90 (MQ=255)
                                             ggCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                                   <  2:963466/90‑1 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCCCTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTaa                                       <  1:824497/85‑1 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGGTTACAAATTTGTCGTCTTaa                                       >  2:824497/1‑85 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                                   >  2:1081968/1‑89 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAg                                   <  1:1081968/89‑1 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                                  >  2:1344131/1‑90 (MQ=255)
                                              gCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGt                                  >  2:457328/1‑90 (MQ=255)
                                                   tggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATAc                             <  2:439911/90‑1 (MQ=255)
                                                   tggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATAc                             <  2:1524668/90‑1 (MQ=255)
                                                    ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                            <  2:83771/90‑1 (MQ=255)
                                                    ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                            <  2:84717/90‑1 (MQ=255)
                                                    ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                            <  2:591834/90‑1 (MQ=255)
                                                    ggtCAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCCCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACa                            <  1:1005687/90‑1 (MQ=255)
                                                       cAGGATTGATTCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGt                         >  2:363680/1‑90 (MQ=255)
                                                            ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc                    >  1:1201665/1‑90 (MQ=255)
                                                            ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc                    <  1:1344131/90‑1 (MQ=255)
                                                            ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc                    <  1:306144/90‑1 (MQ=255)
                                                            ttgattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTgcgc                    <  1:68352/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  2:824989/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  2:1183372/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  2:356919/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  2:1424722/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  1:374876/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  1:30684/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  1:24281/90‑1 (MQ=255)
                                                              gattCTGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGccc                  <  1:1292124/90‑1 (MQ=255)
                                                                   tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc             >  2:244176/1‑90 (MQ=255)
                                                                   tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc             <  2:1418251/90‑1 (MQ=255)
                                                                   tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc             >  1:591535/1‑90 (MQ=255)
                                                                   tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc             >  1:170243/1‑90 (MQ=255)
                                                                   tGCACTTGTTTTCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTc             >  1:100585/1‑90 (MQ=255)
                                                                       cTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCa         >  2:191372/1‑90 (MQ=255)
                                                                             ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc   <  2:1201665/90‑1 (MQ=255)
                                                                             ttCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGAc   <  2:645705/90‑1 (MQ=255)
                                                                              tCCAGTTTTCCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACa  <  1:61518/90‑1 (MQ=255)
                                                                                       |                                                                                
GCACCGCGTGTTCCAGTTTTTTTACCACGATAACCGGTACAACCGGGCCGGTGGTCAGGATTGATTCTGCACTTGTTTTCCAGTTTTTCATCAGAGTTTTCTCTCGCCTGATTACAAATTTGTCGTCTTAAAAAGTGATACAGGTTGCGCCCTGTTCGGCACCGGACA  >  NZ_CP009273/1926923‑1927090

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