Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,135,166 A→G 68.8% intergenic (‑139/+153) pitB ← / ← gss phosphate transporter/fused glutathionylspermidine amidase and glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,135,1660AG68.8% 89.1 / 52.5 64intergenic (‑139/+153)pitB/gssphosphate transporter/fused glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (15/5);  new base G (21/23);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 5.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATA  >  W3110S.gb/3135098‑3135230
                                                                    |                                                                
gAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCaaa                                                                  <  1:1383872/69‑1 (MQ=255)
gAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCaaa                                                                  >  1:3465238/1‑69 (MQ=255)
gAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCaaa                                                                  >  1:42728/1‑69 (MQ=255)
gAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCaaa                                                                  <  1:489854/69‑1 (MQ=255)
  aaTATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGAc                                                                >  1:3119807/1‑69 (MQ=255)
      tCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGAcc                                                               >  1:3325539/1‑66 (MQ=255)
      tCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAAcc                                                               >  1:1662590/1‑66 (MQ=255)
      tCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAAcc                                                               >  1:1757957/1‑66 (MQ=255)
        aTTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTa                                                           >  1:914762/1‑68 (MQ=255)
          tttttAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGa                                                         <  1:3456865/68‑1 (MQ=255)
          tttttAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGa                                                         <  1:378573/68‑1 (MQ=255)
             ttAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAaat                                                      <  1:1984177/68‑1 (MQ=255)
                 cGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATc                                                  <  1:2691042/68‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATc                                                  <  1:3115874/67‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaa                                                <  1:2082819/69‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaa                                                <  1:2157246/69‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaa                                                <  1:1318150/69‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaa                                                <  1:81813/69‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCa                                                 <  1:176161/68‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCa                                                 <  1:2391731/68‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCaa                                                <  1:2751882/69‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCa                                                 <  1:550930/68‑1 (MQ=255)
                  gATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCa                                                 <  1:417768/68‑1 (MQ=255)
                   aTAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaaa                                               <  1:1945226/69‑1 (MQ=255)
                    tAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaaaa                                              <  1:569598/69‑1 (MQ=255)
                     aaaaaGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCaaaa                                              >  1:1929705/1‑68 (MQ=255)
                     aaaaaGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCaaaa                                              >  1:677367/1‑68 (MQ=255)
                           tttgtTATTTTTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGt                                       >  1:3489959/3‑69 (MQ=255)
                               ttATTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   <  1:3189559/69‑1 (MQ=255)
                                 aTTTGTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAgg                                        >  1:17262/1‑62 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCaaa                                                                  <  1:710389/36‑1 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAgg                                        >  1:1018239/1‑62 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTCCa                                 >  1:863888/1‑69 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAgg                                        >  1:2466477/1‑62 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAgg                                        >  1:2454787/1‑62 (MQ=255)
                                 aTTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAgg                                        >  1:1988065/1‑62 (MQ=255)
                                  tttCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:2461815/1‑66 (MQ=255)
                                  tttCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:2623572/1‑66 (MQ=255)
                                  tttCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:3158795/1‑66 (MQ=255)
                                   ttCTTGTGTATCAGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   <  1:3660865/65‑1 (MQ=255)
                                              cAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAgg                                        >  1:2569256/1‑49 (MQ=255)
                                               aGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:1721078/1‑53 (MQ=255)
                                               aGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:3708262/1‑53 (MQ=255)
                                               aGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:2220787/1‑53 (MQ=255)
                                               aGACAATAAGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTc                                   >  1:188384/1‑53 (MQ=255)
                                               aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCa                                                 >  1:826848/1‑39 (MQ=255)
                                               aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCa                                                 >  1:97531/1‑39 (MQ=255)
                                               aGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCa                                                 >  1:2132067/1‑39 (MQ=255)
                                                   aaTAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTcaca               >  1:541399/1‑69 (MQ=255)
                                                      aaGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGt                                       <  1:1320343/42‑1 (MQ=255)
                                                      aaGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGt                                       >  1:1902896/1‑42 (MQ=255)
                                                      aaGCATTCATTCAAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTaa                         <  1:1177423/56‑1 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAgg                                        <  1:1583480/35‑1 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAgg                                        <  1:1602181/35‑1 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAgg                                        <  1:725472/35‑1 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAgg                                        <  1:626386/35‑1 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATatgaaatg      >  1:2007942/1‑69 (MQ=255)
                                                            tcattcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATatgaaatg      >  1:3242729/1‑69 (MQ=255)
                                                             cattcaAGACCTTTAGAAATAATCAAAACACAGGTAt                                     <  1:2787402/37‑1 (MQ=255)
                                                             cattcaAGACCTTTAGAAATAATCAAAACACAGGTAt                                     <  1:1331321/37‑1 (MQ=255)
                                                             cattcaAGACCTTTAGAAATAATCAAAACACAGGTAt                                     <  1:2321836/37‑1 (MQ=255)
                                                             cattcaAGACCTTTAGAAATAATCAAAACACAGGTATTcc                                  >  1:1991506/1‑40 (MQ=255)
                                                             cattcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATatgaaatg      >  1:1772349/1‑68 (MQ=255)
                                                                tcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGaata  >  1:1982090/1‑69 (MQ=255)
                                                                tcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGaata  >  1:2600584/1‑69 (MQ=255)
                                                                tcaAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGaata  >  1:3067467/1‑69 (MQ=255)
                                                                 caAGACCTTTAGAAATAATCAAAACACAGGTATTCCACt                               <  1:2864152/39‑1 (MQ=255)
                                                                 caAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTcaca               >  1:2229433/1‑55 (MQ=255)
                                                                 caAGACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTcaca               >  1:2233667/1‑55 (MQ=255)
                                                                    |                                                                
GAAATATCATTTTTTAACGATAAAAAGAGTGTTATTTCTTGTGTATCAGACAATAAGCATTCATTCAAAACCTTTAGAAATAATCAAAACACAGGTATTCCACTGGTTAATATATTCACATATGAAATGAATA  >  W3110S.gb/3135098‑3135230

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: