Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,532,085 G→C 9.6% P79R (CCC→CGC) ‡ braE ← high‑affinity branched‑chain amino acid ABC transporter permease BraE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,532,0850GC 9.6% 123.6 / 3.1 41P79R (CCC→CGC) ‡braEhigh‑affinity branched‑chain amino acid ABC transporter permease BraE
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base C (2/2);  total (23/18)
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.95e-01

GCAGGTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGG  >  NC_002947/5531939‑5532231
                                                                                                                                                  |                                                                                                                                                  
gcagGTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCggggg                                                                                                                                                 <  1:494528/150‑1 (MQ=255)
    gTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCggggg                                                                                                                                                 <  1:777588/146‑1 (MQ=255)
    gTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCggggg                                                                                                                                                 >  2:777588/1‑146 (MQ=255)
    gTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGt                                                                                                                                             >  2:1078615/1‑150 (MQ=255)
     tACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTc                                                                                                                                            >  2:124651/1‑150 (MQ=255)
        ttgttgGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCggg                                                                                                                                                   >  2:311740/1‑140 (MQ=255)
        ttgttgGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCggg                                                                                                                                                   <  1:311740/140‑1 (MQ=255)
                                   gCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCa                                                                                                                         >  1:373482/1‑139 (MQ=255)
                                   gCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCa                                                                                                                         <  2:373482/139‑1 (MQ=255)
                                     aaTCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCCCGGGCAgcgca                                                                                                                                          >  1:481821/1‑116 (MQ=255)
                                     aaTCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCCCGGGCAgcgca                                                                                                                                          <  2:481821/120‑5 (MQ=255)
                                           gATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCa                                                                                                                                        >  1:99685/1‑116 (MQ=255)
                                           gATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCa                                                                                                                                        <  2:99685/116‑1 (MQ=255)
                                                           atgTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAgg                                                                                                                                >  2:214727/1‑108 (MQ=255)
                                                           atgTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAgg                                                                                                                                <  1:214727/108‑1 (MQ=255)
                                                            ggTAACGCAGGCGCGACTTGTAGTCCGGTGCGGCCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAgg                                                                                     <  1:112396/149‑1 (MQ=255)
                                                                      gcgcGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCCCGGGCAgcgca                                                                                                                                          <  1:287445/87‑5 (MQ=255)
                                                                      gcgcGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCCCGGGCAgcgca                                                                                                                                          >  2:287445/1‑83 (MQ=255)
                                                                           aCTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGt                                                                                                                                             <  2:1057628/79‑1 (MQ=255)
                                                                           aCTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGt                                                                                                                                             >  1:1057628/1‑79 (MQ=255)
                                                                            cTTGTAGTCCGGTGCGGTCACATGCACGCCCGATCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTc                                                                     >  2:531214/1‑150 (MQ=255)
                                                                            cTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCTACCCGCCAGGCCACGCGGCGCGGCTc                                                                     <  2:138901/150‑1 (MQ=255)
                                                                             ttGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCg                                                                    >  1:586773/1‑150 (MQ=255)
                                                                               gTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTg                                                                                                                                          <  2:485121/78‑1 (MQ=255)
                                                                               gTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTg                                                                                                                                          >  1:485121/1‑78 (MQ=255)
                                                                                       gTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGaa                                                                                                                     <  1:276629/91‑1 (MQ=255)
                                                                                       gTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGaa                                                                                                                     >  2:276629/1‑91 (MQ=255)
                                                                                           ggTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGa                                                                                                                            >  2:844836/1‑80 (MQ=255)
                                                                                           ggTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGa                                                                                                                            <  1:844836/80‑1 (MQ=255)
                                                                                                 aTGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTc                                                >  1:1029546/1‑150 (MQ=255)
                                                                                                         ccGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCCGGCACATGCTGAGCAGGAGACGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGcacca                                        >  1:451113/1‑150 (MQ=255)
                                                                                                               cACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCaa                                  >  2:518538/1‑150 (MQ=255)
                                                                                                                aCCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTc                                                                     >  2:309087/1‑114 (MQ=255)
                                                                                                                aCCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTc                                                                     <  1:309087/114‑1 (MQ=255)
                                                                                                                aCCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAAc                                 >  2:523932/1‑150 (MQ=255)
                                                                                                                    cTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGAcc                             <  2:1029546/150‑1 (MQ=255)
                                                                                                                    cTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGAcc                             <  2:586773/150‑1 (MQ=255)
                                                                                                                        cGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGAc                         <  1:725484/150‑1 (MQ=255)
                                                                                                                         gAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGAcc                        >  2:261477/1‑150 (MQ=255)
                                                                                                                            gCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCaaa                     >  2:706233/1‑150 (MQ=255)
                                                                                                                              ccAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGcaccacacc                                    >  2:879428/1‑133 (MQ=255)
                                                                                                                              ccAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGcaccacacc                                    <  1:879428/133‑1 (MQ=255)
                                                                                                                              ccAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAaca                   <  1:523932/150‑1 (MQ=255)
                                                                                                                                        gcgcTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGc         <  1:518538/150‑1 (MQ=255)
                                                                                                                                              cccGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCcagcag   >  2:232027/1‑150 (MQ=255)
                                                                                                                                               ccGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAgg  >  1:257019/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCAGGTACTTGTTGGCGAAGATCGGGAAAACGATGGCAATCACGATCAGCGCGGGGATGATGTAACGCAGGCGCGACTTGTAGTCCGGTGCGGTCACATGCACGCCCGAGCCACCGCTGTCGAAGCCCAGGAGCATGCGCTGCCCGGGGGCAGTCTGCAGGAACAGGCTGAGCAGGAAGCGCCCCACCATCACACCGCCGACCAGCCAGGCCACGCGGCGCGGCTCGGCATTGAAGGTGTAGCCGTCGAGCACCACACCAACGACCGGACCAAACACGATAAGTGCCAGCAGG  >  NC_002947/5531939‑5532231

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