Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,970,940 C→G 9.8% E656Q (GAA→CAA)  speA ← biosynthetic arginine decarboxylase, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,970,9400CG 9.8% 174.6 / 8.7 61E656Q (GAA→CAA) speAbiosynthetic arginine decarboxylase, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base C (22/33);  new base G (3/3);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATC  >  Exported/1970871‑1971006
                                                                     |                                                                  
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)
              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:427020/1‑71 (MQ=255)
                         cGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTg                                          >  1:1680106/1‑71 (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:2080262/1‑71 (MQ=255)
                           cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                        >  1:977673/1‑71 (MQ=255)
                           cTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGcc                                        >  1:1932885/1‑71 (MQ=255)
                            tAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAg                                                            <  1:1605696/50‑1 (MQ=255)
                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCt                                       <  1:3136359/69‑1 (MQ=255)
                              aaTTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCt                                       <  1:2813879/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:3426534/71‑1 (MQ=255)
                                      ggATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTc                             <  1:3374487/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:63040/1‑71 (MQ=255)
                                                 acacAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACtgtt                  <  1:1073160/71‑1 (MQ=255)
                                                  cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                                <  1:1393105/40‑1 (MQ=255)
                                                  cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                                <  1:2024216/40‑1 (MQ=255)
                                                  cacaGGACTTACTCATCTTCAAGATAAGTATAACCGTACa                                                <  1:2051402/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)
                                                                aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCa    <  1:1006162/70‑1 (MQ=255)
                                                                aTCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCAt   >  1:2040661/1‑71 (MQ=255)
                                                                   ttCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGttgttg              >  1:2163371/1‑57 (MQ=255)
                                                                    tCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATc  >  1:3636801/1‑68 (MQ=255)
                                                                     |                                                                  
ACAATAAATTGTGGCGGATTATCACCGCTAAATTAAGCGGATTCAAGTAACACAGGACTTACTCATCTTCAAGATAAGTATAACCGTACAAACCTGCCTCGAACTCTTCAAGGAACTGTTGTTGCAGTTCAGCATC  >  Exported/1970871‑1971006

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: