Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,546,774 A→G 9.7% L191P (CTG→CCG)  pydP ← NCS1 family transporter PydP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,546,7740AG 9.7% 106.3 / 3.4 41L191P (CTG→CCG) pydPNCS1 family transporter PydP
Reads supporting (aligned to +/- strand):  ref base A (19/18);  new base G (2/2);  total (21/20)
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.79e-01

AGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGAATGGCCAGCCC  >  NC_002947/4546627‑4546916
                                                                                                                                                   |                                                                                                                                              
aGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCa                                                                                                                                                <  8:228780/148‑1 (MQ=255)
    cGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCAc                                                                                                                                           >  1:285238/1‑149 (MQ=255)
    cGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCAc                                                                                                                                           <  6:267792/149‑1 (MQ=255)
    cGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCAc                                                                                                                                           <  2:285238/149‑1 (MQ=255)
            aaGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGa                                                                                                                                   <  2:169233/149‑1 (MQ=255)
                      ccACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGcgccg                                                                                                                         <  7:279640/149‑1 (MQ=255)
                       cACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGcgccgc                                                                                                                        >  4:237234/1‑149 (MQ=255)
                       cACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGcgccgc                                                                                                                        >  8:287499/1‑149 (MQ=255)
                         cGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCgcgc                                                                                                                      <  2:2617/149‑1 (MQ=255)
                           cTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCgcag                                                                                                                    <  7:287499/149‑1 (MQ=255)
                              gCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATc                                                                                                                                     <  6:1183/129‑1 (MQ=255)
                              gCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATc                                                                                                                                     >  5:1183/1‑129 (MQ=255)
                                cccTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGcacca                                                                                                               >  3:115334/1‑149 (MQ=255)
                                 ccTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGcaccac                                                                                                              <  4:5443/149‑1 (MQ=255)
                                       ggACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACAAGCAGCGGTGCCGACAGCGTTTGCAGCCACTGGATCGCCTCGGCGCCGCGCCGCACAACCCACAg                                                                                                        >  8:138205/1‑149 (MQ=255)
                                           gCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCGGCCACTTGATCGACTCGGCGCCGCgcagc                                                                                                                   >  8:156575/1‑134 (MQ=255)
                                           gCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCGGCCACTTGATCGACTCGGCGCCGCgcagc                                                                                                                   >  8:161048/1‑134 (MQ=255)
                                           gCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCGGCCACTTGATCGACTCGGCGCCGCgcagc                                                                                                                   <  7:161048/134‑1 (MQ=255)
                                           gCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCGGCCACTTGATCGACTCGGCGCCGCgcagc                                                                                                                   <  7:156575/134‑1 (MQ=255)
                                             ttGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAg                                                                                                  >  8:75509/1‑149 (MQ=255)
                                                   ggcTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCa                                                                                            <  5:138585/149‑1 (MQ=255)
                                                     cTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCag                                                                                           >  2:309041/1‑148 (MQ=255)
                                                           cAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCa                                                                                    >  2:124141/1‑149 (MQ=255)
                                                                  tCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGc                                                                              <  3:237234/148‑1 (MQ=255)
                                                                             cATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGc                                                                  <  1:124141/149‑1 (MQ=255)
                                                                                     gCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCg                                                          >  8:15092/1‑149 (MQ=255)
                                                                                                aaCAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCtt                                               >  2:77247/1‑149 (MQ=255)
                                                                                                        aCCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGt                                       >  6:297266/1‑149 (MQ=255)
                                                                                                         ccGACGCCCACTGCCACCAGCAGCGGTGCCGCCAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTc                                      >  2:231116/1‑149 (MQ=255)
                                                                                                                  aCTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCg                              >  5:135415/1‑148 (MQ=255)
                                                                                                                    tGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAAc                           >  2:257264/1‑149 (MQ=255)
                                                                                                                      ccaccaGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAAccc                         <  6:135415/149‑1 (MQ=255)
                                                                                                                       caccaGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCa                        <  1:309041/149‑1 (MQ=255)
                                                                                                                         ccaGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAg                       >  2:3625/1‑148 (MQ=255)
                                                                                                                            gcagcGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGaa                    >  7:226238/1‑148 (MQ=255)
                                                                                                                            gcagcGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAAc                   <  7:305242/149‑1 (MQ=255)
                                                                                                                               gcGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAgg                <  1:77247/149‑1 (MQ=255)
                                                                                                                                  gTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGa             <  7:75509/149‑1 (MQ=255)
                                                                                                                                     ccGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACGGGTGAAt           <  7:15092/148‑1 (MQ=255)
                                                                                                                                     ccGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGAAt           <  5:297266/148‑1 (MQ=255)
                                                                                                                                       gACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGAATGGc        >  4:158184/1‑149 (MQ=255)
                                                                                                                                             gTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGAATGGCCAGccc  >  2:51726/1‑149 (MQ=255)
                                                                                                                                                   aGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGt               <  3:129121/130‑1 (MQ=255)
                                                                                                                                                   aGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGt               >  4:129121/1‑130 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCAAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAGGAACAGGTGAATGGCCAGCCC  >  NC_002947/4546627‑4546916

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