Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,233,034 C→A 7.6% intergenic (+80/+11) yjbB → / ← pepE predicted transporter/(alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,233,0340CA 7.6% 255.3 / 4.3 78intergenic (+80/+11)yjbB/pepEpredicted transporter/(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base C (44/28);  new base A (2/4);  total (46/32)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAG  >  W3110S.gb/4232969‑4233103
                                                                  |                                                                     
tAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCaa                                                                  >  1:5579175/1‑71 (MQ=255)
  ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                >  1:3184792/1‑71 (MQ=255)
  ggCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                >  1:2920410/1‑71 (MQ=255)
   gCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTc                                                               <  1:1060183/71‑1 (MQ=255)
    ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCg                                                                      >  1:3779133/1‑63 (MQ=255)
    ccTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAAtt                                                                <  1:1974841/69‑1 (MQ=255)
      tGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                            >  1:616697/1‑71 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTc                                                               <  1:3490948/66‑1 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACtt                                                           >  1:1363010/1‑70 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                            <  1:3125710/69‑1 (MQ=255)
        aTAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACt                                                            >  1:7645225/1‑69 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaa                                                       <  1:228958/71‑1 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaa                                                       <  1:7801562/71‑1 (MQ=255)
           aGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaa                                                        <  1:5953320/70‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:5485127/71‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:1735210/71‑1 (MQ=255)
              cgcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:7239894/71‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                   <  1:7575328/71‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAgg                                                                     <  1:3309026/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAgg                                                                     <  1:4842987/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAgg                                                                     <  1:1986592/53‑1 (MQ=255)
               gcgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAgg                                                                     <  1:7314295/53‑1 (MQ=255)
                cgTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGATGCATG‑ATCCgg                                                                     <  1:248578/52‑1 (MQ=18)
                 gTCAAGCGTCGCATCAGGGACTATGAACAAGATGCCGGAT‑CATGCATCAGGCa                                                                   >  1:5545544/1‑53 (MQ=37)
                 gTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCAGGCa                                                                   >  1:7273669/1‑53 (MQ=255)
                  tCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATcc                                                                       <  1:5927415/48‑1 (MQ=255)
                   cAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGAc                                               >  1:5121009/1‑71 (MQ=255)
                   cAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGAc                                               >  1:2410217/1‑71 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCaa                                                                  <  1:1186634/51‑1 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACg                                                    <  1:37096/65‑1 (MQ=255)
                    aaGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGa                                                >  1:945209/1‑69 (MQ=255)
                     aGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCa                                             <  1:3758003/71‑1 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGCCCAg                                            >  1:4653460/1‑71 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            <  1:3141356/71‑1 (MQ=255)
                      gCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            >  1:7427982/1‑71 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTc                                                               <  1:5497408/49‑1 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCa                                             <  1:7346768/67‑1 (MQ=255)
                         tCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCtt                                         <  1:4871179/71‑1 (MQ=255)
                          cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                        >  1:3942660/1‑71 (MQ=255)
                          cGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTc                                        >  1:551897/1‑71 (MQ=255)
                           gCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACtt                                                           >  1:1668234/1‑51 (MQ=255)
                           gCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACtt                                                           >  1:854368/1‑51 (MQ=255)
                            cATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCa                                                              >  1:7570305/1‑47 (MQ=255)
                               cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGc                                    >  1:6657413/1‑70 (MQ=255)
                               cAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                   >  1:913962/1‑71 (MQ=255)
                                   cACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAc                               >  1:4151023/1‑71 (MQ=255)
                                    aCTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAcc                              >  1:5989640/1‑71 (MQ=255)
                                        tGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:3287961/1‑36 (MQ=255)
                                        tGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:2850664/1‑36 (MQ=255)
                                         gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:1521097/1‑35 (MQ=255)
                                         gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:7199933/1‑35 (MQ=255)
                                         gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:7191188/1‑35 (MQ=255)
                                         gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCAc                                                             >  1:6670754/1‑35 (MQ=255)
                                         gAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttc                         >  1:1130121/1‑71 (MQ=255)
                                          aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACgg                                                   >  1:4220177/1‑44 (MQ=255)
                                          aaCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCg                                   >  1:4498948/1‑60 (MQ=255)
                                           aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTaaaaa                                                      <  1:1995980/40‑1 (MQ=255)
                                           aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttctt                       <  1:7482496/71‑1 (MQ=255)
                                           aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttctt                       <  1:3161962/71‑1 (MQ=255)
                                           aCAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                        >  1:6102334/1‑70 (MQ=255)
                                            cAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAg                                            >  1:935167/1‑49 (MQ=255)
                                               gATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGcttct                        >  1:1585079/1‑66 (MQ=255)
                                                 tGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAg                 <  1:6343939/71‑1 (MQ=255)
                                                  gCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGc                >  1:6916156/1‑71 (MQ=255)
                                                   ccGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc               >  1:3383136/1‑71 (MQ=255)
                                                   ccGGAT‑CATCCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGcc               >  1:4148253/1‑71 (MQ=255)
                                                     ggAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAAc                               >  1:3462586/1‑53 (MQ=255)
                                                      gAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCtt             >  1:7234389/1‑70 (MQ=255)
                                                      gAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTa            <  1:4422874/71‑1 (MQ=255)
                                                           aTGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAac         <  1:6229175/70‑1 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCa                                      >  1:86132/1‑38 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCa                                      >  1:3987361/1‑38 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAAcaca      >  1:4067573/1‑70 (MQ=255)
                                                              cATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     <  1:2865334/71‑1 (MQ=255)
                                                               aTCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGa                                 <  1:1459011/42‑1 (MQ=255)
                                                               aTCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     >  1:618495/1‑70 (MQ=255)
                                                               aTCCGGCAATTAACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATa    <  1:7889711/71‑1 (MQ=255)
                                                                tCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCt              >  1:2314834/1‑60 (MQ=255)
                                                                tCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACAt     >  1:5479272/1‑69 (MQ=255)
                                                                 ccGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGc                            >  1:4766842/1‑45 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:5442079/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:7749759/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:2608468/70‑1 (MQ=255)
                                                                  cGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAg  <  1:5742215/70‑1 (MQ=255)
                                                                  |                                                                     
TAGGCCTGATAAGACGCGTCAAGCGTCGCATCAGGCACTATGAACAAGATGCCGGAT‑CATGCATCCGGCAATTCACTTAAAAACGGTGACCAGCTTCCAGCGGAACCGCTTCTTCACCAGCCTTAAACACATAAG  >  W3110S.gb/4232969‑4233103

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