Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,007,077 C→T 34.7% N304N (AAC→AAT ygeX → 2,3‑diaminopropionate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,007,0770CT34.7% 88.1 / 72.9 75N304N (AAC→AATygeX2,3‑diaminopropionate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (28/21);  new base T (17/9);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCCACCATCATGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAAT  >  W3110S.gb/3007010‑3007143
                                                                   |                                                                  
tGGCCACCATCATGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTgc                                                                 >  1:756289/1‑71 (MQ=255)
tGGCCACCATCATGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTg                                                                  <  1:3705075/70‑1 (MQ=255)
  gCCACCATCATGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTgcgc                                                               <  1:951175/71‑1 (MQ=255)
        atcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTgcg                                                                >  1:752984/1‑64 (MQ=255)
        atcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCAccc                                                          >  1:1150757/1‑70 (MQ=255)
        atcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCAccc                                                          >  1:4372790/1‑70 (MQ=255)
        atcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCAccc                                                          >  1:4559775/1‑70 (MQ=255)
        atcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCAccc                                                          >  1:2621546/1‑70 (MQ=255)
         tcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCaa                                                        >  1:4142527/1‑71 (MQ=255)
         tcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCaa                                                        >  1:427842/1‑71 (MQ=255)
         tcatGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCaa                                                        <  1:3333118/71‑1 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAtt                                                      >  1:3168549/1‑71 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAtt                                                      >  1:4279709/1‑71 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAtt                                                      >  1:2672611/1‑71 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAtt                                                      >  1:2984530/1‑71 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAtt                                                      >  1:4494349/1‑71 (MQ=255)
           atGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAAt                                                       >  1:3375460/1‑70 (MQ=255)
              gCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCAt                                                   >  1:580984/1‑71 (MQ=255)
              gCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCa                                                    <  1:692220/70‑1 (MQ=255)
               cAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTgcgc                                                               >  1:3825202/1‑58 (MQ=255)
               cAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCAtc                                                  <  1:4324490/71‑1 (MQ=255)
                    ctggcctgTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTg                                             <  1:627204/71‑1 (MQ=255)
                       gcctgTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAAtt                                                      >  1:3004159/1‑59 (MQ=255)
                       gcctgTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAAtt                                                      >  1:977602/1‑59 (MQ=255)
                       gcctgTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGcc                                           <  1:937287/70‑1 (MQ=255)
                        cctgTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAg                                         <  1:1786883/71‑1 (MQ=255)
                          tgtgGCGAACCTAACCCGCTGGGCTGTGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAg                                         <  1:3246919/69‑1 (MQ=255)
                          tgtgGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAAtt                                                      <  1:1932300/56‑1 (MQ=255)
                            tgGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGa                                       <  1:3531606/69‑1 (MQ=255)
                            tgGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGAc                                      <  1:1979902/70‑1 (MQ=255)
                              gCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACacc                                   >  1:5118713/1‑69 (MQ=255)
                              gCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGc                                   >  1:3619311/1‑71 (MQ=255)
                              gCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGc                                   >  1:2415207/1‑71 (MQ=255)
                              gCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGc                                   >  1:3400642/1‑71 (MQ=255)
                               cGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCg                                  <  1:1254018/71‑1 (MQ=255)
                                 aaCCTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGtt                                >  1:1972103/1‑71 (MQ=255)
                                 aaCCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAg                                         >  1:5079525/1‑62 (MQ=255)
                                   ccTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGt                                 >  1:648151/1‑68 (MQ=255)
                                   ccTAACCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGt                                 >  1:395557/1‑68 (MQ=255)
                                      aaCCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCa                                                    >  1:174927/1‑46 (MQ=255)
                                      aaCCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGt                                 >  1:2097405/1‑65 (MQ=255)
                                      aaCCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCg                            >  1:1635206/1‑70 (MQ=255)
                                      aaCACGCTGGGCTGGGAAATCCTACGTAACTgcgc                                                               >  1:2376084/1‑35 (MQ=255)
                                       aCCCGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCg                            <  1:4131919/69‑1 (MQ=255)
                                       aCCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGc                           <  1:3578044/70‑1 (MQ=255)
                                        cccGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAt                         <  1:3900419/71‑1 (MQ=255)
                                        cccGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAt                         <  1:447309/71‑1 (MQ=255)
                                        cccGCTGGGCTGGGAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAt                         <  1:3412895/71‑1 (MQ=255)
                                        cccGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTc                                                     <  1:972329/43‑1 (MQ=255)
                                        cccGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTc                                                     <  1:4672807/43‑1 (MQ=255)
                                         ccGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAtt                        <  1:167009/71‑1 (MQ=255)
                                         ccGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAt                         <  1:2816502/70‑1 (MQ=255)
                                          cGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCAt                         >  1:2992933/1‑69 (MQ=255)
                                            ctgggctgggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAgg                     <  1:764941/71‑1 (MQ=255)
                                            ctgggctgggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAgg                     <  1:3043408/71‑1 (MQ=255)
                                             tgggctgggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAgg                     <  1:4735290/70‑1 (MQ=255)
                                              gggctgggAAATCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAg                      >  1:939375/1‑68 (MQ=255)
                                              gggctgggAAATCCTACGTAACTGCGCCACCCAATTCa                                                    >  1:2541195/1‑38 (MQ=255)
                                              gggctgggAAATCCTACGTAACTGCGCCACCCAATTCAt                                                   >  1:4861570/1‑39 (MQ=255)
                                                  tgggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGcag                                          >  1:4248010/1‑42 (MQ=255)
                                                  tgggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAg                                         >  1:510704/1‑45 (MQ=255)
                                                   gggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATgcgc              >  1:2228707/1‑71 (MQ=255)
                                                   gggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATgcgc              >  1:3947261/1‑71 (MQ=255)
                                                    ggAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGcgcg             >  1:1611359/1‑71 (MQ=255)
                                                      aaaTCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAgg                     >  1:2893701/1‑61 (MQ=255)
                                                      aaaTCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTg           >  1:943022/1‑71 (MQ=255)
                                                        aTCCTACGTAATTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCt         <  1:2455454/71‑1 (MQ=255)
                                                        aTCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCt         >  1:994824/1‑71 (MQ=255)
                                                         tCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATg                 <  1:2661246/62‑1 (MQ=255)
                                                         tCCTACGAAATTGCGCCACCCAATTCATCTCCTGcc                                           >  1:4757669/1‑36 (MQ=37)
                                                          ccTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTgg       >  1:2816852/1‑71 (MQ=255)
                                                           cTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGt            <  1:2328541/65‑1 (MQ=255)
                                                           cTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGt            <  1:3224952/65‑1 (MQ=255)
                                                              cGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTga                              >  1:1262803/1‑43 (MQ=255)
                                                                tAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTatt  <  1:4785004/70‑3 (MQ=255)
                                                                   |                                                                  
TGGCCACCATCATGGCAGGCCTGGCCTGTGGCGAACCTAACCCGCTGGGCTGGGAAATCCTACGTAACTGCGCCACCCAATTCATCTCCTGCCAGGACAGCGTTGCCGCATTAGGTATGCGCGTGCTGGGTAAT  >  W3110S.gb/3007010‑3007143

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: