Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,970,934 C→T 8.0% E658K (GAG→AAG)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,970,9340CT 8.0% 233.0 / 8.3 75E658K (GAG→AAG) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (26/43);  new base T (3/3);  total (29/46)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

TGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTT  >  Exported/1970867‑1970999
                                                                   |                                                                 
tGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCAt                                                                 <  1:3716046/70‑1 (MQ=255)
cgccacaATTTATTGTGACGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATc                                                                >  1:3671718/4‑71 (MQ=255)
  ccacaATTTATTGTGACGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCt                                                               >  1:200418/2‑70 (MQ=255)
    acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                             >  1:1947603/1‑70 (MQ=255)
    acaATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                             >  1:3052194/1‑70 (MQ=255)
        tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGa                                                         >  1:2152782/1‑70 (MQ=255)
        tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                        <  1:2589048/71‑1 (MQ=255)
        tAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGAt                                                        <  1:3031234/71‑1 (MQ=255)
            ttGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGt                                                    <  1:3711332/71‑1 (MQ=39)
            ttGTGACGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATaa                                                      <  1:425405/69‑1 (MQ=255)
             tgtgGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGta                                                   >  1:1799519/1‑71 (MQ=255)
                 gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGtaagt                                                <  1:106720/70‑4 (MQ=37)
                 gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGATGAGtaagt                                                <  1:2106553/70‑4 (MQ=37)
                 gCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAc                                               <  1:124423/71‑1 (MQ=255)
                  cgtaTTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                              <  1:3028940/68‑1 (MQ=255)
                  cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTc                                                                   <  1:2995851/50‑1 (MQ=255)
                  cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTc                                                                   <  1:1940383/50‑1 (MQ=255)
                  cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                              <  1:3024710/71‑1 (MQ=255)
                  cGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAAcc                                              <  1:1237338/71‑1 (MQ=255)
                   ggATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                             <  1:1411514/71‑1 (MQ=255)
                    gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                             >  1:2502941/1‑70 (MQ=255)
                    gATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCg                                             >  1:1134602/1‑70 (MQ=255)
                             cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                    >  1:195207/1‑70 (MQ=255)
                             cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTg                                   >  1:1680106/1‑71 (MQ=255)
                             cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTg                                   >  1:427020/1‑71 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:2713515/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:3221184/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:2717827/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:3256440/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:3509468/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:3202868/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:3600916/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:1617674/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  <  1:1666465/39‑1 (MQ=255)
                              gCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTc                                                             <  1:247472/44‑1 (MQ=255)
                               ctgaaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                 >  1:2993152/4‑71 (MQ=255)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAt                                                        >  1:2180128/1‑48 (MQ=38)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAt                                                        >  1:932372/1‑48 (MQ=38)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTTATCTTGAAGAt                                                        >  1:2475598/1‑48 (MQ=38)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                 >  1:977673/1‑71 (MQ=255)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                 >  1:1932885/1‑71 (MQ=255)
                               cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                 >  1:2080262/1‑71 (MQ=255)
                                tAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAg                                                     <  1:1605696/50‑1 (MQ=255)
                                 aaaTTAAGCGGATTCAAGTAACACAGGACTTACTc                                                                   >  1:1970180/1‑35 (MQ=255)
                                 aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  >  1:1958656/1‑36 (MQ=255)
                                 aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  >  1:3105445/1‑36 (MQ=255)
                                 aaaTTAAGCGGATTCAAGTAACACAGGACTTACTCa                                                                  >  1:3012636/1‑36 (MQ=255)
                                  aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCt                                <  1:2813879/69‑1 (MQ=255)
                                  aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCt                                <  1:3136359/69‑1 (MQ=255)
                                     tAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAc                           >  1:179689/1‑71 (MQ=255)
                                      aaGCGGATTCAAGTAACACAAGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAAc                           <  1:3215923/70‑1 (MQ=255)
                                         cGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                       >  1:836573/1‑71 (MQ=255)
                                          ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                      <  1:3374487/71‑1 (MQ=255)
                                          ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                      <  1:3426534/71‑1 (MQ=255)
                                          ggATTCAAGTAACACAGGACTTACTCATCTACAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                      <  1:3131050/71‑1 (MQ=255)
                                               cAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCt                                    >  1:126548/1‑52 (MQ=255)
                                                aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTTGAActc                         <  1:480183/62‑1 (MQ=255)
                                                aaGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActc                         <  1:258844/62‑1 (MQ=255)
                                                 aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAActc                         <  1:2566275/61‑1 (MQ=255)
                                                 aGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAAc               <  1:1924486/71‑1 (MQ=255)
                                                  gTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACt              >  1:835745/1‑71 (MQ=255)
                                                     acacAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt           <  1:1073160/71‑1 (MQ=255)
                                                     acacAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt           >  1:63040/1‑71 (MQ=255)
                                                      cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                         <  1:2051402/40‑1 (MQ=255)
                                                      cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                         <  1:2024216/40‑1 (MQ=255)
                                                      cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                         <  1:1393105/40‑1 (MQ=255)
                                                       acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGaa                            <  1:1292070/52‑1 (MQ=255)
                                                       acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttgt         <  1:2133196/71‑1 (MQ=255)
                                                       acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttgt         <  1:1482143/71‑1 (MQ=255)
                                                       acaGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgttgt         <  1:2409955/71‑1 (MQ=255)
                                                        caGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTgttgtt        <  1:411744/71‑1 (MQ=255)
                                                         aGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTgttgtt        >  1:2504112/1‑70 (MQ=255)
                                                          ggACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg       >  1:62342/1‑70 (MQ=255)
                                                            aCTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                       <  1:828026/52‑1 (MQ=255)
                                                            aCTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCtt                       <  1:872425/52‑1 (MQ=255)
                                                            aCTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAg    <  1:2913609/71‑1 (MQ=255)
                                                             cTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGt   >  1:151927/1‑71 (MQ=255)
                                                                aCTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGtt  <  1:1456861/69‑1 (MQ=255)
                                                                   |                                                                 
TGTCACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTT  >  Exported/1970867‑1970999

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: