Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 552,984 C→G 15.1% A60A (GCG→GCC lpxH ← UDP‑2,3‑diacylglucosamine pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb552,9840CG15.1% 73.5 / 5.9 33A60A (GCG→GCClpxHUDP‑2,3‑diacylglucosamine pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (16/12);  new base G (3/2);  total (19/14)
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.88e-01

CAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATA  >  W3110S.gb/552915‑553053
                                                                     |                                                                     
cAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCTCTTTGATcgccgcc                                                                       <  1:3086948/70‑1 (MQ=255)
cAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATcgccgcc                                                                       <  1:1894532/70‑1 (MQ=255)
      aTCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCt                                                                <  1:433104/71‑1 (MQ=255)
      aTCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCt                                                                <  1:3260661/71‑1 (MQ=255)
               gCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATg                                                         >  1:2577939/1‑69 (MQ=255)
                     aatgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTg                                                              >  1:3194369/1‑58 (MQ=255)
                      atgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGt                                                 <  1:3693105/70‑1 (MQ=255)
                       tgaaATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATgg                                                        <  1:341538/62‑1 (MQ=255)
                         aaaTAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTgg                                                   >  1:3661787/1‑65 (MQ=255)
                         aaaTAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTggg                                             <  1:1605388/71‑1 (MQ=255)
                            tAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATgg                                                        <  1:302185/57‑1 (MQ=255)
                             aaCAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAtc                                          <  1:2358767/70‑1 (MQ=255)
                                aGGGAACGCCGGAATCGGACACCGCTTTGATCgccgccgcc                                                                    >  1:3596802/1‑41 (MQ=255)
                                        ccGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAAt                              >  1:2104609/1‑71 (MQ=255)
                                         cGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAtcgtc                                       >  1:2124701/1‑61 (MQ=255)
                                             aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAg                                                      >  1:1517047/1‑42 (MQ=255)
                                             aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAg                                                      >  1:2476390/1‑42 (MQ=255)
                                             aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                         >  1:2250739/1‑71 (MQ=255)
                                             aTCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATg                         >  1:743152/1‑71 (MQ=255)
                                              tCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTggg                                             <  1:3728935/50‑1 (MQ=255)
                                                ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt                                           >  1:906814/1‑50 (MQ=255)
                                                ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt                                           >  1:3370827/1‑50 (MQ=255)
                                                ggACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGAt                                           >  1:514977/1‑50 (MQ=255)
                                                  acacCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                    <  1:3441602/71‑1 (MQ=255)
                                                   tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                    >  1:2131296/2‑70 (MQ=255)
                                                   tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                    >  1:700069/2‑70 (MQ=255)
                                                   tacCGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCa                    >  1:380728/2‑70 (MQ=255)
                                                   cacCGCTTTGAGCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCa                    <  1:1702986/70‑1 (MQ=255)
                                                     ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                 <  1:1991868/71‑1 (MQ=255)
                                                     ccGCTTTGATTGCTGCGGCCATCTGGTGATGAAGTGGGTTGGGATCGTCGTCGCCGATCCAGGCTTCAAAc                 <  1:141476/71‑1 (MQ=255)
                                                     ccGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa                  >  1:1629352/1‑70 (MQ=255)
                                                      cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa                <  1:1679859/71‑1 (MQ=255)
                                                      cGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACa                <  1:235084/71‑1 (MQ=255)
                                                       gCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCaaa                  >  1:1808479/1‑68 (MQ=255)
                                                           tGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCg           >  1:16782/1‑71 (MQ=255)
                                                            gATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCAATGCTACAAACAGATCGc          >  1:3305373/1‑71 (MQ=255)
                                                                    ccgccATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAata  <  1:2444161/71‑1 (MQ=255)
                                                                     |                                                                     
CAGAAAATCACGGTTGCCATGAATGAAATAACAGGGAACGCCGGAATCGGACACCGCTTTGATCGCCGCCGCCATCTTGCGATGGAGTGGGTTGGGATCGTCGTCGCCAATCCATGCTTCAAACAGATCGCCAAGAATA  >  W3110S.gb/552915‑553053

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