Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,128 G→T 20.8% intergenic (+100/+216) 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,1280GT20.8% 185.1 / 23.5 67intergenic (+100/+216)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base T (4/10);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAC  >  minE/337064‑337195
                                                                |                                                                   
ggCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTAc                                                               <  1:2405529/71‑1 (MQ=39)
 gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACa                                                              <  1:3729263/71‑1 (MQ=39)
 gCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGAAGCGACGCTTACCCGTCTTATCAGGCCTACa                                                              <  1:1890017/71‑1 (MQ=255)
  cGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaa                                                             <  1:3936250/71‑1 (MQ=39)
      cACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                          <  1:2951325/70‑1 (MQ=39)
      cACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACt                                                          <  1:2942945/70‑1 (MQ=39)
       aCGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACa                                                              <  1:3926497/65‑1 (MQ=35)
        cgccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACtt                                                         <  1:3781080/69‑1 (MQ=39)
        cgccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                       <  1:3214427/71‑1 (MQ=39)
            gcatccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCt                                                   >  1:1871703/1‑71 (MQ=39)
                ccggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                               >  1:3723918/1‑71 (MQ=39)
                 cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAc                                               <  1:2398076/70‑1 (MQ=40)
                    catccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTAc                                                               >  1:2806241/1‑51 (MQ=35)
                    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:2090379/1‑71 (MQ=255)
                        cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                       >  1:2095066/1‑71 (MQ=255)
                         gAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                       <  1:1588869/54‑1 (MQ=255)
                          aGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                       <  1:2421216/53‑1 (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:2047034/1‑71 (MQ=255)
                                 ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                              >  1:922520/1‑71 (MQ=255)
                                 ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGc                              >  1:1698064/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)
                                        aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                       >  1:3543244/1‑39 (MQ=255)
                                        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:3312398/70‑2 (MQ=11)
                                                     gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag          <  1:2193354/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:3246850/49‑1 (MQ=11)
                                                          ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                           <  1:3485050/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)
                                                                |                                                                   
GGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAC  >  minE/337064‑337195

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: