Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,799,263 G→A 31.4% S132S (AGC→AGT pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,799,2630GA31.4% 46.7 / 25.4 35S132S (AGC→AGTpheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base G (12/12);  new base A (3/8);  total (15/20)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGACGGTTCG  >  W3110S.gb/1799198‑1799323
                                                                 |                                                            
tcaAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACtgtg                                                         <  1:2735693/71‑1 (MQ=255)
 caAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGa                                                        >  1:1556758/1‑71 (MQ=255)
  aaGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGa                                                        >  1:1950340/1‑70 (MQ=255)
  aaGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGAt                                                       >  1:963218/1‑71 (MQ=255)
  aaGTTTAAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTAGATAATGCCGCTGTGAt                                                       >  1:4083/1‑71 (MQ=255)
      ttCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGt                                                    >  1:1448551/1‑70 (MQ=255)
         aGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAg                                                 <  1:1363595/70‑1 (MQ=255)
         aGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGa                                                <  1:1288378/71‑1 (MQ=255)
           gTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGaaa                                              >  1:1430733/1‑71 (MQ=255)
                cACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGcc                                          >  1:1848895/1‑70 (MQ=255)
                        tCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCaga                                 <  1:119274/71‑1 (MQ=255)
                        tCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCag                                  <  1:2470337/70‑1 (MQ=255)
                        tCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCag                                  <  1:530644/70‑1 (MQ=255)
                        tCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATg                                            >  1:1575565/1‑60 (MQ=255)
                          ggTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCagag                                <  1:59039/70‑1 (MQ=255)
                          ggTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCagag                                <  1:1311289/70‑1 (MQ=255)
                            tGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCTGAGAAg                             <  1:959/71‑1 (MQ=255)
                            tGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCAGAGAAg                             <  1:627086/71‑1 (MQ=255)
                               cAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAg                          <  1:2000167/71‑1 (MQ=255)
                               cAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAg                          <  1:2734109/71‑1 (MQ=255)
                               cAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAg                          >  1:245007/1‑71 (MQ=255)
                               cAATCGGCGCATCCGCAGGCAGTTCGATAATGCCACTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAg                          <  1:191949/71‑1 (MQ=255)
                                        cATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATc                 <  1:57090/71‑1 (MQ=255)
                                          tCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCag                                  <  1:450078/52‑1 (MQ=255)
                                              caggcagTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATccc               >  1:727266/1‑67 (MQ=255)
                                              caggcagTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATccc               >  1:837203/1‑67 (MQ=255)
                                                ggcagTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAgg                            <  1:2343227/52‑1 (MQ=255)
                                                 gcagTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGAc        <  1:1419619/71‑1 (MQ=255)
                                                   agTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCAt                  <  1:2731065/59‑1 (MQ=255)
                                                   agTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCAt                  <  1:758001/59‑1 (MQ=255)
                                                   agTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCAt                  <  1:2037615/59‑1 (MQ=255)
                                                       cGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGACGGTTCg  >  1:64815/1‑71 (MQ=255)
                                                       cGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGACGGTTCg  >  1:727706/1‑71 (MQ=255)
                                                              gCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGca                        >  1:2597228/1‑42 (MQ=255)
                                                               ccGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGa         >  1:35615/1‑56 (MQ=255)
                                                                 |                                                            
TCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAGTTCGATAATGCCGCTGTGATCGTCAGAAATGCCCAGTTCAGAGAAGGAGCACAGCATCCCTTCAGACGGTTCG  >  W3110S.gb/1799198‑1799323

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: