Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 940,557 A→G 14.6% P236P (CCA→CCG serS → seryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb940,5570AG14.6% 118.5 / 13.4 48P236P (CCA→CCGserSseryl‑tRNA synthetase, also charges selenocysteinyl‑tRNA with serine
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/5);  total (22/26)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCCAATTAA  >  W3110S.gb/940498‑940625
                                                           |                                                                    
cTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGAAGt                                                           >  1:215361/1‑71 (MQ=255)
   ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCGACGGCAGa                                                              <  1:1145974/65‑1 (MQ=255)
   ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGa                                                              <  1:1114790/65‑1 (MQ=255)
   ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGa                                                              <  1:302852/65‑1 (MQ=255)
   ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGa                                                              <  1:604674/65‑1 (MQ=255)
   ttCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAg                                                            <  1:1735626/67‑1 (MQ=255)
     ccATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGa                                                              <  1:3644259/63‑1 (MQ=255)
     ccATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGc                                                      >  1:550943/1‑71 (MQ=255)
      cATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCt                                                     <  1:2156292/71‑1 (MQ=255)
          cTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGaa                                                             >  1:705876/1‑59 (MQ=255)
          cTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGa                                                   >  1:1705939/1‑69 (MQ=255)
                cGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTg                                           <  1:69543/71‑1 (MQ=255)
                        gaagaagCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCgg                                   >  1:1183979/1‑71 (MQ=255)
                            aagCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCt                                            <  1:1547136/58‑1 (MQ=255)
                              gCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAat                             >  1:3538923/1‑71 (MQ=255)
                                 gACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGtt                                                          >  1:1244965/1‑39 (MQ=255)
                                 gACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAatcat                          >  1:1524440/1‑71 (MQ=255)
                                   cacCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCg                        <  1:896460/71‑1 (MQ=255)
                                   cacCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCg                        <  1:931813/71‑1 (MQ=255)
                                    acCAGTAACTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGa                       >  1:2347720/1‑71 (MQ=255)
                                      cAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGAt                      >  1:3849752/1‑70 (MQ=255)
                                      cAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATg                     <  1:1689290/71‑1 (MQ=255)
                                       aGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGc                                                      >  1:4016986/1‑37 (MQ=255)
                                       aGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGc                                                      >  1:293511/1‑37 (MQ=255)
                                       aGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGc                                                      >  1:2625415/1‑37 (MQ=255)
                                       aGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGc                                                      >  1:2620765/1‑37 (MQ=255)
                                       aGTAACTATGCGCTGATCACAACGGCAGAAGTTCCGc                                                      >  1:2921013/1‑37 (MQ=255)
                                        gTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCt                                            >  1:1831341/1‑46 (MQ=255)
                                         tAACTATGCGCTGATCCCGACTGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAat                             <  1:3123389/60‑1 (MQ=255)
                                           aCTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAgat                >  1:2170742/1‑71 (MQ=255)
                                           aCTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAga                 >  1:491005/1‑70 (MQ=255)
                                            cTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                                <  1:2792783/54‑1 (MQ=255)
                                            cTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                                <  1:775885/54‑1 (MQ=255)
                                            cTATGCGCTGATCCCGACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGa                                <  1:2017663/54‑1 (MQ=255)
                                            cTATGCGCTGATCCCAACGTCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAgatg               <  1:2206087/71‑1 (MQ=255)
                                            cTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAgatg               <  1:2722016/71‑1 (MQ=255)
                                               tGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAAcc                                             <  1:1023324/38‑1 (MQ=255)
                                               tGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAAcc                                             <  1:3305164/38‑1 (MQ=255)
                                               tGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAgatga              <  1:601422/69‑1 (MQ=255)
                                               tGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATc            >  1:3880736/1‑71 (MQ=255)
                                               tGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATc            >  1:12865/1‑71 (MQ=255)
                                                 cgcTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTg                                           >  1:3671475/1‑38 (MQ=255)
                                                      aTCCCAACGTCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAatcatc                         <  1:1320358/51‑1 (MQ=255)
                                                      aTCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAatcatc                         <  1:3636536/51‑1 (MQ=255)
                                                        cccAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCCAAtt    >  1:1893965/1‑70 (MQ=255)
                                                         ccAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCGAATTaa  <  1:8204/71‑1 (MQ=255)
                                                         ccAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCCAATTaa  <  1:2284597/71‑1 (MQ=255)
                                                         ccAACGCCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCCAATTaa  <  1:628535/71‑1 (MQ=255)
                                                           |                                                                    
CTGTTCCATACTCGTCCGCTGGAAGAAGAAGCAGACACCAGTAACTATGCGCTGATCCCAACGGCAGAAGTTCCGCTGACTAACCTGGTGCGCGGTGAAATCATCGATGAAGATGATCTGCCAATTAA  >  W3110S.gb/940498‑940625

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: