Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 215,801 A→G 16.0% E178G (GAG→GGG)  nlpE → lipoprotein involved with copper homeostasis and adhesion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb215,8010AG16.0% 117.3 / 16.3 50E178G (GAG→GGG) nlpElipoprotein involved with copper homeostasis and adhesion
Reads supporting (aligned to +/- strand):  ref base A (23/19);  new base G (2/6);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

ATGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCAGTAGA  >  W3110S.gb/215741‑215864
                                                            |                                                               
aTGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGGGCTGGAGc                                                         <  1:174519/69‑1 (MQ=255)
aTGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGt                                                       <  1:2231957/71‑1 (MQ=255)
 tGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCg                                                        >  1:3168733/1‑69 (MQ=255)
 tGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCg                                                        >  1:3785159/1‑69 (MQ=255)
 tGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCg                                                        >  1:4266440/1‑69 (MQ=255)
  gcggcgACATTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGGGCTGGAGCGTAg                                                     >  1:2448796/1‑71 (MQ=255)
    ggcgACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCt                                                   <  1:709965/71‑1 (MQ=255)
    ggcgACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCt                                                   <  1:980507/71‑1 (MQ=255)
     gcgACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTa                                                  >  1:4365212/1‑71 (MQ=255)
     gcgACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTa                                                  >  1:2685187/1‑71 (MQ=255)
      cgACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGTAGCGTAGCTAc                                                 <  1:4273491/71‑1 (MQ=255)
       gACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGGGCTGGAGCGTAGCTAcc                                                <  1:4079697/71‑1 (MQ=255)
         ccTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCt                                               >  1:4638515/1‑70 (MQ=255)
          cTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTgg                                             >  1:4449624/1‑71 (MQ=255)
            tCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTgg                                             >  1:3473543/1‑69 (MQ=255)
            tCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCt                                           >  1:588270/1‑71 (MQ=255)
            tCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCt                                           >  1:2854530/1‑71 (MQ=255)
             cACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTa                                                  >  1:1852018/1‑63 (MQ=255)
                tGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTAc                                                 >  1:1389753/1‑61 (MQ=255)
                  aTTGCGCGACCGGAAAAGGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTAc                                                 >  1:2837376/1‑59 (MQ=255)
                  aTTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGGGCTGGAGCGTAGCTACCTGGCTgcgcg                                      >  1:3930923/1‑70 (MQ=255)
                  aTTGCGCGACCGGAAAACGTTTCATGGTAGCCAATAACGCAGAGCTg                                                             >  1:2625813/1‑47 (MQ=255)
                    tGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCgg                                   >  1:3577392/1‑71 (MQ=255)
                      cgcgACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGt                                  <  1:2225605/70‑1 (MQ=255)
                       gcgACCGGAACACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTca                                <  1:138795/71‑1 (MQ=255)
                         gACCTGAAAACGTTTCATGGTAGCGAATAACGCAGGGCTg                                                             <  1:1965259/40‑1 (MQ=38)
                           ccGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAg                                                     >  1:489659/1‑46 (MQ=255)
                               aaaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGaaa                        <  1:3077537/71‑1 (MQ=255)
                               aaaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGaaa                        <  1:3969235/71‑1 (MQ=255)
                                aaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTAc                                                 >  1:1955634/1‑45 (MQ=255)
                                aaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTAc                                                 >  1:3859598/1‑45 (MQ=255)
                                aaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGaaaa                       <  1:2114407/71‑1 (MQ=255)
                                aaaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGaaa                        >  1:1239817/1‑70 (MQ=255)
                                 aaCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCt                                                   <  1:327090/42‑1 (MQ=255)
                                  aCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAAc                     <  1:4206069/71‑1 (MQ=255)
                                  aCGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAAc                     >  1:2128004/1‑71 (MQ=255)
                                   cGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGc                                                    >  1:4431667/1‑39 (MQ=255)
                                   cGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAAcc                    <  1:2090114/71‑1 (MQ=255)
                                     tttCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCgg                  <  1:2195375/71‑1 (MQ=255)
                                         aTGGTAGCGAATAACGCAGGGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGtt              <  1:3510695/71‑1 (MQ=255)
                                         aTGGTAGCGAATAACGCAGGGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGtt              <  1:658292/71‑1 (MQ=255)
                                         aTGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGtt              <  1:2824787/71‑1 (MQ=255)
                                          tGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTa             <  1:2503295/71‑1 (MQ=255)
                                            gTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACt           <  1:2467452/71‑1 (MQ=255)
                                             tAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTg          >  1:3029579/1‑71 (MQ=255)
                                                cGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCa       <  1:278055/71‑1 (MQ=255)
                                                    tAACGCAGGGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCAGTAg   <  1:674006/71‑1 (MQ=255)
                                                    tAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTca                                <  1:3434423/42‑1 (MQ=255)
                                                      aCGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCAGTAg   >  1:267575/1‑69 (MQ=255)
                                                      aCGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCAGTAGa  <  1:2811496/70‑1 (MQ=255)
                                                            |                                                               
ATGCGGCGACCTTCACTGATTGCGCGACCGGAAAACGTTTCATGGTAGCGAATAACGCAGAGCTGGAGCGTAGCTACCTGGCTGCGCGCGGTCACAGTGAAAAACCGGTGTTACTGTCAGTAGA  >  W3110S.gb/215741‑215864

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: