Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,323,088 G→A 9.5% A171V (GCC→GTC)  kdgA ← 2‑keto‑3‑deoxygluconate‑6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,323,0880GA 9.5% 140.6 / 3.5 47A171V (GCC→GTC) kdgA2‑keto‑3‑deoxygluconate‑6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (20/22);  new base A (2/2);  total (23/24)
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.97e-01

GCCTTGATTTTCATATTGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCTTGGGAAC  >  BS168_reference_genome/2322951‑2323207
                                                                                                                                         |                                                                                                                       
gCCTTGATTTTCATATTGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTTAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGc                                                                                                                        >  5:105073/1‑139 (MQ=255)
gCCTTGATTTTCATATTGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGc                                                                                                                        >  7:3539/1‑139 (MQ=255)
       ttttCATATTGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGc                                                                                                                 <  7:35288/139‑1 (MQ=255)
               ttGTGGATTCTCCCTTCCTCCTGTTACGAAAAAACAATTCCACTTACACTTGGACAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAgg                                                                                                         <  5:21425/139‑1 (MQ=255)
               ttGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAgg                                                                                                         <  6:105073/139‑1 (MQ=255)
                tgtgGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAgg                                                                                                         <  1:61521/138‑1 (MQ=255)
                tgtgGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAgg                                                                                                         <  4:26537/138‑1 (MQ=255)
                       tctcCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACtt                                                                                                    >  2:96187/1‑136 (MQ=255)
                        ctcCCATCCTCCGGTTATCAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCGGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGa                                                                                                <  6:159338/139‑1 (MQ=255)
                           ccATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGa                                                                                                <  1:96187/136‑1 (MQ=255)
                              tcctccGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGCCGGCGCCCGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATg                                                                                          >  1:290015/1‑139 (MQ=255)
                                ctccGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTa                                                                                        >  4:103066/1‑139 (MQ=255)
                                 tccGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATc                                                                                                            >  7:183550/1‑118 (MQ=255)
                                 tccGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATc                                                                                                            <  8:183550/118‑1 (MQ=255)
                                 tccGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCAAACTGGATGACGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGTCACTACAGACGGATGTa                                                                                        >  7:272577/1‑138 (MQ=255)
                                   cGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTCCCCAACTGGCTGCCTCCACCGACGGCGCCAGCTCCGCCTCTAAGCCAATCAGGCACTTCAGACGGAAGTATcc                                                                                     >  6:192563/1‑139 (MQ=255)
                                           aaaaaaCAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGt                                                                              <  3:103066/138‑1 (MQ=255)
                                               aaCAATTCCACTAACATTTGGAAACCACCCTGCAAATCCTCTTTTCAACAGCTGCCCACCTGGCGGCCGACTCCGACGGCGCAAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGaa                                                                         <  2:277883/139‑1 (MQ=255)
                                                    ttCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGa                                                                                            >  8:75478/1‑115 (MQ=255)
                                                    ttCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGa                                                                                            <  7:75478/115‑1 (MQ=255)
                                                                  ggAAAACAGCCTGCAAATGCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCGAGCTCCGGCTATTAGCCTATCAGGCACTTCAGACGGATGTATCCCGCCTGTCCGAATAAAGCTCACCTGCggcc                                                       >  7:304057/1‑136 (MQ=255)
                                                                  ggAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGa                                                       >  3:158154/1‑138 (MQ=255)
                                                                          gCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCg                                              >  7:58631/1‑139 (MQ=255)
                                                                             tGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCt                                            >  7:44008/1‑138 (MQ=255)
                                                                                              aCAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACgg                           >  7:85474/1‑138 (MQ=255)
                                                                                                aGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAAt                        <  8:46753/139‑1 (MQ=255)
                                                                                                aGCTGCCCAACTGGCTGCCGACGCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAAt                        <  6:164491/139‑1 (MQ=255)
                                                                                                  cTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCAt                                 <  2:65595/128‑1 (MQ=255)
                                                                                                  cTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGc                      >  7:237457/1‑139 (MQ=255)
                                                                                                  cTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGc                      >  1:31241/1‑139 (MQ=255)
                                                                                                  cTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTACAt                                 >  1:65595/1‑128 (MQ=255)
                                                                                                    gCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGaa                                                                         >  8:204649/1‑86 (MQ=255)
                                                                                                    gCCCAACTGGCGGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGaa                                                                         <  7:204649/86‑1 (MQ=255)
                                                                                                        aagtGGCTGCCGCCTCCGACGGCGCGAGCTCAGGCTCTAAGCCTATCAGGCACTTCAGACGGATGTATCCGGCCTGTCGGAATAAAGGTCACCTGCGGGAATGGACCCGCTAAATTTTTCATAAACGGAATGCCAAaca                <  5:192563/136‑1 (MQ=255)
                                                                                                        aaCTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAaca                <  8:237457/139‑1 (MQ=255)
                                                                                                          cTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGa                                                                                            >  7:266886/1‑61 (MQ=255)
                                                                                                          cTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGa                                                                                            <  8:266886/61‑1 (MQ=255)
                                                                                                            ggCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAAGTTTTCATAAACGGAATGCCAAACACAcc            >  8:91035/1‑139 (MQ=255)
                                                                                                             gCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAacacac             <  8:44008/137‑1 (MQ=255)
                                                                                                               tGCCGACCCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCt         >  1:132474/1‑139 (MQ=255)
                                                                                                                     tgccGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCTTGGGaa   <  2:290015/137‑1 (MQ=255)
                                                                                                                     cTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCTTGGGaa   <  4:305201/139‑1 (MQ=255)
                                                                                                                      accgccgGCGCCAGCTCCGGCGCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCTTGGGAAc  <  4:42096/138‑1 (MQ=255)
                                                                                                                              cgcCAGCTCCGACTCTAAGCCAATCAGGCACTTCa                                                                                                  <  4:131669/35‑1 (MQ=255)
                                                                                                                              cgcCAGCTCCGACTCTAAGCCAATCAGGCACTTCa                                                                                                  <  4:126986/35‑1 (MQ=255)
                                                                                                                              cgcCAGCTCCGACTCTAAGCCAATCAGGCACTTCa                                                                                                  >  3:131669/1‑35 (MQ=255)
                                                                                                                              cgcCAGCTCCGACTCTAAGCCAATCAGGCACTTCa                                                                                                  >  3:126986/1‑35 (MQ=255)
                                                                                                                                         |                                                                                                                       
GCCTTGATTTTCATATTGTGGATTCTCCCATCCTCCGGTTATGAAAAAACAATTCCACTTACACTTGGAAAACAGCCTGCAAATCCTCTTTTGAACAGCTGCCCAACTGGCTGCCGACTCCGACGGCGCCAGCTCCGGCTCTAAGCCAATCAGGCACTTCAGACGGATGTATCCCGCCTGTCGGAATAAAGGTCACCTGCGGGAAAGGACCCGCTAAATTTTTCATAAACGGAATGCCAAACACACCGCTTGGGAAC  >  BS168_reference_genome/2322951‑2323207

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