Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,641,127 C→T 8.1% A226T (GCT→ACT)  gabP‑IV ← gamma‑aminobutyrate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,641,1270CT 8.1% 138.6 / 2.5 49A226T (GCT→ACT) gabP‑IVgamma‑aminobutyrate permease
Reads supporting (aligned to +/- strand):  ref base C (25/20);  new base T (2/2);  total (27/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 = 5.54e-01

GACCCTACCGCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAGCCGAA  >  NC_002947/4640982‑4641274
                                                                                                                                                 |                                                                                                                                                   
gACCCTACCGCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCcgc                                                                                                                                                  >  4:13814/1‑149 (MQ=255)
     tACCGCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATgg                                                                                                                                             >  4:106301/1‑149 (MQ=255)
     tACCGCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATgg                                                                                                                                             >  1:297171/1‑149 (MQ=255)
         gCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGAc                                                                                                                                         <  4:149136/149‑1 (MQ=255)
               aGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTc                                                                                                                                   <  6:195584/149‑1 (MQ=255)
                                 gCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGagaa                                                                                                                 <  7:24840/149‑1 (MQ=255)
                                   caccaGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACgag                                                                                                                   >  1:43414/1‑145 (MQ=255)
                                    accaGCGAAACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGAt                                                                                                              >  5:293033/1‑149 (MQ=255)
                                          gACCCCACCAGGAAGAGCGACGCCAGGTAGAAGCTGGCGATGCGGTACACCACCAGGTTGTTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGAt                                                                                                        <  2:338208/149‑1 (MQ=255)
                                          gACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCATCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGAt                                                                                                        <  2:26407/149‑1 (MQ=255)
                                            caccacCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCTGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGAt                                                                                                              <  8:290729/141‑1 (MQ=255)
                                            caccacCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGAt                                                                                                              >  7:290729/1‑141 (MQ=255)
                                             accaccAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAg                                                                                                     >  8:10227/1‑149 (MQ=255)
                                               caccaGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCa                                                                                                   >  4:145347/1‑149 (MQ=255)
                                                  caGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGTCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGc                                                                                                >  5:7859/1‑149 (MQ=255)
                                                  caGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGTCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGc                                                                                                >  5:12640/1‑149 (MQ=255)
                                                  caGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGc                                                                                                >  3:304358/1‑149 (MQ=255)
                                                            gACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCa                                                                                       <  4:304358/148‑1 (MQ=255)
                                                                 cAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTcc                                                                                 >  4:105462/1‑149 (MQ=255)
                                                                    gTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCa                                                                               >  1:280622/1‑148 (MQ=255)
                                                                     tAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCaaa                                                                             <  3:106301/149‑1 (MQ=255)
                                                                             ggCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGccccc                                                                                            <  8:134628/126‑1 (MQ=255)
                                                                             ggCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGccccc                                                                                            >  7:134628/1‑126 (MQ=255)
                                                                                gATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCTAGCACCGTCCCAACCCCATTGGGCa                                                                  >  3:279094/1‑149 (MQ=255)
                                                                                 aTGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCAt                                                                 <  3:145347/149‑1 (MQ=255)
                                                                                  tGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATg                                                                <  2:280622/149‑1 (MQ=255)
                                                                                    cGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCGGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTccc                                                                                <  8:96337/131‑1 (MQ=255)
                                                                                    cGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCGGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTccc                                                                                >  7:96337/1‑131 (MQ=255)
                                                                                      gTACACCACCAGGTTGGTGGCGCGGCGAATCTTGGCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATTGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAac                                                            <  3:13814/149‑1 (MQ=255)
                                                                                          accaccAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGagtgtgta                                                        >  1:29588/1‑142 (MQ=255)
                                                                                              ccaacctggtggCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTgg                                                                     <  5:91396/126‑1 (MQ=255)
                                                                                              ccaacctggtggCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTgg                                                                     >  6:91396/7‑132 (MQ=255)
                                                                                              ccaGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTg                                                    >  4:97097/1‑149 (MQ=255)
                                                                                                aGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTc                                                  <  3:97097/149‑1 (MQ=255)
                                                                                                                gAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCAACGAAGAACGAGAAGATGGTGATCAGCACGCCCCAAAGCAACGTCCCAAACCCATTGGGCATGAAAAAAACGCTGTCCAACAGCCGGATCAcc                                  >  5:171599/1‑149 (MQ=255)
                                                                                                                 aaTCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCg                                                       <  4:150591/127‑1 (MQ=255)
                                                                                                                 aaTCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCg                                                       >  3:150591/1‑127 (MQ=255)
                                                                                                                      tGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGa                             >  1:65811/1‑148 (MQ=255)
                                                                                                                       gTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCg                               <  3:105462/145‑1 (MQ=255)
                                                                                                                        tCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGGGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATAc                          <  2:297171/149‑1 (MQ=255)
                                                                                                                          tttCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACct                        >  4:62134/1‑149 (MQ=255)
                                                                                                                            tCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACctct                      >  1:54050/1‑149 (MQ=255)
                                                                                                                               ggTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGc                    >  2:94862/1‑148 (MQ=255)
                                                                                                                                    gCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCgg              <  3:62134/149‑1 (MQ=255)
                                                                                                                                            gTCAGTCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAg       <  6:7859/148‑1 (MQ=255)
                                                                                                                                            gTCAGTCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAg       <  6:12640/148‑1 (MQ=255)
                                                                                                                                            gTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAg       >  2:10843/1‑148 (MQ=255)
                                                                                                                                            gTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAg       <  1:10843/148‑1 (MQ=255)
                                                                                                                                                gCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAGCCGaa  >  6:337202/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACCCTACCGCTTTGAGCCCGGGGTCGTTCCAGGCCACCAGCGACACCACCAGGAAGATCGACGCCAGGTAGAAGATGGCGATGCGGTACACCACCAGGTTGGTGGCGCGGCGAATCTTGTCTTTCGGGTTGGCGGTTTCGTCAGCCGCGATGGTGACAATTTCGGCACCGAAGAACGAGAAGATGGTGATCAGCACGCCCCCAAGCACCGTCCCAAACCCATTGGGCATGAAACCACCGCTGTCCCACAGCCGGCTCACCCCGGATACCTCTGCCAGCGGCCAGAAGCCGAA  >  NC_002947/4640982‑4641274

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