Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 189,496 T→C 10.4% S331S (TCT→TCC ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported189,4960TC10.4% 175.9 / 9.1 67S331S (TCT→TCCldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base C (5/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCG  >  Exported/189427‑189561
                                                                     |                                                                  
tgatgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTg                                                                   >  1:5642568/1‑71 (MQ=255)
  atgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGcc                                                                 >  1:1373105/1‑71 (MQ=255)
   tgGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCt                                                                <  1:4418760/71‑1 (MQ=255)
     gCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTg                                                               <  1:614753/70‑1 (MQ=255)
       ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            >  1:2428289/1‑71 (MQ=255)
       ttGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            >  1:2674363/1‑71 (MQ=255)
        tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:69485/1‑71 (MQ=255)
        tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:636499/1‑71 (MQ=255)
        tGCTGTACAACACCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTg                                                           >  1:130997/1‑71 (MQ=255)
         gCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGc                                                          >  1:2076418/1‑71 (MQ=255)
          ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                         <  1:5310417/71‑1 (MQ=255)
          ctctACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGcc                                                         >  1:2328671/1‑71 (MQ=255)
            ctACAACACCGACTGGATCAAACAGCCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCg                                                        >  1:2405386/1‑70 (MQ=255)
               caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                    >  1:19933/1‑71 (MQ=255)
               caacaCCGACTTCATCAAGAAAACACTGGATGTGAAATCCATCCACTTTGACTCCGCGTGGGTGCCTTaca                                                    >  1:2258878/1‑71 (MQ=255)
               caacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                    >  1:485429/1‑71 (MQ=255)
                aacaCCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                   >  1:2438796/1‑71 (MQ=255)
                 acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:3466722/61‑1 (MQ=255)
                 acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:302618/61‑1 (MQ=255)
                 acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:4547776/61‑1 (MQ=255)
                 acacCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGt                                                            <  1:5215632/61‑1 (MQ=255)
                  cacCGACTGGACCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACAccc                                                 <  1:5215005/71‑1 (MQ=255)
                   acCGTCTGGATCAAAGAGACGCTGGATGTCCCGTCGATTGACTTCGATTGTGCCTGGGTGCCGTACACCCa                                                <  1:5465383/71‑1 (MQ=255)
                    ccGACTGGAACAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGc                                                                  >  1:770251/1‑52 (MQ=255)
                     cGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTAGATTCTGCCTGGGTGCCGTACACCCAtt                                             >  1:7307/1‑71 (MQ=255)
                              cAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGa                                    >  1:3387964/1‑71 (MQ=255)
                               aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:4588438/1‑45 (MQ=255)
                               aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:801927/1‑45 (MQ=255)
                               aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:2494008/1‑45 (MQ=255)
                               aaaCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTgg                                                              >  1:2431789/1‑45 (MQ=255)
                                 aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTacac                                                   >  1:2049515/1‑54 (MQ=255)
                                 aCAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTaca                                                    >  1:4164848/1‑53 (MQ=255)
                                  cAGACGCTGTATGTCCCGTCGATTCACTTCTATTCTg                                                                   <  1:5132272/37‑1 (MQ=38)
                                  cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                >  1:1274597/1‑71 (MQ=255)
                                  cAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTa                                <  1:3199739/71‑1 (MQ=255)
                                     aCGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:4213022/1‑65 (MQ=255)
                                      cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:1872445/70‑1 (MQ=255)
                                      cGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                            >  1:317984/1‑71 (MQ=255)
                                           gATGTCCCGTCGATTCACTTTGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:5687877/71‑1 (MQ=255)
                                           gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                       <  1:3289155/71‑1 (MQ=255)
                                           gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCTATCTACCAGGGTaa                       <  1:5057101/71‑1 (MQ=255)
                                           gATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:4871286/71‑1 (MQ=255)
                                             tGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaaaa                     <  1:862260/71‑1 (MQ=255)
                                               tCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAg                    <  1:1131259/70‑1 (MQ=255)
                                                  cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:4289751/58‑1 (MQ=255)
                                                  cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:4283153/58‑1 (MQ=255)
                                                  cGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:684431/58‑1 (MQ=255)
                                                    tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:4846197/55‑1 (MQ=255)
                                                    tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:1856686/55‑1 (MQ=255)
                                                    tCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAcc                              <  1:359900/55‑1 (MQ=255)
                                                      gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTa               >  1:4615643/1‑68 (MQ=255)
                                                      gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt              >  1:792084/1‑69 (MQ=255)
                                                      gATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTAt              >  1:2351111/1‑69 (MQ=255)
                                                      gATTCACTTCGATTCTGCCTAGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGa            <  1:4226539/71‑1 (MQ=255)
                                                         tCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAgcg         >  1:2957342/1‑71 (MQ=255)
                                                          cACTTTGACTCCGCGTTGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcg         <  1:2040140/70‑1 (MQ=255)
                                                          cACTTTGACTCCGCGTGGGTGCCTTACACCAACTTCTCACCGATTTACGAAGGTAAATGCGGTATGAgcgg        <  1:609115/71‑1 (MQ=255)
                                                           aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:1966890/55‑1 (MQ=255)
                                                           aCTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTaa                       <  1:5018399/55‑1 (MQ=255)
                                                            cTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTAc                               <  1:3561015/46‑1 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:2358286/1‑39 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGAt                                   >  1:3369274/1‑39 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCa                             <  1:2805639/45‑1 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                           >  1:4457879/1‑47 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAgg                           >  1:5009936/1‑47 (MQ=255)
                                                               cGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAg                            >  1:4993203/1‑46 (MQ=255)
                                                                gATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCg  <  1:5686929/71‑1 (MQ=255)
                                                                     |                                                                  
TGATGGCTTGCTCTACAACACCGACTGGATCAAACAGACGCTGGATGTCCCGTCGATTCACTTCGATTCTGCCTGGGTGCCGTACACCCATTTTCATCCGATCTACCAGGGTAAAAGTGGTATGAGCGGCGAGCG  >  Exported/189427‑189561

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: