Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 189,664 G→A 21.6% E387E (GAG→GAA)  ldcC → lysine decarboxylase 2, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported189,6640GA21.6% 79.1 / 14.6 37E387E (GAG→GAA) ldcClysine decarboxylase 2, constitutive
Reads supporting (aligned to +/- strand):  ref base G (13/16);  new base A (3/5);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGT  >  Exported/189595‑189729
                                                                      |                                                                 
gACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAgg                                                                  >  1:969445/1‑71 (MQ=255)
     acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                                             <  1:1393737/71‑1 (MQ=255)
     acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCttt                                                             <  1:5283223/71‑1 (MQ=255)
     acaAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCtt                                                              >  1:165013/1‑70 (MQ=255)
      caAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTa                                                            <  1:5347006/71‑1 (MQ=255)
        aactGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAAc                                                          <  1:3267072/68‑1 (MQ=255)
         actGCTGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAAc                                                          >  1:2664385/3‑70 (MQ=255)
          aTGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGa                                                        <  1:1993556/71‑1 (MQ=255)
              tGGCGGCGTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGcc                                                    >  1:5052938/1‑71 (MQ=9)
                gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCta                                                  <  1:1625746/71‑2 (MQ=9)
                gcggcgTTCTCTCAGGCTTCCATGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCta                                                  <  1:3187422/71‑2 (MQ=9)
                 cggcgTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCttt                                                 <  1:3467769/71‑1 (MQ=255)
                     gTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatga                                             <  1:1687361/71‑1 (MQ=255)
                        aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                                           >  1:200255/1‑70 (MQ=255)
                        aTCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatgatg                                           >  1:5344576/1‑70 (MQ=255)
                           gCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTat                                               >  1:2311735/1‑63 (MQ=255)
                             aGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                                                   >  1:810037/1‑41 (MQ=255)
                              ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACg                                                         >  1:4697611/1‑50 (MQ=255)
                              ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGa                                                        >  1:3355102/1‑51 (MQ=255)
                              ggCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATacca                                    >  1:5404280/1‑71 (MQ=255)
                                 ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAagag                                                                   >  1:260674/1‑37 (MQ=255)
                                 ttCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATac                                      >  1:1271987/1‑66 (MQ=255)
                                    gCTGATCCACAATAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccaccacc                               <  1:4619047/70‑1 (MQ=255)
                                     cTGATCCACATTAAAGGCGA‑TTATGACGAAGAGGCCTTTAACGAAGCCTTTatgat                                            <  1:224533/56‑1 (MQ=255)
                                     atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                             <  1:4482478/70‑1 (MQ=9)
                                     atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                             <  1:4488805/70‑1 (MQ=9)
                                     atGATCCACGTTAAAGGTGACGTA‑AACGAAGAAACCTTTAACGAAGCCTACATGATGCACACCACCACTTc                             >  1:3772016/2‑71 (MQ=9)
                                         tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGccc                         <  1:5266652/71‑1 (MQ=255)
                                         tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGccc                         <  1:3725837/71‑1 (MQ=255)
                                         tCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGcc                          <  1:1177215/70‑1 (MQ=255)
                                          ccACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTatga                                             >  1:4814425/1‑50 (MQ=255)
                                               ttAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGc                                                     <  1:5376660/37‑1 (MQ=255)
                                                 aaaGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATc                  >  1:5594564/1‑70 (MQ=255)
                                                     gCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAtt             >  1:3416724/1‑71 (MQ=255)
                                                      cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAttg            <  1:3962792/71‑1 (MQ=255)
                                                      cGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAtt             <  1:3617773/70‑1 (MQ=255)
                                                       gA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATaccacca                                 >  1:3986262/1‑49 (MQ=255)
                                                               aCGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCAt              >  1:2169557/1‑61 (MQ=255)
                                                               aCGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTc     <  1:4718442/70‑1 (MQ=255)
                                                                 gAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGt  <  1:2283506/71‑1 (MQ=255)
                                                                      |                                                                 
GACCCACAAAATGCTGGCGGCGTTATCGCAGGCTTCGCTGATCCACATTAAAGGCGA‑GTATGACGAAGAGGCCTTTAACGAAGCCTTTATGATGCATACCACCACCTCGCCCAGTTATCCCATTGTTGCTTCGGT  >  Exported/189595‑189729

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