Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,161,757 G→A 8.5% T127T (ACC→ACT plsC ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,161,7570GA 8.5% 142.4 / 3.8 47T127T (ACC→ACTplsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (2/2);  total (22/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTACAACT  >  W3110S.gb/3161687‑3161821
                                                                      |                                                                
cgCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATg                                                                  >  1:1732235/2‑71 (MQ=255)
     gTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGcc                                                             >  1:1198887/1‑71 (MQ=255)
        ccttccGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGt                                                           >  1:3317849/1‑70 (MQ=255)
         cttccGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTg                                                          <  1:523627/70‑1 (MQ=255)
          ttccGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAg                                                        >  1:1737930/1‑71 (MQ=255)
           tccGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGc                                                       >  1:1608171/1‑71 (MQ=255)
              gggAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCttt                                                    <  1:184408/71‑1 (MQ=255)
              gggAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCtt                                                     >  1:249508/1‑70 (MQ=255)
               ggAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTa                                                   <  1:2946945/71‑1 (MQ=255)
                gAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTa                                                   <  1:2850196/70‑1 (MQ=255)
                gAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATAGTGCCGTGAGCTTTAg                                                  >  1:2989790/1‑71 (MQ=255)
                     tCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGt                                                           <  1:886607/57‑1 (MQ=255)
                     tCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTACCGCAATGGTGCCGTGAGCTTTAGTGCGa                                             >  1:3688720/1‑71 (MQ=255)
                         gATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGt                                         <  1:1744033/71‑1 (MQ=255)
                         gATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGt                                         >  1:717948/1‑71 (MQ=255)
                         gATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGt                                         <  1:3649242/71‑1 (MQ=255)
                                 tACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCg                                 <  1:4172644/71‑1 (MQ=255)
                                    ggCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGa                                >  1:3351819/1‑69 (MQ=255)
                                     gCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATc                              >  1:4065067/1‑70 (MQ=255)
                                      cGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCaa                            <  1:1483999/71‑1 (MQ=255)
                                        ttttttGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTc                                  <  1:870381/63‑1 (MQ=255)
                                        ttttttGAAGTGATTCACTACTTCAGCAATGGTGCCGTGAGCtt                                                     <  1:837740/44‑1 (MQ=255)
                                         tttttGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAAt                           >  1:2158071/1‑69 (MQ=255)
                                          ttttGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAg                        >  1:1118981/1‑71 (MQ=255)
                                          ttttGAAGTGATTCACTACTTCCGCAATAGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAg                        >  1:4229211/1‑71 (MQ=255)
                                           tttGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAg                                                  <  1:3862379/44‑1 (MQ=255)
                                               aaGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGc                    <  1:1908811/70‑1 (MQ=255)
                                                 gTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCgg                 <  1:1034027/71‑1 (MQ=255)
                                                 gTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCgg                 >  1:1906842/1‑71 (MQ=255)
                                                  tGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGAt                               <  1:1294050/56‑1 (MQ=255)
                                                  tGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCgg                 <  1:4035691/70‑1 (MQ=255)
                                                   gATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGa                                             >  1:3155147/1‑41 (MQ=255)
                                                   gATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCgg                 >  1:2253569/1‑69 (MQ=255)
                                                      tCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAAc            >  1:2369790/1‑71 (MQ=255)
                                                       cACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAAc            >  1:1029370/1‑70 (MQ=255)
                                                        actactTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATaa                         >  1:3914047/1‑56 (MQ=255)
                                                            ctTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGttt                                       >  1:760287/1‑38 (MQ=255)
                                                            ctTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGttt                                       >  1:1332935/1‑38 (MQ=255)
                                                            ctTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGttt                                       >  1:651856/1‑38 (MQ=255)
                                                            ctTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTAc      <  1:3858859/71‑1 (MQ=255)
                                                             ttCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATa                          <  1:75559/50‑1 (MQ=255)
                                                               ccGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTACaa    <  1:100428/70‑1 (MQ=255)
                                                                cGCAATAGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTACAACt  <  1:1803263/71‑1 (MQ=255)
                                                                cGCAATAGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTACAACt  <  1:1808176/71‑1 (MQ=255)
                                                                  cAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAAt                           <  1:3234698/44‑1 (MQ=255)
                                                                     tGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGtt               <  1:3370026/53‑1 (MQ=255)
                                                                     tGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGtt               <  1:1423288/53‑1 (MQ=255)
                                                                     tGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCa          <  1:2783724/58‑1 (MQ=255)
                                                                      |                                                                
TGCGGGTTCCTTCCGGGAACATCCAGATGGAAATACGGCGTTTTTTGAAGTGATTCACTACTTCCGCAATGGTGCCGTGAGCTTTAGTGCGATTGTTTCTGTCGATCAATAAGTTGCCGGTTAACCAGTACAACT  >  W3110S.gb/3161687‑3161821

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: