Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,125,813 C→T 35.3% 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,8130CT35.3% 74.0 / 65.5 68P152L (CCT→CTT) rpeD‑ribulose‑5‑phosphate 3‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (26/18);  new base T (13/11);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAG  >  W3110S.gb/4125746‑4125881
                                                                   |                                                                    
aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCAttcttc                                                                    >  1:1006326/1‑70 (MQ=255)
aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCa                                                                   <  1:150102/71‑1 (MQ=255)
aTGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTc                                                                    >  1:4040971/1‑70 (MQ=255)
   gATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCaaa                                                                 <  1:1059307/70‑1 (MQ=255)
   gATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCaaa                                                                 >  1:1671282/1‑70 (MQ=255)
    aTAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAaca                                                               <  1:3141587/71‑1 (MQ=255)
     tAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAaca                                                               <  1:3388535/70‑1 (MQ=255)
      aaGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACt                                                             >  1:414135/1‑71 (MQ=255)
            gATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATa                                                        >  1:2248215/1‑70 (MQ=255)
            gATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAaca                                                               >  1:2910649/1‑63 (MQ=255)
            gATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAaca                                                               >  1:3936247/1‑63 (MQ=255)
            gATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATaa                                                       <  1:3853747/71‑1 (MQ=255)
            gATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:3709870/70‑1 (MQ=255)
                tgATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATaaa                                                      >  1:241813/1‑68 (MQ=255)
                tgATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTg                                                   >  1:2347721/1‑71 (MQ=255)
                     ctgctgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAacac                                                              <  1:492376/55‑1 (MQ=255)
                     ctgctgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGcgcg                                               <  1:2043203/70‑1 (MQ=255)
                        ctgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            <  1:3227241/70‑1 (MQ=255)
                         tgATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTa                                          >  1:822912/1‑71 (MQ=255)
                             gTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTc                                      >  1:217383/1‑71 (MQ=255)
                             gTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCaaa                                                                 >  1:3187313/1‑44 (MQ=255)
                             gTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCaaa                                                                 >  1:651708/1‑44 (MQ=255)
                                cGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGcc                                   >  1:623965/1‑71 (MQ=255)
                                cGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGcc                                   >  1:2965065/1‑71 (MQ=255)
                                 gTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:2153402/49‑1 (MQ=255)
                                 gTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:1877801/49‑1 (MQ=255)
                                 gTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:1474222/49‑1 (MQ=255)
                                  tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGt                                 <  1:3093418/71‑1 (MQ=255)
                                  tCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGt                                 <  1:11781/71‑1 (MQ=255)
                                    aaCCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTAt                               >  1:746538/1‑71 (MQ=255)
                                    aaCCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTAt                               >  1:2591666/1‑71 (MQ=255)
                                      cccTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcg                             <  1:2807194/71‑1 (MQ=255)
                                        cTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAaca                                                               >  1:497168/1‑35 (MQ=255)
                                        cTGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACt                                                             >  1:1593583/1‑37 (MQ=255)
                                        cTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                           <  1:1093850/71‑1 (MQ=255)
                                        cTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                           <  1:1559391/71‑1 (MQ=255)
                                        cTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                           <  1:3987414/71‑1 (MQ=255)
                                         tGGTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                           >  1:3666950/1‑70 (MQ=255)
                                         tGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:1970680/41‑1 (MQ=255)
                                         tGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACAATGGATa                                                        <  1:527325/41‑1 (MQ=255)
                                          ggTTTCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgac                           >  1:3196070/1‑69 (MQ=255)
                                            tttCGGCGGTCAGTCTTTCATTCTTCAAACACTGGATAAACTGcgcg                                               <  1:3230591/47‑1 (MQ=255)
                                            tttCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATcgacga                         <  1:3008927/69‑1 (MQ=255)
                                             ttCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTc                      >  1:511992/1‑71 (MQ=255)
                                              tCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATa                                                        <  1:969552/36‑1 (MQ=255)
                                                   ggTCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCtt                <  1:3467618/71‑1 (MQ=255)
                                                     tCAGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTg              <  1:514640/71‑1 (MQ=255)
                                                       aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa             >  1:3186658/1‑70 (MQ=255)
                                                       aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa             <  1:299288/70‑1 (MQ=255)
                                                       aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa             >  1:1593375/1‑70 (MQ=255)
                                                       aGTCTTTCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGa             >  1:692075/1‑70 (MQ=255)
                                                        gTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            <  1:1431271/38‑1 (MQ=255)
                                                        gTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            <  1:2642473/38‑1 (MQ=255)
                                                         tCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACAt          <  1:1822350/71‑1 (MQ=255)
                                                          cTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            >  1:2970608/1‑36 (MQ=255)
                                                          cTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            >  1:1868758/1‑36 (MQ=255)
                                                          cTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            >  1:1745166/1‑36 (MQ=255)
                                                          cTTTCATTCCTCAAACACTGGATAAACTGCGCGAAg                                            >  1:1595194/1‑36 (MQ=255)
                                                          cTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACAtt         >  1:1734715/1‑71 (MQ=255)
                                                           tttCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACAtt         >  1:262493/1‑70 (MQ=255)
                                                           tttCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc        >  1:524103/1‑71 (MQ=255)
                                                           tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTg                    >  1:1096742/1‑59 (MQ=255)
                                                           tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTg                    >  1:1930719/1‑59 (MQ=255)
                                                           tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTg                    >  1:1432677/1‑59 (MQ=255)
                                                           tttCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTc        >  1:161220/1‑71 (MQ=255)
                                                             tCATTCTTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACa           >  1:3459135/1‑66 (MQ=255)
                                                             tCATTCCTCAAACACTGGATAAACTGCGCGAAGTa                                          >  1:412353/1‑35 (MQ=255)
                                                             tCATTCCTCAAACACTGGATAAACTGCGCGAAGTa                                          >  1:975183/1‑35 (MQ=255)
                                                             tCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGa      >  1:1842880/1‑71 (MQ=255)
                                                                ttcttcAAACACTGGATAAACTGCGCGAAGTACGTc                                      <  1:62488/36‑1 (MQ=255)
                                                                ttcttcAAACACTGGATAAACTGCGCGAAGTACGTc                                      <  1:813500/36‑1 (MQ=255)
                                                                  ccTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAg  <  1:1210677/70‑1 (MQ=255)
                                                                   |                                                                    
ATGGATAAGCTGGATGTGATCCTGCTGATGTCCGTCAACCCTGGTTTCGGCGGTCAGTCTTTCATTCCTCAAACACTGGATAAACTGCGCGAAGTACGTCGCCGTATCGACGAGTCTGGCTTTGACATTCGACTAG  >  W3110S.gb/4125746‑4125881

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: