Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,157,230 C→T 8.9% H240Y (CAC→TAC)  yihS → aldose‑ketose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,157,2300CT 8.9% 124.6 / 2.5 45H240Y (CAC→TAC) yihSaldose‑ketose isomerase
Reads supporting (aligned to +/- strand):  ref base C (21/20);  new base T (2/2);  total (23/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 = 6.48e-01

AGACGTCACCGGCAACGGCCTGTGGCTGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCGCGCATTGTTC  >  NC_002947/1157085‑1157376
                                                                                                                                                 |                                                                                                                                                  
aGACGTCACCGGCAACGGCCTGTGGCTGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGTCTATCGGGTTATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACAcc                                                                                                                                                 <  1:335539/149‑1 (MQ=255)
    gTCACCGGCAACGGCCTGTGGCTGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCccgc                                                                                                                                             <  8:14916/149‑1 (MQ=255)
              aCGGCCTGTGGCTGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTc                                                                                                                                   <  1:160852/149‑1 (MQ=255)
                         cTGCAACGTGCCTTTCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCTCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATgg                                                                                                                        <  7:289044/149‑1 (MQ=255)
                          tGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATgg                                                                                                                        >  2:186858/1‑148 (MQ=255)
                              aCGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACt                                                                                                                    >  8:133551/1‑148 (MQ=255)
                               cGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTac                                                                                                                  <  6:68444/149‑1 (MQ=255)
                                gTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTaca                                                                                                                 >  1:27046/1‑149 (MQ=255)
                                gTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTaca                                                                                                                 <  3:331360/149‑1 (MQ=255)
                                     ttGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAATACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCgg                                                                                                            >  3:307944/1‑149 (MQ=255)
                                            tCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCt                                                                                                                             <  6:306604/125‑1 (MQ=255)
                                            tCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCt                                                                                                                             >  5:306604/1‑125 (MQ=255)
                                                 gAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGa                                                                                                <  6:350370/149‑1 (MQ=255)
                                                 gAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGa                                                                                                <  3:145164/149‑1 (MQ=255)
                                                      catcatCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTggg                                                                                           >  2:352283/1‑149 (MQ=255)
                                                          atcCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAATACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGgcgc                                                                                        <  4:307944/148‑1 (MQ=255)
                                                           tcCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAATACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCgg                                                                                      >  5:145229/1‑149 (MQ=255)
                                                            ccATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGt                                                                                     >  7:81596/1‑149 (MQ=255)
                                                             cATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGACCACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGtt                                                                                    <  3:225638/149‑1 (MQ=255)
                                                                    aTGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTg                                                                             <  7:133551/149‑1 (MQ=255)
                                                                    aTGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTg                                                                             >  5:360932/1‑149 (MQ=255)
                                                                         aCCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCt                                                                        >  5:90668/1‑149 (MQ=255)
                                                                                  ggCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGaagcaag                                                               <  1:186858/149‑1 (MQ=255)
                                                                                       tCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCtt                                                           >  6:261166/1‑148 (MQ=255)
                                                                                          ggTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAAc                                                       >  5:265535/1‑149 (MQ=255)
                                                                                              aTCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGggc                                                   <  1:352283/149‑1 (MQ=255)
                                                                                                 gAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGggccg                                                 <  6:360932/148‑1 (MQ=255)
                                                                                                  atcATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGggccggc                                               <  1:108556/147‑1 (MQ=255)
                                                                                                  aGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTg                                                                                                  >  2:144526/1‑98 (MQ=255)
                                                                                                  aGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTg                                                                                                  <  1:144526/98‑1 (MQ=255)
                                                                                                         cGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACg                                        >  8:203928/1‑149 (MQ=255)
                                                                                                               ccATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCg                                   >  6:223613/1‑148 (MQ=255)
                                                                                                                 aTTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTACGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCaa                                >  4:52362/1‑149 (MQ=255)
                                                                                                                     gCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATgg                             <  3:52362/148‑1 (MQ=255)
                                                                                                                     gCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGc                            >  8:145506/1‑149 (MQ=255)
                                                                                                                     gCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGc                            <  5:223613/149‑1 (MQ=255)
                                                                                                                     gCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGc                            >  2:251324/1‑149 (MQ=255)
                                                                                                                      cAGCCGCTGCCGGATTACAACCTCGAATACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCt                           <  6:145229/149‑1 (MQ=255)
                                                                                                                                 ggATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTg                >  8:315126/1‑149 (MQ=255)
                                                                                                                                        aaCCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCcgcg          >  4:126321/1‑148 (MQ=255)
                                                                                                                                        aaCCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCgcgc         <  5:261166/149‑1 (MQ=255)
                                                                                                                                          ccTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCACAATGGCTGCCGGACAGTGCCCGCGCAt       >  6:43403/1‑149 (MQ=255)
                                                                                                                                          ccTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAACCCTTGAACGGGCCGGCCTGTAAGCCCCGAAATGGCTGCCGGACAGTCCCCGCGCAt       >  4:195928/1‑149 (MQ=255)
                                                                                                                                              gAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCGCGCATTGtt   <  1:251324/149‑1 (MQ=255)
                                                                                                                                               aaCACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCGCGCATTGTTc  >  8:336977/1‑149 (MQ=255)
                                                                                                                                                acacacCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCGCGCATTGTTc  >  2:224295/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AGACGTCACCGGCAACGGCCTGTGGCTGCAACGTGCCTTGCGTATCGCCGAACGCATCATCCATACCCATGCCACCGCCAACGGCTATCGGGTGATCGAGCATTTCGATGTCCATTGGCAGCCGCTGCCGGATTACAACCTCGAACACCCCGCCGACCCCTTCCGGCCTTATGGCACTACACCCGGGCATGCCCTGGAGTGGGCGCGGTTGCTGTTGCACCTGGAAGCAAGCCTTGAACGGGCCGGCCTGTACGCCCCGCAATGGCTGCCGGACAGTGCCCGCGCATTGTTC  >  NC_002947/1157085‑1157376

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