Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,475,513 A→G 37.5% pseudogene (8408/8622 nt)
pseudogene (3110/3324 nt)
ydbA →
ydbA →
ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein
ECK1398:JW1402:b1405; predicted outer membrane protein, C‑ter fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,475,5130AG37.5% 8.7 / 12.6 16pseudogene (8408/8622 nt)
pseudogene (3110/3324 nt)
ydbA
ydbA
ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein
ECK1398:JW1402:b1405; predicted outer membrane protein, C‑ter fragment
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base G (4/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGT  >  W3110S.gb/1475454‑1475580
                                                           |                                                                   
cGCAGGTAAAGACGGCTGGAGAGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGc                                                           >  1:954369/1‑70 (MQ=255)
  cAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTTGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAg                                                        <  1:2706953/71‑1 (MQ=255)
     gTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTAcgcg                                                                 <  1:1497473/59‑1 (MQ=255)
      tAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTa                                                    >  1:1058138/1‑71 (MQ=255)
      tAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTa                                                    <  1:1462756/71‑1 (MQ=255)
       aaaGCCGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTGCGCGAAGAGCCAGCGTAc                                                   >  1:1593636/1‑71 (MQ=255)
       aaaGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTGCGCGAAGAGCCAGCGTAc                                                   >  1:1809955/1‑71 (MQ=255)
       aaaGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTGCGCGAAGAGCCAGCGTAc                                                   >  1:1919134/1‑71 (MQ=255)
       aaaGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTAc                                                   >  1:2243020/1‑71 (MQ=255)
       aaaGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTAc                                                   >  1:3110025/1‑71 (MQ=255)
          gACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTGCGAGAAGAGCCAGCGTACGGc                                                <  1:1125900/71‑1 (MQ=255)
                       ctagcgctaCGCTGGAAGGCGGGCCGAACCTGAGCTGCGCGAAGAGCCAGCGTACGGCAAGCCTGgcagg                                    <  1:296191/70‑1 (MQ=255)
                                      aaGGCGGGCCGAACCTGAGCTGCGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCAc                    >  1:2253505/1‑71 (MQ=255)
                                               cGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTc           <  1:310363/71‑1 (MQ=255)
                                                    cTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTa                <  1:1243465/61‑1 (MQ=255)
                                                          tACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCTATGACGGt  >  1:2379284/1‑69 (MQ=255)
                                                           |                                                                   
CGCAGGTAAAGACGGCTGGAGCGCTAGCGCTACGCTGGAAGGCGGGCCGAACCTGAGCTACGCGAAGAGCCAGCGTACGGCAAGCCTGGCAGGCGCAGGCAGTCAGCACTTTAACGTCGATGACGGT  >  W3110S.gb/1475454‑1475580

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: