Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,128 G→T 17.0% 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,1280GT17.0% 133.2 / 10.8 47intergenic (+100/+216)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base T (3/5);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.52e-01

GGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAT  >  minE/337058‑337198
                                                                      |                                                                      
ggATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAg                                                                        <  1:2520453/71‑1 (MQ=39)
       gCGTTCACTCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACa                                                                 <  1:1385624/71‑1 (MQ=255)
          ttCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACaa                                                                <  1:1414843/69‑1 (MQ=40)
            cACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACtt                                                            >  1:208545/1‑71 (MQ=39)
              cgccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                          <  1:2774818/71‑1 (MQ=39)
               gccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                          <  1:2290489/70‑1 (MQ=40)
                ccgcATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGc                                                        <  1:1908381/71‑1 (MQ=39)
                         gcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCg                                                >  1:1304859/1‑70 (MQ=40)
                             ccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                           >  1:1036875/1‑71 (MQ=39)
                              cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                          >  1:3246959/1‑55 (MQ=255)
                              cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                          >  1:300237/1‑71 (MQ=255)
                              cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                           <  1:286048/70‑1 (MQ=255)
                                aGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCgg                                        <  1:2826427/71‑1 (MQ=255)
                                 gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAcc                                                 >  1:4136104/1‑61 (MQ=255)
                                 gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCg                                                >  1:1792256/1‑62 (MQ=255)
                                 gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                       <  1:4024984/71‑1 (MQ=255)
                                  cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtg                                                           >  1:782387/1‑50 (MQ=255)
                                  cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtg                                                           >  1:1598147/1‑50 (MQ=255)
                                  cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                      <  1:1640671/71‑1 (MQ=255)
                                   aCCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGc                                           <  1:3707386/65‑1 (MQ=255)
                                      aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                  >  1:3480069/1‑71 (MQ=255)
                                      aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                                   >  1:494883/1‑70 (MQ=255)
                                       ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTAc                                                                  >  1:3853705/1‑38 (MQ=40)
                                       ttGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                       <  1:2518887/65‑1 (MQ=255)
                                        ttccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgt                                                          <  1:3213249/43‑1 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGcc                                                       >  1:2589232/1‑47 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGcc                                                       >  1:1331733/1‑47 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                              >  1:4001905/1‑56 (MQ=255)
                                         gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATa                                     >  1:2111306/1‑65 (MQ=255)
                                          ccTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                              <  1:2578733/71‑1 (MQ=39)
                                           cTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                              >  1:123200/1‑70 (MQ=40)
                                             gATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTtgtg                                                         <  1:441281/41‑1 (MQ=255)
                                                 cGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGtt                              <  1:4044008/64‑1 (MQ=40)
                                                   aCGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                     <  1:2595084/71‑1 (MQ=39)
                                                        tACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                >  1:2480968/1‑71 (MQ=39)
                                                        tACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                >  1:3004984/1‑71 (MQ=39)
                                                            cGTCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTATGCCgcatccggcagt            <  1:3735837/71‑3 (MQ=11)
                                                              tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgccg                       <  1:4045133/58‑1 (MQ=11)
                                                              tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgccg                       <  1:541952/58‑1 (MQ=11)
                                                              tCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCg                                         <  1:1005731/40‑1 (MQ=18)
                                                                ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATaa                                    >  1:3022205/1‑43 (MQ=14)
                                                                ttATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGa                                       >  1:129888/1‑40 (MQ=18)
                                                                ttATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAg        >  1:2326016/1‑71 (MQ=255)
                                                                 tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGttc                             >  1:2213155/1‑48 (MQ=11)
                                                                   tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCAcgcc                        <  1:1979013/52‑1 (MQ=11)
                                                                   tCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAc     <  1:2318689/71‑1 (MQ=255)
                                                                    cAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGc       <  1:2027619/68‑1 (MQ=255)
                                                                     aGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCAcc    <  1:2605495/70‑1 (MQ=255)
                                                                      ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa      >  1:2810483/1‑67 (MQ=255)
                                                                      ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCa      >  1:2805773/1‑67 (MQ=255)
                                                                      ggCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAt  <  1:4016339/71‑1 (MQ=255)
                                                                      |                                                                      
GGATAAGGCGTTCACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAT  >  minE/337058‑337198

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: