Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,258,286 T→A 8.8% T591T (ACA→ACT plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,258,2860TA 8.8% 149.6 / 6.7 57T591T (ACA→ACTplsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (35/17);  new base A (3/2);  total (38/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.70e-01

CGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAG  >  W3110S.gb/4258217‑4258356
                                                                     |                                                                      
cgcAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCTAtt                                                                        <  1:3273060/70‑1 (MQ=255)
 gcAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATtgt                                                                      >  1:2764670/1‑71 (MQ=255)
  cAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATtgtg                                                                     >  1:1472862/1‑71 (MQ=255)
   aaTGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATtgtg                                                                     >  1:1233407/1‑70 (MQ=255)
   aaTGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATtgtg                                                                     >  1:1327706/1‑70 (MQ=255)
   aaTGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTgtgt                                                                    >  1:2279926/1‑71 (MQ=255)
   aaTGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTgtgt                                                                    >  1:1884232/1‑71 (MQ=255)
     tgttgtTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCt                                                                  <  1:1321541/71‑1 (MQ=255)
       ttgttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCttt                                                                >  1:2561773/1‑71 (MQ=255)
       ttgttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCttt                                                                <  1:1845441/71‑1 (MQ=255)
       ttgttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCtt                                                                 >  1:1687133/1‑70 (MQ=255)
       ttgttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCtt                                                                 >  1:171170/1‑70 (MQ=255)
         gtttCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCTATTGTGTCTTtc                                                               <  1:1955105/70‑1 (MQ=255)
         gttgCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTtct                                                              >  1:1957676/1‑71 (MQ=255)
            gCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCg                                                            <  1:3093492/70‑1 (MQ=255)
            gCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGa                                                           >  1:1744883/1‑71 (MQ=255)
               aTAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACtt                                                        <  1:1364721/71‑1 (MQ=255)
                  gtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCa                                                      >  1:2032116/1‑70 (MQ=255)
                  gtagGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCa                                                      >  1:1223453/1‑70 (MQ=255)
                     ggTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTTGACTTCaaa                                                    >  1:809772/1‑69 (MQ=255)
                     ggTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAAc                                                   <  1:1592174/70‑1 (MQ=255)
                      gTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTtc                                                               <  1:1391928/57‑1 (MQ=255)
                      gTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTtc                                                               <  1:2453540/57‑1 (MQ=255)
                          tcaGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACt                                                  <  1:3109129/66‑1 (MQ=255)
                          tcaGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACt                                                  <  1:263335/66‑1 (MQ=255)
                           caGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCaaa                                                    >  1:1858144/1‑63 (MQ=255)
                           caGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCaaa                                                    >  1:2757383/1‑63 (MQ=255)
                           caGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCaaa                                                    >  1:3043470/1‑63 (MQ=255)
                            agcagcGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTtctc                                                             >  1:2806844/1‑53 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCaa                                                     >  1:2581852/1‑61 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCAt                                           >  1:1803866/1‑71 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCAt                                           >  1:879725/1‑71 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATAGTGTCTTTCTCAACTTCAAACTTGTTCAt                                           >  1:2727510/1‑71 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATAGTGTCTTTCTCAACTTCAAACTTGTTCAt                                           >  1:2732134/1‑71 (MQ=255)
                            aGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATAGTGTCTTTCTCAACTTCAAACTTGTTCAt                                           >  1:1075616/1‑71 (MQ=255)
                               aCCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTg                                        <  1:2061870/71‑1 (MQ=255)
                               aCCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTg                                        <  1:3474423/71‑1 (MQ=255)
                                ccGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTg                                        >  1:1401807/1‑70 (MQ=255)
                                 cGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTg                                        >  1:1184066/1‑69 (MQ=255)
                                 cGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATAGTGTCTTtctc                                                             <  1:3540968/48‑1 (MQ=255)
                                 cGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATAGTGTCTTtctc                                                             <  1:2894892/48‑1 (MQ=255)
                                  gCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAg                                     >  1:2120942/1‑71 (MQ=255)
                                  gCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAg                                     <  1:525508/71‑1 (MQ=255)
                                  gCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAg                                     >  1:2934137/1‑71 (MQ=255)
                                   cTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAgc                                    >  1:1171240/1‑71 (MQ=255)
                                     tGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAgcgc                                  >  1:3650985/1‑71 (MQ=255)
                                          gcgcgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTg                               <  1:2243418/69‑1 (MQ=255)
                                             cgGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGa                          <  1:2816432/71‑1 (MQ=255)
                                                cAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATaa                       <  1:3403249/71‑1 (MQ=255)
                                                 aGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAg                      >  1:953906/1‑71 (MQ=255)
                                                 aGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAg                      >  1:1719367/1‑71 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCaaa                                                    >  1:1334938/1‑38 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCaaa                                                    >  1:1180667/1‑38 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAgcg                                   >  1:2441881/1‑55 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTc                   >  1:2567131/1‑71 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTc                   >  1:2013196/1‑71 (MQ=255)
                                                    atgatgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTc                   >  1:3661496/1‑71 (MQ=255)
                                                       atgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTcgct                <  1:987120/71‑1 (MQ=255)
                                                       atgatgTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTcgct                >  1:3349803/1‑71 (MQ=255)
                                                             tCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGa                          >  1:1109941/1‑55 (MQ=255)
                                                             tCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTcgctggcgct          >  1:2447624/1‑71 (MQ=255)
                                                                  gATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCgctggct      >  1:2804010/1‑70 (MQ=255)
                                                                     tgtgTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAg  >  1:2757910/1‑71 (MQ=255)
                                                                     |                                                                      
CGCAATGTTGTTGCGATAGTAGGTCATCAGCACCGCTTGCTCGCGCGGCAGAATGATGATGTCGCCGATTGTGTCTTTCTCGACTTCAAACTTGTTCATTTGCAGCGCGTGATCGATAAGCTCGCTGGCGCTGGCTGAAG  >  W3110S.gb/4258217‑4258356

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: