Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 413,754 C→T 40.9% L143L (CTG→CTA lnt ← apolipoprotein N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE413,7540CT40.9% 4.2 / 11.4 22L143L (CTG→CTAlntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base T (3/6);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ATCGGTGCCAGCCCTTTTAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGCAGGGGCGGCAATCGCTACGCGCAGCCA  >  minE/413686‑413822
                                                                    |                                                                    
aTCGGTGCCAGCCCTTTTAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTc                                                                      >  1:365569/1‑69 (MQ=255)
aTCGGTGCCAGCCCTTTTAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTc                                                                      >  1:2000204/1‑69 (MQ=255)
                ttAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACcgcg                                                        >  1:2211810/1‑67 (MQ=255)
                ttAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCgcgc                                                       >  1:1315239/1‑68 (MQ=255)
                ttAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCgcgc                                                       >  1:3377681/1‑68 (MQ=255)
                 tAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAg                                                                    <  1:235062/54‑1 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGa                                                    <  1:349955/58‑1 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:1126059/1‑69 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:808474/1‑69 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:782014/1‑69 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:405499/1‑69 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:2049875/1‑69 (MQ=255)
                             aaTCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCAc                                         >  1:3092122/1‑69 (MQ=255)
                               tCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACtt                                       <  1:2753943/69‑1 (MQ=255)
                                   gCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACt                                        <  1:3050283/64‑1 (MQ=255)
                                   gCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACt                                        <  1:1629706/64‑1 (MQ=255)
                                     tataGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:1133920/67‑1 (MQ=255)
                                     tataGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:1582003/67‑1 (MQ=255)
                                     tataGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:1472601/67‑1 (MQ=255)
                                     tataGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:1360923/67‑1 (MQ=255)
                                     tataGCCGAACTGTAACCACGGGAAGCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:1252460/67‑1 (MQ=255)
                                         gCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGc                             >  1:851023/1‑69 (MQ=255)
                                                        cGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGCAGGGGCGGCAATCGc              <  1:1582643/69‑1 (MQ=255)
                                                              gCCTGTTAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCa                                   <  1:2731093/42‑1 (MQ=38)
                                                                    cAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGCAGGGGCGGCAATCGCTACGCGCAGCCa  >  1:2323385/1‑69 (MQ=255)
                                                                    cAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGCAGGGGCGGCAATCGCTACGCGCAGCCa  >  1:864879/1‑69 (MQ=255)
                                                                    |                                                                    
ATCGGTGCCAGCCCTTTTAACGGGCCATCAATCTGGCTATAGCCGAACTGTAACCACGGGAAGCCGGTCAGTACCCAACCGCGCAGAAACTCGGTCACTTGCCAGAGGGCAGGGGCGGCAATCGCTACGCGCAGCCA  >  minE/413686‑413822

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: