Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,546,755 C→T 8.1% P197P (CCG→CCA)  pydP ← NCS1 family transporter PydP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,546,7550CT 8.1% 136.8 / 2.2 49P197P (CCG→CCA) pydPNCS1 family transporter PydP
Reads supporting (aligned to +/- strand):  ref base C (19/26);  new base T (2/2);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

GCCCACCATGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAG  >  NC_002947/4546611‑4546895
                                                                                                                                                  |                                                                                                                                            
gCCCACTATGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTg                                                                                                                                          <  3:159557/149‑1 (MQ=255)
  ccaccaTGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGCACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGcc                                                                                                                                        >  1:83493/1‑149 (MQ=255)
   caccaTGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGcc                                                                                                                                        >  3:190729/1‑148 (MQ=255)
   caccaTGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGcc                                                                                                                                        >  3:184609/1‑148 (MQ=255)
    accaTGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGa                                                                                                                                      <  3:96644/149‑1 (MQ=255)
            cGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGttt                                                                                                                              <  1:109209/149‑1 (MQ=255)
             gTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACACCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTc                                                                                                                             >  4:96358/1‑149 (MQ=255)
             gTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACACCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTc                                                                                                                             <  3:96358/149‑1 (MQ=255)
                          aGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATc                                                                                                                <  1:392350/149‑1 (MQ=255)
                             aGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCAACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGa                                                                                                              >  4:122790/1‑148 (MQ=255)
                              gTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACt                                                                                                            >  4:280465/1‑149 (MQ=255)
                                   tGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGCTGAGGGCCCAGAACAGCAGACAGACGCCCACTGCCACCAGCAGCGGTGCCGCCAGCGTTTCCAGCCACTGGATCGACTCgg                                                                                                         >  1:319992/1‑149 (MQ=255)
                                        aCGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACTGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTTATCGACTCGGCGCcgcg                                                                                                  <  4:225576/149‑1 (MQ=255)
                                                cccTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGcacca                                                                                          >  6:290087/1‑149 (MQ=255)
                                                     ccGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCTGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCAccaccc                                                                                       <  1:308430/147‑1 (MQ=255)
                                                       ggACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCAccaccc                                                                                       <  8:220826/145‑1 (MQ=255)
                                                       ggACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCAccaccc                                                                                       >  7:220826/1‑145 (MQ=255)
                                                        gACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACTTGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAgg                                                                                  <  4:184609/149‑1 (MQ=255)
                                                        gACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAgg                                                                                  <  4:190729/149‑1 (MQ=255)
                                                           gCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGtt                                                                                >  6:19480/1‑148 (MQ=255)
                                                            cTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCa                                                                              <  8:33487/149‑1 (MQ=255)
                                                             ttGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCgcag                                                                                               <  5:369482/131‑1 (MQ=255)
                                                             ttGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCgcag                                                                                               >  6:369482/1‑131 (MQ=255)
                                                                ggcggcTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGc                                                                           <  2:319992/148‑1 (MQ=255)
                                                                        cgcCAACATTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCagaag                                                                   <  6:27053/148‑1 (MQ=255)
                                                                          ccAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATc                                                                <  3:280465/149‑1 (MQ=255)
                                                                            aaCAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCAt                                                              <  5:290087/149‑1 (MQ=255)
                                                                                gTTCGGTAATCGACATGTGAGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCAACTGCCACCAGCAGCGGTGCAGACAGCGTTTCCAGCCACTTGATCGACTCGGCGACGCGCAGCACCAACCACAGGTTCAGGCACAAGAAGATCATGAAg                                                          >  6:297128/1‑149 (MQ=255)
                                                                                           gACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCt                                                <  7:121105/148‑1 (MQ=255)
                                                                                           gACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCACACTGCCACCAGAAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTc                                               <  3:290904/149‑1 (MQ=255)
                                                                                            aCATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCg                                              >  3:7808/1‑149 (MQ=255)
                                                                                            aCATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCg                                              >  3:31098/1‑149 (MQ=255)
                                                                                                   cGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGc                                                         <  7:344881/131‑1 (MQ=255)
                                                                                                   cGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGc                                                         >  8:344881/1‑131 (MQ=255)
                                                                                                    ggC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGCCAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTccc                                      >  3:68081/1‑149 (MQ=255)
                                                                                                    ggC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGCGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTccc                                      <  6:7463/149‑1 (MQ=255)
                                                                                                      c‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGTGGTGCCGACAGCGTTTCTAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCAGCAGAAGATCATGAAGCCAATCACCTCGc                                             <  1:222313/140‑1 (MQ=255)
                                                                                                      c‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGTGGTGCCGACAGCGTTTCTAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGc                                             >  2:222313/1‑140 (MQ=255)
                                                                                                      c‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGTGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGc                                             >  6:167827/1‑140 (MQ=255)
                                                                                                      c‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGTGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGc                                             <  5:167827/140‑1 (MQ=255)
                                                                                                             cccATAACAGCAGAACGACGCCCACTGCCACCAGAAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCATAAGATCATGAAGCCAATCACCTAGCCTGTCCCGCCCatc                               <  5:297128/149‑3 (MQ=255)
                                                                                                             cccAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCTCAGGTTCAGGCGCCAGACGCTCATGAAGCGAATCACCTCGCCTGTCCCGCCCGgc                               >  6:325187/1‑149 (MQ=255)
                                                                                                             cccAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAgc                               >  7:140810/1‑149 (MQ=255)
                                                                                                                     agcagACCGACTCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCa                       <  1:369517/149‑1 (MQ=255)
                                                                                                                     agcagACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCa                       <  5:19480/149‑1 (MQ=255)
                                                                                                                        agACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGcc                    <  4:31098/149‑1 (MQ=255)
                                                                                                                              aCGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGaa              >  5:368664/1‑149 (MQ=255)
                                                                                                                                        ccaccaGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAAc      <  3:363501/147‑1 (MQ=255)
                                                                                                                                          accaGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAg  <  4:7808/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCCCACCATGGCCGTAAGCCCGGCGCAGAAGTAACTGACCACGCTCGCCCCTTCCGGACGCTTGGGCGGCTGCGCCAACAGTTCGGTCATCGACATGTGCGGC‑‑AAGGCCCAGAACAGCAGACCGACGCCCACTGCCACCAGCAGCGGTGCCGACAGCGTTTCCAGCCACTTGATCGACTCGGCGCCGCGCAGCACCACCCACAGGTTCAGGCACCAGAAGATCATGAAGCCAATCACCTCGCCTGTCCCGCCCAGCGCCTTCCAGCCGTCGAACACCGAACCCAG  >  NC_002947/4546611‑4546895

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: ‑