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,4950GA85.7% 48.6 / 2.7 21S135S (TCG→TCAddgpalmitoleoyl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (2/16);  total (2/19)
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.38e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACA  >  W3110S.gb/2501430‑2501557
                                                                 |                                                              
gacaaCCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:1589218/68‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:1458185/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:910713/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:591060/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:49370/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:1746517/71‑1 (MQ=255)
gATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGa                                                           <  1:165663/71‑1 (MQ=255)
 aTAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTAGGCGTCCATTTTATGTCACTGGa                                                           <  1:1944823/70‑1 (MQ=255)
               gcACAAATGCAAAATCGCGGCGTAATTGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGt                                            <  1:1644572/71‑1 (MQ=255)
               gcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGt                                            <  1:1911444/71‑1 (MQ=255)
               gcACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGt                                            <  1:908632/71‑1 (MQ=255)
                        cAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGAc                                    <  1:1858277/70‑1 (MQ=255)
                        cAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGAc                                    <  1:376881/70‑1 (MQ=255)
                        cAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGAc                                    <  1:76396/70‑1 (MQ=255)
                        cAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGGGAc                                    <  1:131953/70‑1 (MQ=255)
                         aaaaTCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGTGATGTGAc                                    <  1:156162/69‑1 (MQ=255)
                                           ttGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGt                                            <  1:205942/43‑1 (MQ=255)
                                           ttGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGt                                            <  1:1638171/43‑1 (MQ=255)
                                            tGTCGGCGTCCATTTTATGTCACTGGAACTGGGCGGCCGCGt                                            <  1:224289/42‑1 (MQ=255)
                                                          ttATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCcaca  >  1:174765/1‑70 (MQ=255)
                                                          ttATGTCACTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCcaca  >  1:895434/1‑70 (MQ=255)
                                                                 |                                                              
GATAACCTTAAACGCGCACAAATGCAAAATCGCGGCGTAATGGTTGTCGGCGTCCATTTTATGTCGCTGGAACTGGGCGGCCGCGTGATGGGACTGTGCCAACCAATGATGGCTACCTATCGTCCACA  >  W3110S.gb/2501430‑2501557

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: