Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,171 T→C 21.6% intergenic (+143/+173) 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,1710TC21.6% 114.4 / 23.6 51intergenic (+143/+173)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base T (24/16);  new base C (4/7);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACC  >  minE/337102‑337240
                                                                     |                                                                     
tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTa                                                                      >  1:3633491/1‑71 (MQ=39)
 gATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                                                                     >  1:862347/1‑71 (MQ=39)
  aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcg                                                                    <  1:3194471/71‑1 (MQ=39)
  aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAc                                                                     <  1:2836967/70‑1 (MQ=40)
       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)
               gcGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCgcatccggcag                                                       >  1:3188136/2‑70 (MQ=11)
                 gTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcagtt                                                     <  1:2350782/71‑4 (MQ=11)
                   cTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggcatcc                                                     >  1:75634/1‑69 (MQ=40)
                    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)
                            ccTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTg                                          <  1:70359/71‑1 (MQ=255)
                               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:999915/1‑71 (MQ=255)
                                  aaCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGAt                                    >  1:115674/1‑71 (MQ=255)
                                       gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAc                               >  1:2354274/1‑71 (MQ=255)
                                       gtgCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGAc                               >  1:633180/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:3248760/71‑1 (MQ=39)
                                                     ggCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTa                 <  1:1966153/71‑1 (MQ=39)
                                                        cGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTAt                >  1:785950/1‑69 (MQ=40)
                                                         ggATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAg             >  1:3810200/1‑71 (MQ=39)
                                                         ggATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAg             <  1:1935317/71‑1 (MQ=39)
                                                          gATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGt                     <  1:503095/62‑1 (MQ=40)
                                                          gATAAGGCGTTTAAGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAgg            >  1:2663028/1‑71 (MQ=255)
                                                             aaGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt                                       >  1:3910981/1‑41 (MQ=38)
                                                                 cGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACa      >  1:3522680/1‑70 (MQ=40)
                                                                    ttACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAAcc  >  1:1558527/1‑71 (MQ=39)
                                                                     tACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGc           >  1:1275691/1‑61 (MQ=37)
                                                                     tACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGc           >  1:2809601/1‑61 (MQ=37)
                                                                     |                                                                     
TGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACAAACC  >  minE/337102‑337240

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: