Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,655,114 C→T 28.0% 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,1140CT28.0% 36.7 / 14.6 25intergenic (‑100/+78)pepB/yfhJaminopeptidase B/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base T (2/5);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAT  >  W3110S.gb/2655045‑2655181
                                                                     |                                                                   
cAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGcccc                                                                    <  1:451941/71‑1 (MQ=255)
                     ttCCGCTAAATGATTCGCGTTGCAGGAATGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtgt                                                <  1:995701/70‑1 (MQ=255)
                     ttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTAtgt                                                <  1:2903872/70‑1 (MQ=255)
                     ttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTAtgt                                                <  1:3443801/70‑1 (MQ=255)
                     ttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTAtgt                                                <  1:1723413/70‑1 (MQ=255)
                     ttCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTAtgt                                                <  1:977670/70‑1 (MQ=255)
                         gCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTg                                           <  1:901668/71‑1 (MQ=255)
                          cTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTgg                                          <  1:1824996/71‑1 (MQ=255)
                             aaTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGcccc                                                                    <  1:3990747/42‑1 (MQ=255)
                             aaTGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGcccc                                                                    <  1:2623219/42‑1 (MQ=255)
                               tGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTATGTGa                                              >  1:73127/1‑62 (MQ=255)
                               tGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCTCAGGAGCATAGATAACTATGTGa                                              >  1:3063101/1‑62 (MQ=255)
                               tGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAg                                                                  >  1:2887625/1‑42 (MQ=255)
                               tGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGACCCAg                                                                  >  1:2489845/1‑42 (MQ=255)
                                 aTTCTCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATa                                                           <  1:1085636/47‑1 (MQ=255)
                                            aGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgcgc                             <  1:1220512/66‑1 (MQ=255)
                                               aaGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGaa                                    >  1:2159526/1‑56 (MQ=255)
                                                aGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCa                   >  1:2384228/1‑71 (MQ=255)
                                                aGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCa                   >  1:2379522/1‑71 (MQ=255)
                                                          tgagtgaAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTgcggcg          >  1:2135826/1‑71 (MQ=255)
                                                              tgaAGCCTCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgcg                              <  1:319456/47‑1 (MQ=255)
                                                              tgaAGCCTCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAgcg                              <  1:1924045/47‑1 (MQ=255)
                                                              tgaAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCa                   <  1:3874684/58‑1 (MQ=255)
                                                              tgaAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCa                   <  1:3880176/58‑1 (MQ=255)
                                                              tgaAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCa                   <  1:2024372/58‑1 (MQ=255)
                                                               gaAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAg     >  1:2234661/1‑71 (MQ=255)
                                                                aaGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAg     >  1:389190/1‑70 (MQ=255)
                                                                  gCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAt  >  1:2981539/1‑71 (MQ=255)
                                                                     |                                                                   
CAGTTGGCAGCCCCAGGATTTTTCCGCTAAATGATTCGCGTTGCAGGAAGGCGGCAAGTGAGTGAAGCCCCAGGAGCATAGATAACTATGTGACTGGGGTGAACGAGCGCAGCCAACGCATCTGCGGCGTGAAGCAT  >  W3110S.gb/2655045‑2655181

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: