Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,501,495 G→A 100% S135S (TCG→TCA ddg → palmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,501,4950GA87.5% 79.3 / 5.0 32S135S (TCG→TCAddgpalmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (9/19);  total (10/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACATAAT  >  W3110S.gb/2501430‑2501561
                                                                 |                                                                  
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTgg                                                                <  1:1703996/70‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTgg                                                                <  1:2754895/70‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                               <  1:902784/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                               <  1:1719010/71‑1 (MQ=255)
    aCCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACt                                                            <  1:3570427/70‑1 (MQ=255)
       ttAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGggc                                                        >  1:3423192/1‑71 (MQ=255)
         aaaCGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGggcg                                                       <  1:83082/70‑1 (MQ=255)
              cgcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACt                                                            <  1:1038624/60‑1 (MQ=255)
               gcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGt                                                <  1:942785/71‑1 (MQ=255)
                   aaaTGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATg                                            >  1:1000736/1‑71 (MQ=255)
                                    gTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat                           <  1:3277750/71‑1 (MQ=255)
                                     tAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat                           <  1:3303343/70‑1 (MQ=255)
                                       aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat                        >  1:2628096/1‑71 (MQ=255)
                                       aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat                        >  1:749727/1‑71 (MQ=255)
                                       aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatgat                        >  1:1715156/1‑71 (MQ=255)
                                       aTGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAatga                         >  1:1643787/1‑70 (MQ=255)
                                       aTGGATGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGGCAACCAatgct                        >  1:708437/1‑69 (MQ=255)
                                             gTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTa                   <  1:1668088/70‑1 (MQ=255)
                                                  cGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATgg                                           <  1:2670228/41‑1 (MQ=255)
                                                   gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat                           <  1:936895/56‑1 (MQ=255)
                                                   gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat                           <  1:2483506/56‑1 (MQ=255)
                                                   gTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAat                           <  1:1721976/56‑1 (MQ=255)
                                                    tCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGccaacc                              <  1:2034521/52‑1 (MQ=255)
                                                    tCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGccaacc                              <  1:258681/52‑1 (MQ=255)
                                                      cATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTc          >  1:2976872/1‑70 (MQ=255)
                                                      cATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTc          >  1:3704029/1‑70 (MQ=255)
                                                        ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg                                           <  1:2037899/35‑1 (MQ=255)
                                                        ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg                                           <  1:1783406/35‑1 (MQ=255)
                                                        ttttATGTCACTGGAACTGGGCGGCCGCGTGATgg                                           <  1:672379/35‑1 (MQ=255)
                                                          ttATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCcaca      >  1:2980910/1‑70 (MQ=255)
                                                             tGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACataat  <  1:379510/71‑1 (MQ=255)
                                                             tGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACataat  <  1:2908696/71‑1 (MQ=255)
                                                                 |                                                                  
GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACATAAT  >  W3110S.gb/2501430‑2501561

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: