Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 2,972,729 T→C 13.1% S648G (AGC→GGC)  aas ← fused 2‑acylglycerophospho‑ethanolamine acyl transferase and acyl‑acyl carrier protein synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090482,972,7290TC13.1% 155.4 / 14.2 61S648G (AGC→GGC) aasfused 2‑acylglycerophospho‑ethanolamine acyl transferase and acyl‑acyl carrier protein synthetase
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (5/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGA  >  AP009048/2972660‑2972798
                                                                     |                                                                     
gCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGct                                                                       <  1:1064996/70‑1 (MQ=255)
  aaCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGccct                                                                     >  1:3676215/1‑70 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                    >  1:1858814/1‑70 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                    >  1:776767/1‑70 (MQ=255)
   aCTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCtt                                                                    >  1:782154/1‑70 (MQ=255)
      tATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGccct                                                                     >  1:3406812/1‑66 (MQ=255)
      tATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGccct                                                                     >  1:1663337/1‑66 (MQ=255)
       aTCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAAt                                                                 <  1:401480/69‑1 (MQ=255)
        tCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCg                                                               <  1:2097213/70‑1 (MQ=255)
        tCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCg                                                               >  1:2945268/1‑70 (MQ=255)
        tCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGc                                                              <  1:2876330/71‑1 (MQ=255)
        tCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCg                                                               <  1:4030766/70‑1 (MQ=255)
           cgcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGt                                                           <  1:3652117/71‑1 (MQ=255)
            gcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTg                                                          <  1:4598974/71‑1 (MQ=255)
            gcgTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTg                                                          >  1:1690467/1‑71 (MQ=255)
                cAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGc                                                              >  1:1754790/1‑63 (MQ=255)
                cAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCAt                                                      >  1:100891/1‑71 (MQ=255)
                 aGTTCGTTATCTGTGTTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCAt                                                      <  1:3401501/70‑1 (MQ=255)
                  gTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATgg                                                    <  1:676157/71‑1 (MQ=255)
                   ttCGTTATCTGTGGTGAAAAGCACGAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGc                                                              <  1:44254/60‑1 (MQ=255)
                   ttCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCa                                                       >  1:3659348/1‑67 (MQ=255)
                   ttCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGa                                                   >  1:2815956/1‑71 (MQ=255)
                    tCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTg                                                          >  1:974871/1‑63 (MQ=255)
                     cGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACt                                                 >  1:236035/1‑71 (MQ=255)
                     cGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACt                                                 >  1:2471598/1‑71 (MQ=255)
                     cGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACt                                                 >  1:3774635/1‑71 (MQ=255)
                       ttATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGctct                                                                     <  1:4391534/49‑1 (MQ=255)
                         aTCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTg                                                          >  1:350653/1‑58 (MQ=255)
                         aTCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTAt                                             >  1:3064893/1‑71 (MQ=255)
                         aTCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTAt                                             >  1:3215610/1‑71 (MQ=255)
                          tCTGTGGTGAAGCGCACCAGTGCCTCGCCTTTACTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATc                                            >  1:2535308/1‑71 (MQ=255)
                              tggtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGAATGGACTTTATCt                                           <  1:4443923/68‑1 (MQ=255)
                                gtgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGa                                      <  1:1461557/71‑1 (MQ=255)
                                 tgAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGaa                                     <  1:934594/71‑1 (MQ=255)
                                  gAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATg                                                     <  1:981673/54‑1 (MQ=255)
                                  gAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                    >  1:1774080/1‑71 (MQ=255)
                                   aaaaGCACCAGTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTAt                                             >  1:2578001/1‑61 (MQ=255)
                                      aGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCt                                           <  1:1956455/60‑1 (MQ=255)
                                        cACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCaa                              <  1:3007208/71‑1 (MQ=255)
                                           cAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGa                                                   <  1:2837844/47‑1 (MQ=255)
                                           cAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTg                           >  1:3492050/1‑71 (MQ=255)
                                            aGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGa                                      >  1:4548998/1‑59 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAaca                                  >  1:2925521/1‑62 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCAGGGACTTTATCTgg                                         >  1:2194087/1‑55 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                    >  1:3345922/1‑60 (MQ=255)
                                             gTGCCTCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                    >  1:1562712/1‑60 (MQ=255)
                                              tGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCa                        <  1:3476092/71‑1 (MQ=255)
                                                ccTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAg                       <  1:2602484/70‑1 (MQ=255)
                                                ccTCGCCCTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTa                                              >  1:2452517/1‑47 (MQ=255)
                                                  tCGCCTTTGCTGGCATCGCCCTTAATCGCAGTGGCATGGACtt                                                <  1:1214671/43‑1 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGTCAgttgt                   >  1:796943/1‑71 (MQ=255)
                                                   cGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgt                      <  1:3308411/68‑1 (MQ=255)
                                                    gCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                  >  1:1675078/1‑71 (MQ=255)
                                                    gCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                  <  1:3223346/71‑1 (MQ=255)
                                                    gCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                  >  1:1341735/1‑71 (MQ=255)
                                                    gCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAgttgtt                  <  1:1042812/71‑1 (MQ=255)
                                                     ccTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTc                 <  1:1171417/71‑1 (MQ=255)
                                                      cTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTc                 >  1:626503/1‑70 (MQ=255)
                                                         tGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTTTTccac              <  1:10428/70‑1 (MQ=255)
                                                         tGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTccac              >  1:2524084/1‑70 (MQ=255)
                                                              gCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTccac              >  1:2871461/1‑65 (MQ=255)
                                                                  cGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAg     <  1:4330234/70‑1 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                    >  1:1094450/1‑38 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGaaa                                    >  1:1019214/1‑38 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGc    >  1:3083458/1‑70 (MQ=255)
                                                                   gCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCg   <  1:3736554/71‑1 (MQ=255)
                                                                    ctctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGc                          <  1:1253193/47‑1 (MQ=255)
                                                                    ctctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCAtt          <  1:3129551/63‑1 (MQ=255)
                                                                    ctctTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGa  >  1:3119504/1‑71 (MQ=255)
                                                                     |                                                                     
GCAACTTATCGCGCGTCAGTTCGTTATCTGTGGTGAAAAGCACCAGTGCCTCGCCTTTGCTGGCATCGCTCTTAATCGCAGTGGCATGGACTTTATCTGGCGAAACACCAAGTGCCAGTTGTTCCACCATTTCCAGCGA  >  AP009048/2972660‑2972798

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: