Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,655,114 C→T 34.6% intergenic (‑100/+78) pepB ← / ← yfhJ aminopeptidase B/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,655,1140CT34.6% 25.5 / 20.1 26intergenic (‑100/+78)pepB/yfhJaminopeptidase B/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/8);  new base T (3/6);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTG  >  W3110S.gb/2655051‑2655175
                                                               |                                                             
gCAGCCCCAGGATTTTTCAGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGc                                                        <  1:1739042/71‑1 (MQ=255)
      ccAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGAt                                                  <  1:891202/71‑1 (MQ=255)
      ccAGGACTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                             >  1:10001/1‑60 (MQ=255)
             ttttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCa                                                             >  1:4115541/1‑53 (MQ=255)
             ttttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAg                                                            >  1:4305273/1‑54 (MQ=255)
              tttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtgt                                          >  1:4164401/1‑71 (MQ=255)
               ttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTAtgt                                          <  1:4113747/70‑1 (MQ=255)
                tCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCAt                                                      >  1:2314761/1‑57 (MQ=255)
                 ccGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTATGTGAc                                       >  1:3986958/1‑71 (MQ=255)
                   gCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGa                                                          <  1:1168820/50‑1 (MQ=255)
                   gCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGa                                                          <  1:3568518/50‑1 (MQ=255)
                     tAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGAt                                                  <  1:1882740/56‑1 (MQ=255)
                      aaaTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgggg                                  >  1:2326471/1‑71 (MQ=255)
                       aaTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTATGTGACTgggg                                  >  1:1863704/1‑70 (MQ=255)
                         tGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGt                                 >  1:295079/1‑69 (MQ=255)
                              cgcgTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAg                                                            >  1:1460776/1‑37 (MQ=255)
                                 gTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgcgc                       <  1:2026915/71‑1 (MQ=255)
                                 gTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgcgc                       <  1:730254/71‑1 (MQ=255)
                                   tGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTATGTGACTGGGGTGAACGGGCGCAg                     >  1:4437513/1‑71 (MQ=255)
                                     cAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGcc                   <  1:192596/71‑1 (MQ=255)
                                      aGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGcc                   <  1:187328/70‑1 (MQ=255)
                                        gAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGAt                                                  <  1:3481228/37‑1 (MQ=255)
                                        gAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtgt                                          <  1:4147158/45‑1 (MQ=255)
                                         aaGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgggg                                  >  1:3615400/1‑52 (MQ=255)
                                           ggcggcAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGc                    <  1:2774157/64‑1 (MQ=255)
                                                      agtgaAGCCTCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTg  <  1:526524/71‑1 (MQ=255)
                                                               |                                                             
GCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTG  >  W3110S.gb/2655051‑2655175

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: