Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,151 C→G 25.8% intergenic (+123/+193) 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,1510CG25.8% 135.2 / 41.2 66intergenic (+123/+193)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base C (29/20);  new base G (6/11);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGT  >  minE/337084‑337221
                                                                   |                                                                      
catccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                                                     <  1:2670650/71‑1 (MQ=39)
 atccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAg                                                                    <  1:3467092/71‑1 (MQ=39)
  tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAgg                                                                   <  1:1714184/71‑1 (MQ=255)
    cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                                 >  1:2090379/1‑71 (MQ=255)
    cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                                 >  1:2095066/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: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:3485050/49‑1 (MQ=11)
                                      ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                                     <  1:2229522/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)
                                                                   cGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGt  <  1:1431032/71‑1 (MQ=39)
                                                                   |                                                                      
CATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGT  >  minE/337084‑337221

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: