Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,328,254 C→T 14.3% D46D (GAC→GAT atoD → acetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,328,2540CT14.3% 208.2 / 22.8 77D46D (GAC→GATatoDacetyl‑CoA:acetoacetyl‑CoA transferase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base T (6/5);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GGCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGTCAATGGTCGAG  >  W3110S.gb/2328186‑2328324
                                                                    |                                                                      
ggCGGATTTATGGGGATTTGCACTCCATCCCGCCTGGTTGTAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:2799324/71‑1 (MQ=255)
ggCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:925609/71‑1 (MQ=255)
ggCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:1869305/71‑1 (MQ=255)
ggCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:3004970/71‑1 (MQ=255)
ggCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:1560010/71‑1 (MQ=255)
ggCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCt                                                                      <  1:3433367/71‑1 (MQ=255)
  cGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:1127012/71‑1 (MQ=255)
  cGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:2023007/71‑1 (MQ=255)
  cGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:3605170/71‑1 (MQ=255)
    gATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACa                                                                  <  1:1739121/71‑1 (MQ=255)
     aTTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACAt                                                                 <  1:2465052/71‑1 (MQ=255)
     aTTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACAt                                                                 <  1:3650258/71‑1 (MQ=255)
      tttATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACAt                                                                 <  1:970875/70‑1 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:3536713/1‑71 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:156290/1‑71 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:149841/1‑71 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:1434492/1‑71 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:960207/1‑71 (MQ=255)
        tATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGa                                                              >  1:2299405/1‑71 (MQ=255)
                  ggCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATa                                                            >  1:534859/1‑63 (MQ=255)
                      cTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCaa                                                       <  1:3112262/64‑1 (MQ=255)
                      cTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCaa                                                       <  1:2854200/64‑1 (MQ=255)
                      cTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCaa                                                       <  1:825555/64‑1 (MQ=255)
                      cTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCaa                                                       <  1:2549575/64‑1 (MQ=255)
                      cTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCaa                                                       <  1:2417793/64‑1 (MQ=255)
                        ccATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATg                                                     >  1:260031/1‑64 (MQ=255)
                        ccATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACcgc                                              >  1:2731920/1‑71 (MQ=255)
                         aaTCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACa                                                                  >  1:993056/2‑50 (MQ=255)
                           tCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGa                                                                    <  1:2241109/46‑1 (MQ=255)
                           tCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:2134710/46‑1 (MQ=255)
                           tCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:2660046/46‑1 (MQ=255)
                           tCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGa                                                                    <  1:1235251/46‑1 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:1634444/1‑46 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:1588714/1‑46 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:32698/1‑46 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:66267/1‑46 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:1975453/1‑46 (MQ=255)
                             ccGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTttg                                         >  1:3062063/1‑71 (MQ=255)
                             acGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACa                                                                  >  1:126465/2‑46 (MQ=255)
                               gCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTttgtt                                       <  1:836581/71‑1 (MQ=255)
                               gCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTttgtt                                       <  1:2557182/71‑1 (MQ=255)
                                   ggTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATa                                   >  1:2258550/1‑71 (MQ=255)
                                   ggTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATa                                   >  1:2876487/1‑71 (MQ=255)
                                    gTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACATTGATAGCCAATGATACCGCGttt                                          <  1:2562103/63‑1 (MQ=255)
                                    gTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACATTGATAGCCAATGATACCGCGttt                                          <  1:1238941/63‑1 (MQ=255)
                                    gTTGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACATTGATAGCCAATGATACCGCGttt                                          <  1:2392838/63‑1 (MQ=255)
                                     ttGAAGCATTACTGGAATCTGGTGTTCGCGATCTGACATTGATAGCCAATGATACCGCGttt                                          <  1:1234368/62‑1 (MQ=255)
                                      tGAAGCATTAGTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATa                                                  >  1:41898/1‑53 (MQ=255)
                                      tGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATa                                                  >  1:3099320/1‑53 (MQ=255)
                                      tGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATa                                                  >  1:157188/1‑53 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCgg                               <  1:2440920/70‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCgg                               <  1:544067/70‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCgg                               <  1:3501240/70‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGc                              <  1:1336501/71‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGc                              <  1:1991678/71‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGc                              <  1:3306895/71‑1 (MQ=255)
                                        aaGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGACAATGATACCGCGTTTGTTGATACCGGc                              <  1:728255/71‑1 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGa                                                    >  1:2916479/1‑42 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGa                                                    >  1:67641/1‑42 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGa                                                    >  1:3156494/1‑42 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGa                                                    >  1:2561674/1‑42 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGa                                                    >  1:203671/1‑42 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTc                       >  1:2862500/1‑71 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTc                       >  1:1718201/1‑71 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTc                       >  1:1640386/1‑71 (MQ=255)
                                               aCTGGAATCTGGTGTTCGCGACCTGACATAGATAGCCAATGa                                                    >  1:1049205/1‑42 (MQ=255)
                                                 tGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAAt                                                      <  1:3149179/38‑1 (MQ=255)
                                                 tGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAAt                                                      <  1:3018868/38‑1 (MQ=255)
                                                            gtTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGa                                     >  1:3712144/1‑44 (MQ=255)
                                                            gtTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGa                                     >  1:153215/1‑44 (MQ=255)
                                                            gtTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGa                                     >  1:3129387/1‑44 (MQ=255)
                                                            gtTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGa                                     >  1:2356656/1‑44 (MQ=255)
                                                             ttCGCGACCTGACATTGATAGCCAATGATACCGCGTttgtt                                       <  1:2690048/41‑1 (MQ=255)
                                                             ttCGCGACCTGACATTGATAGCCAATGATACCGCGTttgtt                                       <  1:1126776/41‑1 (MQ=255)
                                                             ttCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGTCAAt         >  1:1621342/1‑71 (MQ=255)
                                                             ttCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGACAAt         >  1:3387529/1‑71 (MQ=255)
                                                               cgcgACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCg              >  1:2288437/1‑64 (MQ=255)
                                                                    ccTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGTCAATGGTCGAg  <  1:382341/71‑1 (MQ=255)
                                                                    ccTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGTCAATGGTCGAg  <  1:2992232/71‑1 (MQ=255)
                                                                    |                                                                      
GGCGGATTTATGGGGATTGGCACTCCATCCCGCCTGGTTGAAGCATTACTGGAATCTGGTGTTCGCGACCTGACATTGATAGCCAATGATACCGCGTTTGTTGATACCGGCATCGGTCCGCTCATCGTCAATGGTCGAG  >  W3110S.gb/2328186‑2328324

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: