Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,417,434 C→G 17.9% A100P (GCG→CCG)  gudD ← bifunctional D‑glucarate dehydratase/L‑idarate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,417,4340CG17.9% 65.3 / 9.0 28A100P (GCG→CCG) gudDbifunctional D‑glucarate dehydratase/L‑idarate epimerase
Reads supporting (aligned to +/- strand):  ref base C (10/13);  new base G (2/3);  total (12/16)
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.40e-01

CGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTA  >  NC_002947/5417296‑5417576
                                                                                                                                          |                                                                                                                                              
cgcgTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGt                                                                                                                                     <  2:545691/150‑1 (MQ=255)
cgcgTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGAGGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTTTGCAGGCCGCGACTCGGCGGGTCGCggc                                                                                                                                     <  1:228241/150‑2 (MQ=255)
          gCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGt                                                                                                                           >  1:739552/1‑150 (MQ=255)
            ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACTCGGCGGGTCGCggc                                                                                                                                     <  1:222312/138‑2 (MQ=255)
            ccTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACTCGGCGGGTCGCggc                                                                                                                                     >  2:222312/1‑137 (MQ=255)
                      cAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGc                                                                                                               >  1:600323/1‑150 (MQ=255)
                      cAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCCAGTCGCGGTTGGCGAAGGTCTGGCGCATGGc                                                                                                               <  2:527313/150‑1 (MQ=255)
                                           cAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATg                                                                                          <  1:820126/150‑1 (MQ=255)
                                             ggTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGc                                                                                        <  1:133707/150‑1 (MQ=255)
                                                       gAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTg                                                                              <  2:409414/150‑1 (MQ=255)
                                                         gCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGc                                                                            <  1:420153/150‑1 (MQ=255)
                                                            ggTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGa                                                                         >  1:454001/1‑150 (MQ=255)
                                                                 aGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGATCACGCGCTGGTAATgg                                                                                         <  2:549893/129‑1 (MQ=255)
                                                                 aGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGATCACGCGCTGGTAATgg                                                                                         >  1:549893/1‑129 (MQ=255)
                                                                          gccTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGCCCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                           <  2:352777/148‑1 (MQ=255)
                                                                          gACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCgg                                                           <  2:454001/150‑1 (MQ=255)
                                                                                      gTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAgg                                               <  1:678208/150‑1 (MQ=255)
                                                                                             cATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGc                                        >  2:166093/1‑150 (MQ=255)
                                                                                                gTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGaa                                     >  2:467812/1‑150 (MQ=255)
                                                                                                    ggTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGTCGCATCGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtctt                                 >  2:761484/1‑150 (MQ=255)
                                                                                                                gTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                               <  2:158325/140‑1 (MQ=255)
                                                                                                                gTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                               >  1:158325/1‑140 (MQ=255)
                                                                                                                            cAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAg         >  2:572080/1‑150 (MQ=255)
                                                                                                                               gCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGcc      >  1:306916/1‑150 (MQ=255)
                                                                                                                                 cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                               >  2:35329/1‑123 (MQ=255)
                                                                                                                                 cgcgACTCGGCGGGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAAtcttct                               <  1:35329/123‑1 (MQ=255)
                                                                                                                                 cgcgACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCgg    <  2:600323/150‑1 (MQ=255)
                                                                                                                                   cgACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTa  <  1:166093/150‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                              
CGCGTTGCTGGCCCTCGCCCAGCAATGCTGCCATCGGCACATTCAGGTGCTGGCCGAGCAGGTCGAGCAACGCCGACTCGATGGCCGTTACCGCATGTACGGTAATGCGCAGGTCGAAGGTCTGCAGGCCGCGACCCGCCGAGTCGCGGTTGGCGAAGGTCTGGCGCATGGCATTGAGCACGCGCTGGTAATGGCCGATCGGCTGGCCGACCACCAGGCTGCGGGCGTCCTCCAGGGTCTGGCGAATCTTCTCGCCGCCGGGCACTTCGCCCAGGCCGGTA  >  NC_002947/5417296‑5417576

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