Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,517,009 G→A 9.5% G124S (GGC→AGC)  PP_2210 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,517,0090GA 9.5% 114.3 / 2.5 42G124S (GGC→AGC) PP_2210LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base A (2/2);  total (20/22)
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.00e+00

AAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTGG  >  NC_002947/2516869‑2517155
                                                                                                                                            |                                                                                                                                                  
aagCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCt                                                                                                                                            >  6:95519/1‑149 (MQ=255)
  gcgcAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCt                                                                                                                                            >  4:331706/1‑147 (MQ=255)
   cgcAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAg                                                                                                                                         <  7:8847/149‑1 (MQ=255)
     cAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGga                                                                                                                                       <  8:297187/149‑1 (MQ=255)
     cAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCACCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGga                                                                                                                                       <  8:54733/149‑1 (MQ=255)
          cATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCt                                                                                                                                  >  2:46628/1‑149 (MQ=255)
            tGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTg                                                                                                                                <  3:2090/149‑1 (MQ=255)
              cTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgc                                                                                                                              <  1:46628/149‑1 (MQ=255)
                               gTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGc                                                                                                             <  2:102991/149‑1 (MQ=255)
                                      aCCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGa                                                                                                      >  4:23736/1‑149 (MQ=255)
                                       ccGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGAc                                                                                                     <  5:95519/149‑1 (MQ=255)
                                        cGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcac                                                                                                                                     <  7:197028/116‑4 (MQ=255)
                                        cGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcac                                                                                                                                     >  8:197028/1‑113 (MQ=255)
                                           gcgcCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTgg                                                                                                 <  3:23736/149‑1 (MQ=255)
                                               cTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcgc                                                                                                                                     <  4:308720/109‑4 (MQ=255)
                                               cTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTCGTGTCACCGTCAACTTGCGCCTGAGCATTGCCCAggcgc                                                                                                                                     >  3:308720/1‑106 (MQ=255)
                                                  gAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTg                                                                                            <  5:2525/147‑1 (MQ=255)
                                                     gTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgctgc                                                                                                                           >  2:130912/1‑113 (MQ=255)
                                                     gTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGgctgc                                                                                                                           <  1:130912/113‑1 (MQ=255)
                                                       aGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGAcc                                                                                      >  3:4584/1‑148 (MQ=255)
                                                                     cGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCcgcg                                                                        >  2:60596/1‑148 (MQ=255)
                                                                                 acccGATGCGGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCAt                                                           <  8:297615/147‑1 (MQ=255)
                                                                                   ccGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGACGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCt                                                         >  8:238068/1‑149 (MQ=255)
                                                                                    cGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGc                                                                                           <  7:211686/114‑1 (MQ=255)
                                                                                    cGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGc                                                                                           >  8:211686/1‑114 (MQ=255)
                                                                                              gctgctGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATctct                                              <  7:238068/149‑1 (MQ=255)
                                                                                                  ctgcGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCtgt                                            >  5:341562/1‑147 (MQ=255)
                                                                                                  ctgcGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGa                                          <  2:266509/149‑1 (MQ=255)
                                                                                                       tGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTcgc                                     <  6:341562/149‑1 (MQ=255)
                                                                                                           cGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCt                                 <  1:34028/149‑1 (MQ=255)
                                                                                                            ggTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCtg                                >  6:36720/1‑149 (MQ=255)
                                                                                                                 ccTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTgc                           >  2:40220/1‑149 (MQ=255)
                                                                                                                 ccTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTgc                           <  1:60596/149‑1 (MQ=255)
                                                                                                                  cTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCGGACCGAAATCGAGCCGCGCAAGGGGCTGCCTCTGCAAAATCTCTGGGAGTCGCGCCTGTGTgcg                          >  1:321799/1‑149 (MQ=255)
                                                                                                                        tCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTg                                                        >  3:229996/1‑113 (MQ=255)
                                                                                                                        tCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTg                                                        <  4:229996/113‑1 (MQ=255)
                                                                                                                          aCCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTAcc                  >  3:273119/1‑149 (MQ=255)
                                                                                                                                aaCTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTgg             <  4:273119/148‑1 (MQ=255)
                                                                                                                                       gcCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGcc     >  1:86884/1‑149 (MQ=255)
                                                                                                                                        ccTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCt    >  2:216804/1‑149 (MQ=255)
                                                                                                                                          tGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGa                                          >  2:149264/1‑109 (MQ=255)
                                                                                                                                          tGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGa                                          <  1:149264/109‑1 (MQ=255)
                                                                                                                                            ggCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTgg  <  3:331706/147‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
AAGCGCAGGCCATGCTTGAGGCCAACCGCCAGTTGCTGACCGGGCGCCTGGAGGTAGCCGCCGACGGCCCGCATCGGGTCATGCCGATGCTGGCGCTGCTGCGTGAGCGGTACCCTGGTGTCACCGTCAACTTGCGCCTGGGCAATGCTCAGGAGACCTTGGCTGCGCTGCTCGGCGAGCATGCCGACGTGGCGGTGCTGACCGAAATCGAGCCGCGCAAGGGGCTGCATCTGCAAAATCTCTGTGAGTCGCGCCTGTGTGCGCTGTTACCGGTGGGGCACGCCTGG  >  NC_002947/2516869‑2517155

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