Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,659,609 A→T 5.3% intergenic (+25/‑174) ynfD → / → ynfE hypothetical protein/oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,659,6090AT 5.3% 215.8 / 3.4 76intergenic (+25/‑174)ynfD/ynfEhypothetical protein/oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base T (2/2);  total (35/41)
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.31e-01

CCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATG  >  W3110S.gb/1659541‑1659677
                                                                    |                                                                    
ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGTatat                                                                     <  1:5490349/70‑1 (MQ=255)
ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                     >  1:2013496/1‑70 (MQ=255)
ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                     >  1:2113141/1‑70 (MQ=255)
ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                     <  1:2287893/70‑1 (MQ=255)
ccGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGata                                                                      <  1:3612878/69‑1 (MQ=255)
   cgcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGAc                                                                  >  1:1571254/1‑70 (MQ=255)
    gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                 >  1:4592458/1‑70 (MQ=255)
    gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                 >  1:260662/1‑70 (MQ=255)
      aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                              <  1:3397337/71‑1 (MQ=255)
      aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                              <  1:3116916/71‑1 (MQ=255)
        tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                            >  1:635609/1‑71 (MQ=255)
        tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                            >  1:4124584/1‑71 (MQ=255)
         cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                           >  1:6769075/1‑71 (MQ=255)
         cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                           >  1:2762958/1‑71 (MQ=255)
          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                          >  1:5752596/1‑71 (MQ=255)
          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                          >  1:6438097/1‑71 (MQ=255)
          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                          >  1:6447306/1‑71 (MQ=255)
          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                          >  1:3602005/1‑71 (MQ=255)
          tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                          >  1:2946915/1‑71 (MQ=255)
            gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                           <  1:2303064/68‑1 (MQ=38)
            gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                             <  1:297128/66‑1 (MQ=255)
            gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                             <  1:5692329/66‑1 (MQ=255)
            gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                           <  1:3496815/68‑1 (MQ=255)
             ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                 <  1:7320326/61‑1 (MQ=255)
              cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                           >  1:5551894/1‑66 (MQ=255)
               aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                     >  1:1312980/1‑71 (MQ=255)
               aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                     >  1:7750757/1‑71 (MQ=255)
                ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                             >  1:6250750/1‑62 (MQ=255)
                ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                    >  1:6264807/1‑71 (MQ=255)
                     gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGatat                                                                     <  1:663712/49‑1 (MQ=255)
                     gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                <  1:5361663/70‑1 (MQ=255)
                     gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                               >  1:6994427/1‑71 (MQ=255)
                      cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                               >  1:2898138/1‑70 (MQ=255)
                       gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTg                                                   <  1:2857769/65‑1 (MQ=255)
                        cgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                            <  1:367781/71‑1 (MQ=255)
                        cgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                            <  1:6228825/71‑1 (MQ=255)
                         gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                            >  1:3141447/1‑70 (MQ=255)
                            tGCCGATCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaac                                         <  1:4611776/70‑1 (MQ=255)
                            tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATg                                                                    <  1:7764482/43‑1 (MQ=255)
                            tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATg                                                                    <  1:2730468/43‑1 (MQ=255)
                            tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                            <  1:3853946/51‑1 (MQ=255)
                            tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                            <  1:2242572/51‑1 (MQ=255)
                            tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                        >  1:2764904/1‑71 (MQ=255)
                             gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                      >  1:5674882/1‑56 (MQ=255)
                              ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                        <  1:4965419/69‑1 (MQ=255)
                              ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                        <  1:2885520/69‑1 (MQ=255)
                                gAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                           >  1:7017804/1‑48 (MQ=25)
                                      cAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                  >  1:2981791/1‑51 (MQ=255)
                                        gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                            <  1:5381865/71‑1 (MQ=255)
                                        gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                            <  1:7088610/71‑1 (MQ=255)
                                        gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCt                             >  1:2740863/1‑70 (MQ=255)
                                         tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                    <  1:748844/46‑1 (MQ=255)
                                           acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                         <  1:4439356/71‑1 (MQ=255)
                                           acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                         <  1:3862338/71‑1 (MQ=255)
                                           acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                          <  1:4177389/70‑1 (MQ=255)
                                           acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                          <  1:7410301/70‑1 (MQ=255)
                                           acacATTATCGACTGAACGCCGGATAAGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                          >  1:3523984/1‑70 (MQ=255)
                                              caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                       <  1:4263754/54‑1 (MQ=255)
                                              caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                         <  1:7610873/68‑1 (MQ=255)
                                                ttATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATc                                   >  1:4969510/1‑56 (MQ=255)
                                                ttATCGACTGAACGCCGGATATGAAAAATCCGGCGATTTGAACATACAACATAATCCCACCttatta                        >  1:2404060/1‑67 (MQ=255)
                                                  aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                   <  1:7021548/70‑1 (MQ=255)
                                                   tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                  <  1:5709111/70‑1 (MQ=255)
                                                   tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                 >  1:3895719/1‑71 (MQ=255)
                                                       cTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt             >  1:4468933/1‑71 (MQ=255)
                                                          aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATaa                                     >  1:5493212/1‑44 (MQ=255)
                                                          aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc            <  1:5800203/69‑1 (MQ=255)
                                                          aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc            <  1:4140405/69‑1 (MQ=255)
                                                           aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa          <  1:4793747/70‑1 (MQ=255)
                                                           aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAt         <  1:2855376/71‑1 (MQ=255)
                                                            cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTATTCATACCCTTCTAtt        >  1:7247639/1‑71 (MQ=255)
                                                            cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt        <  1:4912921/71‑1 (MQ=255)
                                                              ccGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                      <  1:5161990/55‑1 (MQ=37)
                                                              ccGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa      <  1:3412591/71‑1 (MQ=255)
                                                               cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                         <  1:5670731/51‑1 (MQ=255)
                                                               cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                         <  1:600856/51‑1 (MQ=255)
                                                               cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat     <  1:5232603/71‑1 (MQ=255)
                                                               cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat     <  1:1922881/71‑1 (MQ=255)
                                                                ggATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata    >  1:7396876/1‑71 (MQ=255)
                                                                ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata    >  1:6745855/1‑71 (MQ=255)
                                                                ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTAatc    <  1:4418867/71‑2 (MQ=255)
                                                                 gATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                   >  1:6160229/1‑55 (MQ=255)
                                                                  atatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg  >  1:4827057/1‑71 (MQ=255)
                                                                   tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa          >  1:4234647/1‑62 (MQ=255)
                                                                   tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg  >  1:1667073/1‑70 (MQ=255)
                                                                    |                                                                    
CCGCGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATG  >  W3110S.gb/1659541‑1659677

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: