Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,125,813 C→T 22.4% P152L (CCT→CTT)  rpe → D‑ribulose‑5‑phosphate 3‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,125,8130CT22.4% 136.5 / 35.3 67P152L (CCT→CTT) rpeD‑ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (28/24);  new base T (8/7);  total (36/31)
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.25e-01

TGATGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAGA  >  W3110S.gb/4125744‑4125882
                                                                     |                                                                     
tGATGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCt                                                                      >  1:2164892/1‑71 (MQ=255)
  aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTc                                                                     <  1:2047813/70‑1 (MQ=255)
  aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCa                                                                    <  1:469203/71‑1 (MQ=255)
  aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCa                                                                    <  1:4588347/71‑1 (MQ=255)
   tGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCa                                                                    >  1:1740291/1‑70 (MQ=255)
    ggATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCaaa                                                                  >  1:3455987/1‑71 (MQ=255)
     gATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCaaa                                                                  <  1:4439088/70‑1 (MQ=255)
      aTAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAaca                                                                >  1:2565559/1‑71 (MQ=255)
      aTAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAaca                                                                <  1:3496162/71‑1 (MQ=255)
            tGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGa                                                           >  1:3448811/1‑70 (MQ=255)
             ggATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGAt                                                          >  1:945471/1‑70 (MQ=255)
                    aTCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATa                                                         <  1:736268/64‑1 (MQ=255)
                    aTCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTgcg                                                  >  1:1407248/1‑71 (MQ=255)
                      ccTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTgcgc                                                 >  1:124484/1‑70 (MQ=255)
                       ctgctgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGcgcg                                                <  1:3853965/70‑1 (MQ=255)
                         gctgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                             >  1:3288731/1‑71 (MQ=255)
                             aTGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTAc                                          <  1:825688/70‑1 (MQ=255)
                              tGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCaaa                                                                  <  1:3837104/45‑1 (MQ=255)
                                 ccGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCa                                                                    <  1:1317778/40‑1 (MQ=255)
                                   gTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCaaa                                                                  >  1:1159043/1‑40 (MQ=255)
                                    tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAAc                                                      <  1:264121/51‑1 (MQ=255)
                                    tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAAc                                                      <  1:1471552/51‑1 (MQ=255)
                                    tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                             >  1:3625114/1‑60 (MQ=255)
                                    tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTACTCAAACACTGGATAAACTGcgcg                                                <  1:1716061/57‑1 (MQ=255)
                                     cAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTa                                 >  1:1621784/1‑71 (MQ=255)
                                      aaCCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTAt                                >  1:2203428/1‑71 (MQ=255)
                                      aaCCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTa                                 >  1:3305553/1‑70 (MQ=255)
                                        cccTGGTTTCGGCTGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGcgcg                                                <  1:999435/53‑1 (MQ=255)
                                         ccTGGTTTCGTCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcga                             >  1:344909/1‑71 (MQ=255)
                                          cTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                            <  1:3560148/71‑1 (MQ=255)
                                           tGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTa                                           >  1:2901659/1‑55 (MQ=255)
                                           tGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                            >  1:3051251/1‑70 (MQ=255)
                                           tGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                            >  1:469905/1‑70 (MQ=255)
                                            ggTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                            >  1:3565810/1‑69 (MQ=255)
                                             gTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTg                                                    >  1:725884/1‑44 (MQ=255)
                                               ttCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTc                       >  1:3516190/1‑71 (MQ=255)
                                               ttCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTc                       >  1:1877696/1‑71 (MQ=255)
                                                tCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAAc                                                      >  1:2624732/1‑39 (MQ=255)
                                                tCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCt                      >  1:898616/1‑71 (MQ=255)
                                                   gcggTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                             <  1:3442118/45‑1 (MQ=255)
                                                      gTCAGTCTTTCATTCCTCAAACACTGGATAAACTGcgcg                                                <  1:2972847/39‑1 (MQ=255)
                                                      gTCAGTCTTTCATTCCTCAAACACTGGATAAACTGcgcg                                                <  1:2707590/39‑1 (MQ=255)
                                                      gTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcga                             >  1:1174775/1‑58 (MQ=255)
                                                      gTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCttt                <  1:2856406/71‑1 (MQ=255)
                                                      gTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCttt                <  1:3908355/71‑1 (MQ=255)
                                                       tCAGTCTTTCATTCCTCAAACACTGGATAAACTGcgcg                                                <  1:519976/38‑1 (MQ=255)
                                                       tCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTg               <  1:3971204/71‑1 (MQ=255)
                                                        cAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGt                                        >  1:2990005/1‑45 (MQ=255)
                                                        cAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa              >  1:768063/1‑71 (MQ=255)
                                                         aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa              >  1:693349/1‑70 (MQ=255)
                                                         aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGAc             <  1:345005/71‑1 (MQ=255)
                                                         aGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTAc                                          <  1:1648539/42‑1 (MQ=255)
                                                         aGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa              >  1:2416528/1‑70 (MQ=255)
                                                         aGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa              <  1:4350823/70‑1 (MQ=255)
                                                          gTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                             <  1:4060074/38‑1 (MQ=255)
                                                          gTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACa            <  1:2769719/71‑1 (MQ=255)
                                                             tttCATTCTTCAAACACTGGATAAACTGCGCGAAGTAc                                          <  1:3422406/38‑1 (MQ=255)
                                                             tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcga                             >  1:3544888/1‑51 (MQ=255)
                                                             tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc         >  1:2416914/1‑71 (MQ=255)
                                                             tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc         >  1:1141756/1‑71 (MQ=255)
                                                             tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc         >  1:1555702/1‑71 (MQ=255)
                                                               tCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTAt                                >  1:4629609/1‑46 (MQ=255)
                                                               tCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGa       >  1:3357877/1‑71 (MQ=255)
                                                               tCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGa       >  1:83728/1‑71 (MQ=255)
                                                                 aTTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGAc      <  1:2213906/70‑1 (MQ=255)
                                                                  ttcttcAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGAc      >  1:1760407/1‑69 (MQ=255)
                                                                  ttCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgacga                          <  1:544886/49‑1 (MQ=255)
                                                                  ttCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAg                         <  1:914307/50‑1 (MQ=255)
                                                                    ccTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTATGAc             >  1:4157636/1‑60 (MQ=255)
                                                                     cTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATc                               <  1:424816/41‑1 (MQ=255)
                                                                     cTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc         <  1:3642581/63‑1 (MQ=255)
                                                                     cTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAGa  <  1:1329069/70‑1 (MQ=255)
                                                                     |                                                                     
TGATGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAGA  >  W3110S.gb/4125744‑4125882

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: