Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,090,967 G→A 9.3% V340M (GTG→ATG)  gudP → glucarate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,090,9670GA 9.3% 114.9 / 2.4 43V340M (GTG→ATG) gudPglucarate transporter
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base A (2/2);  total (20/23)
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.79e-01

CGGCGGCGTCCTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACATCGCGTCGATCAC  >  NC_002947/4090822‑4091109
                                                                                                                                                 |                                                                                                                                              
cggcggCGTCCTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTAGATGACGAGTGGAtggtgg                                                                                                                                             >  3:203519/1‑149 (MQ=255)
 ggcggcGTCCTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGAtggtgg                                                                                                                                             <  1:113406/148‑1 (MQ=255)
  gcggcgTCCTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGtt                                                                                                                                           >  6:74135/1‑149 (MQ=255)
          cTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCAt                                                                                                                                   >  6:327778/1‑149 (MQ=255)
           tGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCAt                                                                                                                                   >  1:144251/1‑148 (MQ=255)
                  tgtgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTgg                                                                                                                           >  6:134611/1‑149 (MQ=255)
                  tgtgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGATGGTTGGCtt                                                                                                                                      >  6:200930/1‑138 (MQ=255)
                  tgtgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGATGGTTGGCtt                                                                                                                                      <  5:200930/138‑1 (MQ=255)
                  tgtgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGATGGTTGGCtt                                                                                                                                      >  6:196364/1‑138 (MQ=255)
                  tgtgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGATGGTTGGCtt                                                                                                                                      <  5:196364/138‑1 (MQ=255)
                    tgATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGc                                                                                                                          <  6:210934/148‑1 (MQ=255)
                      aTCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGAtggtg                                                                                                                                              <  7:63906/126‑1 (MQ=255)
                      aTCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGAtggtg                                                                                                                                              >  8:63906/1‑126 (MQ=255)
                           cGACTGGCTTCTGCGTCGCGGCAACTGGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGttcttc                                                                                                                   <  4:203519/148‑1 (MQ=255)
                           cGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCg                                                                                                                  <  3:313282/149‑1 (MQ=255)
                              cTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCa                                                                                                               <  5:74135/149‑1 (MQ=255)
                                 gctgctGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAgg                                                                                                            >  3:219514/1‑149 (MQ=255)
                                 gctgctGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAgg                                                                                                            >  8:305131/1‑149 (MQ=255)
                                   tgctgcGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAggg                                                                                                           >  7:112524/1‑148 (MQ=255)
                                          ccgcGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGgcg                                                                                                   <  1:285049/149‑1 (MQ=255)
                                              ggCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCg                                                                                                                  <  7:96655/130‑1 (MQ=255)
                                              ggCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCg                                                                                                                  >  8:96655/1‑130 (MQ=255)
                                                               ctcGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTcg                                                                              >  2:326431/1‑149 (MQ=255)
                                                                  gcgcAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTcgcg                                                                            <  1:6477/148‑1 (MQ=255)
                                                                    gcAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGa                                                                          >  2:136306/1‑148 (MQ=255)
                                                                      aaaCTGCCGATTGTCTGTCGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGAGGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGTGCGGTGGTCGCGGacac                                                                       <  2:24161/149‑1 (MQ=255)
                                                                          tGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACtgggctgggc                                                                                      <  6:206928/130‑1 (MQ=255)
                                                                          tGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACtgggctgggc                                                                                      >  5:206928/1‑130 (MQ=255)
                                                                             cGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCatc                                                                 <  7:305131/148‑3 (MQ=255)
                                                                               aTTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCaa                                                              <  2:74226/149‑1 (MQ=255)
                                                                                ttGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAg                                                             <  4:285015/149‑1 (MQ=255)
                                                                                      gtgGCCTGCTGCTGTCCACCACCATGGTCTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGAtt                                                       >  1:185183/1‑149 (MQ=255)
                                                                                                    tCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTTGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCgg                                         >  6:48390/1‑149 (MQ=255)
                                                                                                    tCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCgg                                         >  2:256155/1‑149 (MQ=255)
                                                                                                      caccaccaTGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGt                                        <  5:134611/148‑1 (MQ=255)
                                                                                                      caccaccaTGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGt                                        <  2:185183/148‑1 (MQ=255)
                                                                                                       accaccaTGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTg                                       >  8:241811/1‑148 (MQ=255)
                                                                                                        cgaccaTGGGTTTCTGCAACTACGTCGATGCCGAGTGGATGGTGTTTGGCTTCTTGACGCTGGCGTTCTTCGGCAAGGGCATTGGCTCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCGAATCAGATTGCCGGGCTGTCCGGTGGa                                     <  4:219514/147‑1 (MQ=255)
                                                                                                        ccaccaTGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGGTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGa                                     <  5:327778/149‑1 (MQ=255)
                                                                                                            caTGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGt                                 <  2:144251/149‑1 (MQ=255)
                                                                                                                   ttCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAaca                           <  1:136306/148‑1 (MQ=255)
                                                                                                                               gTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACAt              >  8:47155/1‑149 (MQ=255)
                                                                                                                               gTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTAGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACATTCGGCAACAt              >  7:235628/1‑149 (MQ=255)
                                                                                                                                        gAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACATCGCGTCGAt     <  7:241811/149‑1 (MQ=255)
                                                                                                                                         cgTGGAGTGTGGGTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGCTGGACTGTTCAACACCTTCGGCAACATCGCGTCGATc    <  8:112524/148‑1 (MQ=255)
                                                                                                                                           tGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACATCGCGTCGATca   <  1:256155/148‑1 (MQ=255)
                                                                                                                                            ggATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACATCGCGTCGATcac  >  6:114338/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGGCGGCGTCCTGGGTGGTGTGATCTCCGACTGGCTGCTGCGCCGCGGCAACTCGCTGACCTTCTCGCGCAAACTGCCGATTGTCTGTGGCCTGCTGCTGTCCACCACCATGGTGTTCTGCAACTACGTCGATGCCGAGTGGATGGTGGTTGGCTTCATGACGCTGGCGTTCTTCGGCAAGGGCATTGGCGCACTGGGCTGGGCGGTGGTCGCGGACACTTCACCCAAGCAGATTGCCGGGCTGTCCGGTGGACTGTTCAACACCTTCGGCAACATCGCGTCGATCAC  >  NC_002947/4090822‑4091109

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