Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 413,316 A→G 100% R289R (CGT→CGC lnt ← apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE413,3160AG79.3% 30.6 / 2.1 24R289R (CGT→CGClntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (1/18);  total (3/21)
Fisher's exact test for biased strand distribution p-value = 9.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTATCGTAGCGGTTCTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGAACGGTTGCT  >  minE/413248‑413353
                                                                    |                                     
gTATCGTAGCGGTTCTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCa                                       >  1:692732/1‑69 (MQ=255)
 tATCGTAGCGGTTCTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCa                                       >  1:2505720/1‑68 (MQ=255)
 tATCGTAGCGGTTCTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCa                                       >  1:3129310/1‑68 (MQ=255)
  aTCGTAGCGGTTCTGTTTATTGAGACGCGCGTCGACAATCCCGG‑‑ACCAGCGAGCTACCTTTATCGCGCa                                   >  1:748042/1‑69 (MQ=255)
              cTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTc                         <  1:363713/69‑1 (MQ=255)
                            gcgcgTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGa           >  1:2018344/1‑69 (MQ=255)
                            gcgcgTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGa           >  1:2294955/1‑69 (MQ=255)
                               cgTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGAACg        <  1:1415370/69‑1 (MQ=255)
                               cgTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGAACg        <  1:2089790/69‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCTCCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:2648144/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:919501/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1032285/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:790791/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:344507/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:3430976/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:3419718/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:2682686/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:2006545/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1944563/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1907775/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1771658/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1414616/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1367520/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1121884/47‑1 (MQ=255)
                                                           aCCTTTATCGCGCAACTCACAGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:1107102/47‑1 (MQ=255)
                                                             cTTTATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:525198/45‑1 (MQ=255)
                                                                tATCGCGCAACTCACCGTCCAGAGCTTTGAGGAACGGTTGCt  <  1:395650/42‑1 (MQ=255)
                                                                    |                                     
GTATCGTAGCGGTTCTGCTTATTGAGACGCGCGTCGACAATCCCGGTTACCAGCGAGCTACCTTTATCACGCAACTCACCGTCCAGTGCTTTGAGGAACGGTTGCT  >  minE/413248‑413353

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: