Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,149 A→T 9.0% intergenic (+121/+195) 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,1490AT 9.0% 161.3 / 9.3 67intergenic (+121/+195)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/4);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

CCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGC  >  minE/337080‑337219
                                                                     |                                                                      
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)
        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:1698064/1‑71 (MQ=255)
                 ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                                      >  1:922520/1‑71 (MQ=255)
                 ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                                                      >  1:2047034/1‑71 (MQ=255)
                  tGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCg                                                     >  1:2394883/1‑71 (MQ=255)
                   gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                                    >  1:2412190/1‑71 (MQ=255)
                   gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                                    >  1:2434305/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:3485050/49‑1 (MQ=11)
                                          ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                                   <  1:3246850/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)
                                                                     |                                                                      
CCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGC  >  minE/337080‑337219

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: