Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,244,991 T→A 7.7% I286F (ATC→TTC)  xylE ← D‑xylose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,244,9910TA 7.7% 220.9 / 9.3 78I286F (ATC→TTC) xylED‑xylose transporter
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base A (4/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

GCCCCCAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTGCGGC  >  W3110S.gb/4244921‑4245059
                                                                      |                                                                    
gCCCCCAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAt                                                                      >  1:860852/1‑71 (MQ=255)
    ccAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATg                                                                     >  1:11382/1‑68 (MQ=255)
    ccAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATg                                                                     >  1:2761984/1‑68 (MQ=255)
     cAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAAGga                                                                   >  1:4520762/1‑69 (MQ=255)
     cAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAAGga                                                                   >  1:2205996/1‑69 (MQ=255)
         gTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGc                                                               >  1:4920715/1‑69 (MQ=255)
         gTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCAt                                                             >  1:1305237/1‑71 (MQ=255)
         gTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCAt                                                             <  1:4149465/71‑1 (MQ=255)
         gTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCAt                                                             >  1:2740475/1‑71 (MQ=255)
          ttttGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCAt                                                             >  1:3942499/1‑70 (MQ=255)
            ttGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTAc                                                          >  1:662898/1‑71 (MQ=255)
            ttGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTAc                                                          >  1:2473787/1‑71 (MQ=255)
            ttGAACACATCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTa                                                           >  1:2059808/1‑70 (MQ=255)
             tGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTAc                                                          <  1:385501/70‑1 (MQ=255)
              gAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGc                                                        >  1:2622052/1‑71 (MQ=255)
                acacTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCg                                                      <  1:1755395/71‑1 (MQ=255)
                  acTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAAGGagag                                                                <  1:2603462/59‑1 (MQ=255)
                      ccGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTAc                                                 >  1:1423425/1‑70 (MQ=255)
                      ccGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTAc                                                 >  1:2256063/1‑70 (MQ=255)
                      ccGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAAGGAGAGCATTACGCCGATTAc                                                 >  1:4097993/1‑70 (MQ=255)
                       cGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACaa                                               <  1:4364443/71‑1 (MQ=255)
                       cGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACaa                                               <  1:1398645/71‑1 (MQ=255)
                        ggCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAAt                                              <  1:4728989/71‑1 (MQ=255)
                        ggCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGAAGGAGAGCATTACGCCGATTACAAt                                              <  1:1976760/71‑1 (MQ=255)
                                gtaCAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATcacgccca                                      >  1:617795/1‑71 (MQ=255)
                                 taCAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCAt                                                             <  1:691296/47‑1 (MQ=255)
                                 taCAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATcacgcc                                        <  1:4024253/68‑1 (MQ=255)
                                   cAGCACCACATTGATGCCGACAAATAGCTGGAAGATGGAGAGCATTACGCCGATTACAAAcacgcccacgc                                   >  1:4994069/1‑71 (MQ=255)
                                     gCACCACATTGATGCCGACAAATTGCTGGAAGATGGagag                                                                <  1:149450/40‑1 (MQ=255)
                                      caccacATTGATGCCGACAAATTGCTGGAAGATGGagag                                                                <  1:2140783/39‑1 (MQ=255)
                                       accacATTGATGCCGACAAATTGCTGGAAGATGGAGAGCa                                                              >  1:1814479/1‑40 (MQ=255)
                                       accacATTGATGCCGACAAATTGCTGGAAGATGGAGAGCa                                                              >  1:1446395/1‑40 (MQ=255)
                                       accacATTGATGCCGACAAATTGCTGGAAGATGGAGAGCa                                                              >  1:3400585/1‑40 (MQ=255)
                                       accacATTGATGCCGACAAATTGCTGGAAGATGGAGAGCa                                                              >  1:3601288/1‑40 (MQ=255)
                                         cacaTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCg                                                      <  1:4549450/46‑1 (MQ=255)
                                         cacaTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCg                                                      <  1:4243496/46‑1 (MQ=255)
                                         cacaTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCg                                                      <  1:4025107/46‑1 (MQ=255)
                                         cacaTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCg                                                      <  1:4054026/46‑1 (MQ=255)
                                          acaTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACAt                            >  1:5094649/1‑71 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACAt                            <  1:358133/70‑1 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACAt                            <  1:3522620/70‑1 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATc                           <  1:3992115/71‑1 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATc                           <  1:4454913/71‑1 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATc                           <  1:117486/71‑1 (MQ=255)
                                           caTTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATc                           <  1:4087905/71‑1 (MQ=255)
                                             ttGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcag                         >  1:2498297/1‑71 (MQ=255)
                                             ttGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATc                           <  1:3528243/69‑1 (MQ=255)
                                              tGATGCCGACAAATTGCTGGAAGATGGAGAGCATTa                                                           >  1:2681606/1‑36 (MQ=255)
                                              tGATGCCGACAAATTGCTGGAAGATGGAGAGCATTa                                                           >  1:1729055/1‑36 (MQ=255)
                                              tGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcag                         >  1:799049/1‑70 (MQ=255)
                                              tGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcag                         >  1:452947/1‑70 (MQ=255)
                                               gATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcag                         >  1:158785/1‑69 (MQ=255)
                                               gATGCCGACAAATTGCTGGAAGAAGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcag                         >  1:4548559/1‑69 (MQ=255)
                                                aTGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCaaa                               >  1:3772665/1‑62 (MQ=255)
                                                aTGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCaaa                               >  1:3943663/1‑62 (MQ=255)
                                                aTGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagca                       <  1:4135222/70‑1 (MQ=255)
                                                aTGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagca                       <  1:843780/70‑1 (MQ=255)
                                                aTGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagca                       <  1:4554311/70‑1 (MQ=255)
                                                 tGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATcacgccca                                      <  1:4380180/54‑1 (MQ=255)
                                                  gCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAgac                    >  1:1898915/1‑71 (MQ=255)
                                                   ccGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAgac                    <  1:3905208/70‑1 (MQ=255)
                                                    cGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAgacga                  >  1:457011/1‑71 (MQ=255)
                                                     gACAAATTGCTGGAAGATGGAGAGCCTTACGCCGATTACa                                                <  1:4609779/40‑1 (MQ=255)
                                                     gACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACa                                                <  1:1020782/40‑1 (MQ=255)
                                                      aCAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACa                                                <  1:1228528/39‑1 (MQ=255)
                                                      aCAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagc                        <  1:927527/63‑1 (MQ=255)
                                                      aCAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagc                        <  1:1621598/63‑1 (MQ=255)
                                                       cAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGCCGacca               >  1:2730388/1‑71 (MQ=255)
                                                       cAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGacca               >  1:571826/1‑71 (MQ=255)
                                                        aaaTTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAAc                              <  1:839097/55‑1 (MQ=255)
                                                        aaaTTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAga                     >  1:224721/1‑64 (MQ=255)
                                                              cTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGttt        <  1:2384681/71‑1 (MQ=255)
                                                                  aaGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGCCCACCGGTTTTg      >  1:628661/1‑69 (MQ=255)
                                                                  aaGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTg      >  1:3384548/1‑69 (MQ=255)
                                                                  aaGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTg      >  1:4145738/1‑69 (MQ=255)
                                                                   aGATGGAGAGCATTACGCGGATTACAATCACGCCCACGCGAAAc                              <  1:4087060/44‑1 (MQ=38)
                                                                    gATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTGCGGc  <  1:220129/71‑1 (MQ=255)
                                                                    gATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTGCGGc  <  1:3493212/71‑1 (MQ=255)
                                                                     aTGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATcagcag                      <  1:396909/50‑1 (MQ=255)
                                                                      |                                                                    
GCCCCCAGCGTTTTGAACACTTCCGGCGCGTAGTACAGCACCACATTGATGCCGACAAATTGCTGGAAGATGGAGAGCATTACGCCGATTACAATCACGCCCACGCCAAACATCAGCAGACGACCACCGGTTTTGCGGC  >  W3110S.gb/4244921‑4245059

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: ‑