Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,298,651 Δ1 bp coding (432/1359 nt) glpT ← sn‑glycerol‑3‑phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,298,6510A.93.5% 129.2 / 4.0 31coding (432/1359 nt)glpTsn‑glycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (12/17);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

CCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACC  >  REL606/2298580‑2298720
                                                                       |                                                                     
ccGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACg                                                                  <  1:437285/76‑1 (MQ=255)
    cGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACca                                                              <  1:1276543/76‑1 (MQ=255)
     gTTGTGCGCATAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACcac                                                             <  2:1270942/76‑1 (MQ=255)
      ttGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACcac                                                             <  1:1053917/76‑1 (MQ=255)
      ttGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACcaca                                                            <  1:638408/76‑1 (MQ=255)
       tgtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCCTGTACGACCacac                                                           <  1:1222795/76‑1 (MQ=255)
        gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACAcg                                                          >  1:1529418/1‑76 (MQ=255)
         tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACAcgg                                                         >  1:804057/1‑76 (MQ=255)
         ggcgcACCGCTCCCCGCCGACACAACGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACTcga                                                         <  2:1478723/75‑2 (MQ=255)
          gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACAcggc                                                        <  1:394393/76‑1 (MQ=255)
            gcACACTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACAcggcgg                                                      <  2:1813775/76‑1 (MQ=255)
                aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCAc                                                  <  2:1321849/76‑1 (MQ=255)
                  ttCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCAccc                                                >  2:1472258/1‑76 (MQ=255)
                   tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCAcccc                                               <  1:615444/76‑1 (MQ=255)
                    ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCa                                              >  2:722378/1‑76 (MQ=255)
                       cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATcc                                           >  2:513345/1‑76 (MQ=255)
                           gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTg                                       <  2:470824/76‑1 (MQ=255)
                               caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGgac                                   >  1:1705021/1‑76 (MQ=255)
                                 aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGaacca                                 <  1:526081/76‑1 (MQ=255)
                                  tGCCGCCACGTTCTNTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGaaccaa                                >  1:994701/1‑76 (MQ=255)
                                   gccgccACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCCGGaacccac                               >  1:1152992/1‑76 (MQ=255)
                                    ccgccACGTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGaaccaacc                              <  2:384412/76‑1 (MQ=255)
                                           gTTCTTTCTGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAgg                       <  2:674481/76‑1 (MQ=255)
                                                ttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaaca                   >  1:1117796/1‑76 (MQ=255)
                                                ttctGCGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaacaa                  <  1:200123/76‑1 (MQ=255)
                                                tgctGCGACCACCAGTGCGCCATGTCCGACCACACGGCGGCCTCCCCATCCCCTGGAACCAACCGCAGAGGaacaa                  <  2:196995/74‑1 (MQ=255)
                                                    gcgcCCACCAGTGCACCATGTACGACCACACGGCGTCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGt              <  1:512941/76‑1 (MQ=255)
                                                     cGACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCCTCCCCTGGAACCAACCGCAGAGGAACAACAgaa             >  2:348532/1‑74 (MQ=255)
                                                      gACCACCAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTAc            >  1:1308743/1‑76 (MQ=255)
                                                       accaccAGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACa           <  2:349434/76‑1 (MQ=255)
                                                            caGTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt      <  1:347483/76‑1 (MQ=255)
                                                              gTGCATCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa    <  2:1817038/76‑1 (MQ=255)
                                                              gTGCACCATGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa    >  1:1400762/1‑76 (MQ=255)
                                                               tGCACCATGTACGNCCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCAc   >  1:1234184/1‑76 (MQ=255)
                                                                gCACCATGTACGGCCACACGGCGCCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCAcc  <  2:921868/76‑1 (MQ=255)
                                                                       |                                                                     
CCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACC  >  REL606/2298580‑2298720

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: