Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,157 C→A 16.6% intergenic (+129/+187) 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,1570CA16.6% 181.0 / 17.0 66intergenic (+129/+187)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base C (31/24);  new base A (4/7);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

CGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTA  >  minE/337088‑337225
                                                                     |                                                                    
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:2047034/1‑71 (MQ=255)
         ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                                            >  1:1698064/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:2849537/1‑65 (MQ=255)
            ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                               >  1:3215591/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:2559889/67‑1 (MQ=11)
                           gcgcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                                           <  1:2135614/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:3485050/49‑1 (MQ=11)
                                  ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                                         <  1:3246850/49‑1 (MQ=11)
                                  ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                                         <  1:2229522/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:633180/1‑71 (MQ=255)
                                                     gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAc                >  1:2354274/1‑71 (MQ=255)
                                                           gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGAt                     <  1:2182144/60‑1 (MQ=255)
                                                             aCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTAcgc        >  1:2557321/1‑71 (MQ=40)
                                                               cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGt      <  1:1431032/71‑1 (MQ=39)
                                                                 tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCt    >  1:3941135/1‑71 (MQ=39)
                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa  >  1:3660796/1‑71 (MQ=39)
                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa  <  1:1966153/71‑1 (MQ=39)
                                                                   ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa  <  1:3248760/71‑1 (MQ=39)
                                                                     |                                                                    
CGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTA  >  minE/337088‑337225

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: