Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,146 T→G 9.1% intergenic (+118/+198) ybaQ → / ← copA predicted DNA‑binding transcriptional regulator/copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE337,1460TG 9.1% 157.5 / 11.4 66intergenic (+118/+198)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (33/27);  new base G (2/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAC  >  minE/337076‑337211
                                                                      |                                                                 
gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCt                                                                   >  1:1871703/1‑71 (MQ=39)
    ccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                                               >  1:3723918/1‑71 (MQ=39)
     cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                                               <  1:2398076/70‑1 (MQ=40)
        catccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                                           <  1:2670650/71‑1 (MQ=39)
         atccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAg                                                          <  1:3467092/71‑1 (MQ=39)
          tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAgg                                                         <  1:1714184/71‑1 (MQ=255)
           ccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCt                                                                   >  1:2401576/1‑60 (MQ=255)
            cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                       >  1:2095066/1‑71 (MQ=255)
            cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                       >  1:2090379/1‑71 (MQ=255)
               gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                                    >  1:1527728/1‑71 (MQ=255)
                cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                   <  1:1668501/71‑1 (MQ=255)
                 aCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGtata                                                  <  1:1462145/71‑1 (MQ=255)
                 aCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATa                                                  <  1:2735455/71‑1 (MQ=255)
                   cATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                                                >  1:972171/1‑71 (MQ=255)
                     ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                              >  1:922520/1‑71 (MQ=255)
                     ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                              >  1:1698064/1‑71 (MQ=255)
                     ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                              >  1:2047034/1‑71 (MQ=255)
                      tGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCg                                             >  1:2394883/1‑71 (MQ=255)
                       gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                            >  1:2434305/1‑71 (MQ=255)
                       gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                            >  1:2412190/1‑71 (MQ=255)
                        ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                 >  1:3215591/1‑65 (MQ=255)
                        ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                 >  1:2849537/1‑65 (MQ=255)
                        ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                           <  1:1626123/71‑1 (MQ=39)
                          tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                         >  1:3633491/1‑71 (MQ=39)
                           gATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                            >  1:2083478/1‑67 (MQ=255)
                           gATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                                        >  1:862347/1‑71 (MQ=39)
                            aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcg                                       <  1:3194471/71‑1 (MQ=39)
                            aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                                        <  1:2836967/70‑1 (MQ=40)
                               cGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                   <  1:1553886/56‑1 (MQ=255)
                                 aCGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                  <  1:2175251/71‑1 (MQ=39)
                                     ttACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatc                                >  1:644828/1‑69 (MQ=40)
                                       gcgcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                             <  1:2135614/67‑1 (MQ=11)
                                       gcgcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                             <  1:2559889/67‑1 (MQ=11)
                                       aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                             >  1:1864280/1‑70 (MQ=40)
                                         gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                          <  1:2193354/70‑2 (MQ=11)
                                         gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                          <  1:3312398/70‑2 (MQ=11)
                                         gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatc                                >  1:2263847/2‑65 (MQ=11)
                                         gcGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCgg                                                     >  1:108496/2‑44 (MQ=11)
                                         gcGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCgcatccggcag                          >  1:3188136/2‑70 (MQ=11)
                                           gTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcagtt                        <  1:2350782/71‑4 (MQ=11)
                                             cTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggcatcc                        >  1:75634/1‑69 (MQ=40)
                                              ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                           <  1:2229522/49‑1 (MQ=11)
                                              ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                           <  1:3246850/49‑1 (MQ=11)
                                              ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                           <  1:3485050/49‑1 (MQ=11)
                                              ttATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAg                     >  1:245724/1‑71 (MQ=255)
                                               tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgccgc                                   <  1:2620764/56‑1 (MQ=11)
                                               tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatcc                               <  1:2685157/60‑1 (MQ=11)
                                                aTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGa                      <  1:1594343/68‑1 (MQ=255)
                                                aTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAg                     <  1:814902/69‑1 (MQ=255)
                                                aTCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa                   >  1:3421221/1‑71 (MQ=255)
                                                 tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcat                                 >  1:3670952/1‑56 (MQ=11)
                                                 tCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGa                      <  1:2008335/67‑1 (MQ=255)
                                                 tCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa                   <  1:1098202/70‑1 (MQ=255)
                                                  cAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAc                  >  1:3127716/1‑70 (MQ=255)
                                                   atgCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCg                                             <  1:1631977/40‑1 (MQ=14)
                                                      ccTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTg             <  1:70359/71‑1 (MQ=255)
                                                        tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                                           <  1:1105033/39‑1 (MQ=40)
                                                         aCAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCTAg                     >  1:491038/1‑60 (MQ=255)
                                                         aCAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAcc                 <  1:783144/64‑1 (MQ=255)
                                                          cAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                          >  1:2898298/1‑53 (MQ=11)
                                                          cAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                          >  1:1196675/1‑53 (MQ=11)
                                                          cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                         >  1:1801825/1‑39 (MQ=40)
                                                          cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggca                           <  1:1610931/53‑1 (MQ=35)
                                                           aaaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                                        >  1:3417190/1‑39 (MQ=40)
                                                            aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGAt       >  1:115674/1‑71 (MQ=255)
                                                            aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGAt       >  1:999915/1‑71 (MQ=255)
                                                                 gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAc  >  1:2354274/1‑71 (MQ=255)
                                                                 gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAc  >  1:633180/1‑71 (MQ=255)
                                                                      |                                                                 
GCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAC  >  minE/337076‑337211

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: