Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,497,222 A→G 25.0% Y59H (TAT→CAT)  oxc ← predicted oxalyl‑CoA decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,497,2220AG25.0% 60.8 / 20.8 36Y59H (TAT→CAT) oxcpredicted oxalyl‑CoA decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (15/12);  new base G (5/4);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCTGGCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGCCTGGGCATG  >  W3110S.gb/2497153‑2497279
                                                                     |                                                         
aTCCTGGCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAg                                                          <  1:602282/71‑1 (MQ=255)
aTCCTGGCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAg                                                          <  1:2213929/71‑1 (MQ=255)
   cTGGCGCAGAAACTGTCAGGCAGATCCCTGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATGGCCt                                                       >  1:501858/1‑71 (MQ=255)
   cTGGCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATGGCCt                                                       >  1:493688/1‑71 (MQ=255)
     ggCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCCTGGCCTGc                                                     >  1:3934257/1‑71 (MQ=255)
     ggCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGc                                                     >  1:337815/1‑71 (MQ=255)
       cgcAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGcc                                                    <  1:3129191/70‑1 (MQ=255)
         cAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGcc                                                    <  1:118620/68‑1 (MQ=255)
             aaCTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCt                                             >  1:3572201/1‑71 (MQ=255)
             aaCTGTCAGGCAGATCCACGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCt                                             >  1:2322376/1‑71 (MQ=255)
                    aGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTAATGAc                                       >  1:691113/1‑70 (MQ=255)
                        aGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGaaaa                                  >  1:1326067/1‑71 (MQ=255)
                         gATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGaa                                    <  1:1459465/68‑1 (MQ=255)
                         gATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGaa                                    <  1:843001/68‑1 (MQ=255)
                          aTCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATGGCCTGCCGACTGCTCATGACGAAAAcc                                <  1:94846/71‑1 (MQ=255)
                            ccccGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGTCATAGCCTGCCGACTGCTCATGACGAAAACCaa                              >  1:1599185/1‑71 (MQ=255)
                            ccccGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCaa                              >  1:632622/1‑71 (MQ=255)
                            ccccGGTTTTTGGGTAAGAAAACCGCTCGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCaa                              >  1:3376149/1‑71 (MQ=255)
                                    tttGGGTAAGAAAACCGCTTGCCGCAGCGGCATGGCCTGCCGACTGCTCATGACGAAAACCAATATAACg                       >  1:1899978/1‑70 (MQ=255)
                                    tttGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACg                       <  1:1209434/70‑1 (MQ=255)
                                         gTAAGAAAACCGCTTGCCGCAGCGGCATGGCCTGCCGACTGCTCATGACGAAAACCAatat                           <  1:1867255/61‑1 (MQ=255)
                                         gTAAGAAAACCGCTTGCCGCAGCGGCATGGCCTGCCGACTGCTCATGACGAAAACCAATATAACGaa                     <  1:737382/67‑1 (MQ=255)
                                              aaaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATg                                         >  1:2165976/1‑42 (MQ=255)
                                              aaaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGaa                     >  1:636419/1‑62 (MQ=255)
                                              aaaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGaa                     >  1:2164264/1‑62 (MQ=255)
                                              aaaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGaa                     >  1:1890040/1‑62 (MQ=255)
                                              aaaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTcc             <  1:2340301/70‑1 (MQ=255)
                                               aaaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCa                                           >  1:693225/1‑39 (MQ=255)
                                                aaCCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGcc          <  1:2768673/71‑1 (MQ=255)
                                                    gCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGCCTggg      <  1:3612274/71‑1 (MQ=255)
                                                        gCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGCCTGGGCATg  >  1:2534259/1‑71 (MQ=255)
                                                           gCAGCGGCATGGCCTGCCGACTGCTCATGACGAAAACCAATATAACg                       <  1:296357/47‑1 (MQ=255)
                                                             aGCGGCATGGCCTGCCGACTGCTCATGACGAAAACCAATATAACGaa                     >  1:888618/1‑47 (MQ=255)
                                                              gCGGCATAGCCTGCCGACTGCTCATGACGAAAAcc                                <  1:2134802/35‑1 (MQ=255)
                                                                   aTAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGCCt         >  1:3711002/1‑53 (MQ=255)
                                                                    tAGCCTGCCGACTGCTCATGACGAAAACCAATATAACg                       <  1:3015251/38‑1 (MQ=255)
                                                                     |                                                         
ATCCTGGCGCAGAAACTGTCAGGCAGATCCCCGGTTTTTGGGTAAGAAAACCGCTTGCCGCAGCGGCATAGCCTGCCGACTGCTCATGACGAAAACCAATATAACGAATGCCTTCCGCCTGGGCATG  >  W3110S.gb/2497153‑2497279

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: