Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,819,760 G→T 10.8% A152E (GCG→GAG)  chbF ← cryptic phospho‑beta‑glucosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,819,7600GT10.8% 141.7 / 3.3 46A152E (GCG→GAG) chbFcryptic phospho‑beta‑glucosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base T (3/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ATATTACACACGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTACGATGTCAAAA  >  W3110S.gb/1819694‑1819828
                                                                  |                                                                    
atatTACACACGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCccgccg                                                                   <  1:1708854/70‑1 (MQ=255)
   ttACACACGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCTCCGGGt                                                                >  1:1580718/1‑70 (MQ=255)
       acacGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAg                                                             <  1:3074812/69‑1 (MQ=255)
        cacGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGt                                                                <  1:3560858/65‑1 (MQ=255)
         acGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTg                                                           <  1:3081208/69‑1 (MQ=255)
          cGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAg                                                        <  1:404898/71‑1 (MQ=255)
           gCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGaa                                                         >  1:2184185/1‑69 (MQ=255)
             cGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGtt                                                      <  1:1270712/70‑1 (MQ=255)
             cGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTa                                                     <  1:1527132/71‑1 (MQ=255)
               aTAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGtt                                                      <  1:853448/68‑1 (MQ=255)
                 aaaGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCa                                                 >  1:542433/1‑71 (MQ=255)
                 aaaGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCa                                                 >  1:4179110/1‑71 (MQ=255)
                 aaaGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCa                                                 >  1:430401/1‑71 (MQ=255)
                 aaaGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCa                                                 >  1:2046443/1‑71 (MQ=255)
                 aaaGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCa                                                 >  1:1930856/1‑71 (MQ=255)
                   aGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGaa                                                         >  1:2560412/1‑61 (MQ=255)
                   aGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCAcc                                               <  1:1840144/71‑1 (MQ=255)
                       tttAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTaa                                                    <  1:3851877/62‑1 (MQ=255)
                       tttAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATg                                           <  1:2058416/71‑1 (MQ=255)
                         tAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTaa                                                    >  1:3794481/1‑60 (MQ=255)
                           aaTCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCAt                                        <  1:3339381/70‑1 (MQ=255)
                               cGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGa                                   >  1:2610043/1‑71 (MQ=255)
                                ggTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTGGt                                                            <  1:1154692/45‑1 (MQ=255)
                                    tGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGt                                                                <  1:3361706/37‑1 (MQ=255)
                                     gACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGa                                                          <  1:143810/42‑1 (MQ=255)
                                      aCGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGtt                            <  1:3179348/71‑1 (MQ=255)
                                         aTAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTc                                      >  1:1594313/1‑58 (MQ=255)
                                          tAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCg                                     <  1:106714/58‑1 (MQ=255)
                                          tAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGttcttc                        <  1:15972/71‑1 (MQ=255)
                                           aaaCGGCTTCAGTGACCATTCCCGCCGGGTTAGTGa                                                          >  1:930384/1‑36 (MQ=255)
                                           aaaCGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATg                                           >  1:588991/1‑51 (MQ=255)
                                            aaCGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACa                                 >  1:2106124/1‑60 (MQ=255)
                                             aCGGCTTCAGTGACCATTCCCTCCGGGTTAGTGAAGTTAATCACCCATGCATTCgg                                    <  1:3326927/56‑1 (MQ=255)
                                             aCGGCTTCAGTGACCATTCCCTCCGGGTTAGTGAAGTTAATCACCCATGCATTCgg                                    <  1:821137/56‑1 (MQ=255)
                                               ggCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACAt                   <  1:32236/71‑1 (MQ=255)
                                               ggCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACAt                   <  1:3608436/71‑1 (MQ=255)
                                                    cAGTGACCATTCCCTCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCttt               >  1:3285882/1‑70 (MQ=255)
                                                    cAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCttt               >  1:1239854/1‑70 (MQ=255)
                                                    cAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCttt               >  1:1557950/1‑70 (MQ=255)
                                                    cAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCttt               >  1:1660835/1‑70 (MQ=255)
                                                    cAGTGAACATTCCCTCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCtt                >  1:856016/1‑69 (MQ=255)
                                                      gTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTAc            <  1:2474106/71‑1 (MQ=255)
                                                      gTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTAc            <  1:253821/71‑1 (MQ=255)
                                                       tGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGttc                           <  1:2483527/55‑1 (MQ=255)
                                                           cATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTACGATg        <  1:2499469/70‑1 (MQ=255)
                                                               cccGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCtttt              >  1:895307/1‑60 (MQ=255)
                                                               cccGCCGGGTTAGTGAAGTTAATCACCCATGCATTAGGACAAAGTTCTTCGACATCTTTTa             >  1:850885/1‑61 (MQ=255)
                                                               cacgccgGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGttct                          >  1:2610141/3‑48 (MQ=255)
                                                                ccgccgGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTACGATGTCaaaa  <  1:3022205/71‑1 (MQ=255)
                                                                  |                                                                    
ATATTACACACGCCGATAAAGCGTTTAAATCCGGTATGACGATAAACGGCTTCAGTGACCATTCCCGCCGGGTTAGTGAAGTTAATCACCCATGCATTCGGACAAAGTTCTTCGACATCTTTTACGATGTCAAAA  >  W3110S.gb/1819694‑1819828

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: