Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,286,501 C→A 30.3% A261S (GCG→TCG)  yejK ← nucleotide associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,286,5010CA30.3% 119.0 / 33.1 56A261S (GCG→TCG) yejKnucleotide associated protein
Reads supporting (aligned to +/- strand):  ref base C (15/24);  new base A (7/10);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

AACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCG  >  W3110S.gb/2286437‑2286571
                                                                |                                                                      
aaCTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTc                                                                  >  1:3812180/1‑71 (MQ=255)
aaCTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACt                                                                   >  1:3330055/1‑70 (MQ=255)
         aGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCt                                                         <  1:1356989/71‑1 (MQ=255)
           tgtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTg                                                              <  1:3590264/64‑1 (MQ=255)
           tgtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTg                                                              <  1:3433589/64‑1 (MQ=255)
           tgtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCt                                                         <  1:3278891/69‑1 (MQ=255)
            gtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCaaa                                                                     >  1:3656612/1‑56 (MQ=255)
            gtgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCaaa                                                                     >  1:3030780/1‑56 (MQ=255)
             tgCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACt                                                     >  1:4361834/1‑71 (MQ=255)
              gCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACtt                                                    >  1:2970315/1‑71 (MQ=255)
               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTc                                                   >  1:965214/1‑71 (MQ=255)
               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTg                                                              <  1:4424543/60‑1 (MQ=255)
               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTc                                                   <  1:3258304/71‑1 (MQ=255)
               cTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTc                                                   <  1:297818/71‑1 (MQ=255)
                 aCGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGaa                                                            >  1:2143440/1‑60 (MQ=255)
                  cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:1489437/71‑1 (MQ=255)
                  cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:2995670/71‑1 (MQ=255)
                  cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:2335631/71‑1 (MQ=255)
                  cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:748243/71‑1 (MQ=255)
                  cGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:2779678/71‑1 (MQ=255)
                   ggTCTGCCGGGAGGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                >  1:4624324/1‑70 (MQ=255)
                         ccGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCt                                         >  1:2589619/1‑71 (MQ=255)
                         ccGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCt                                         >  1:849511/1‑71 (MQ=255)
                              aaGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCa                                     >  1:3960335/1‑70 (MQ=255)
                               aGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGc                                   <  1:3485996/71‑1 (MQ=255)
                                gCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGct                                  <  1:14126/71‑1 (MQ=255)
                                 ctctctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:880056/56‑1 (MQ=255)
                                  tctctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGctct                                >  1:555330/1‑71 (MQ=255)
                                   ctctTCCAGTTCGTAGCCGTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCtt                               >  1:3300818/1‑71 (MQ=255)
                                    tctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:1191782/53‑1 (MQ=255)
                                    tctTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCt                                                <  1:359489/53‑1 (MQ=255)
                                    tctTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCg                                                                <  1:4290031/37‑1 (MQ=255)
                                    tctTCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCg                                                                <  1:1173658/37‑1 (MQ=255)
                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTa                                              >  1:3084791/1‑53 (MQ=255)
                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCtt                               <  1:394378/68‑1 (MQ=255)
                                      ttCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCtt                               <  1:2916065/68‑1 (MQ=255)
                                      ttCCAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTg                                                              >  1:1736825/1‑37 (MQ=255)
                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAg                                                           >  1:2284608/1‑38 (MQ=255)
                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             <  1:1100217/68‑1 (MQ=255)
                                        ccAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             <  1:4708200/68‑1 (MQ=255)
                                        ccAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             <  1:1065850/68‑1 (MQ=255)
                                        ccAGTTCGTAGCCTTTTTCAGCCGAAAACTCCGTGAACCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             <  1:1915761/68‑1 (MQ=255)
                                        acAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAg                                                           >  1:1901619/2‑38 (MQ=255)
                                         cAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTTAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             <  1:2295713/67‑1 (MQ=255)
                                            ttCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTAc                                             <  1:2699728/48‑1 (MQ=255)
                                              cGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGTCTGCCAGCTCTTTCGACAGCGa                     <  1:4525430/70‑1 (MQ=255)
                                                tAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt                  >  1:1046678/1‑71 (MQ=255)
                                                tAGCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt                  <  1:1766537/71‑1 (MQ=255)
                                                  gCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt                  >  1:3939804/1‑69 (MQ=255)
                                                  gCCTTTTTCAGCCGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGAttt                  >  1:4721423/1‑69 (MQ=255)
                                                   ccTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCg                                                  >  1:494773/1‑36 (MQ=255)
                                                     tttttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTc                 <  1:3037260/67‑1 (MQ=255)
                                                        ttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCg           <  1:2159487/70‑1 (MQ=255)
                                                        ttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCg           <  1:4563093/70‑1 (MQ=255)
                                                        ttCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCg           <  1:129831/70‑1 (MQ=255)
                                                             ccGAAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATctct     <  1:3803286/71‑1 (MQ=255)
                                                                cAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTc                             >  1:3968135/1‑44 (MQ=255)
                                                                cAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCg  >  1:1458620/1‑71 (MQ=255)
                                                                |                                                                      
AACTGACGCAGTGTGCTACGGTCTGCCGGGAAGCTCTCTTCCAGTTCGTAGCCTTTTTCAGCCGCAAACTCCGTGAAGCTCACTTCGCTTACGCCTGCCAGCTCTTTCGACAGCGATTTCAGTTCGATCTCTTCG  >  W3110S.gb/2286437‑2286571

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