Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 337,151 C→G 46.1% 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,1510CG46.1% 62.5 / 84.5 67intergenic (+123/+193)ybaQ/copApredicted DNA‑binding transcriptional regulator/copper transporter
Reads supporting (aligned to +/- strand):  ref base C (15/21);  new base G (7/23);  total (22/45)
Fisher's exact test for biased strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

CGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCT  >  minE/337081‑337214
                                                                      |                                                               
cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAcc                                                                 >  1:1160387/1‑71 (MQ=39)
cggcatccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAAcc                                                                 >  1:1167226/1‑71 (MQ=39)
   catccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                                              <  1:4643591/71‑1 (MQ=39)
     tccgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAgg                                                            <  1:4124728/71‑1 (MQ=39)
       cgAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGcc                                                          >  1:3645992/1‑71 (MQ=39)
          gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                                       >  1:2769352/1‑71 (MQ=255)
          gCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGa                                                       >  1:2203401/1‑71 (MQ=255)
           cACCGTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                      <  1:841036/71‑1 (MQ=255)
           cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                      <  1:4611801/71‑1 (MQ=255)
           cACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                      <  1:3458981/71‑1 (MQ=255)
               aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                      <  1:1673968/67‑1 (MQ=255)
               aTTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                  >  1:3156300/1‑71 (MQ=255)
                  gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATaa                                                    >  1:1919618/1‑66 (MQ=255)
                  gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                                                   >  1:5474537/1‑67 (MQ=255)
                  gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAg                                                   >  1:3132075/1‑67 (MQ=255)
                  gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCg                                                <  1:1366499/70‑1 (MQ=40)
                  gCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGt                                               >  1:1283872/1‑71 (MQ=39)
                     tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                                              >  1:5067364/1‑53 (MQ=255)
                     tGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTa                                                              >  1:2726734/1‑53 (MQ=255)
                       aTGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGAt                                                      <  1:31433/59‑1 (MQ=255)
                          cGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTAcgccg                                       <  1:2408332/71‑1 (MQ=39)
                             cGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcat                                    <  1:1080292/71‑1 (MQ=39)
                               cTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatcc                                  <  1:277113/71‑1 (MQ=39)
                               cTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatcc                                  >  1:4314238/1‑71 (MQ=39)
                                  gcgcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                                <  1:709523/67‑1 (MQ=11)
                                  aCCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccgg                                >  1:2097000/1‑70 (MQ=40)
                                    gcGTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAg                                                             <  1:3841883/38‑1 (MQ=11)
                                      gTCTTATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                                           >  1:2300586/1‑38 (MQ=14)
                                       tCTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                             <  1:5343446/68‑2 (MQ=11)
                                       tCTTATCATGCCTACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGttt                                             <  1:2583956/52‑1 (MQ=11)
                                       tCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                     <  1:4141989/60‑1 (MQ=35)
                                        cTTATCATGCCTACAAACTTGTGCCTGAACGGTAg                                                             <  1:2725590/35‑1 (MQ=14)
                                        cTTATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCg                                                >  1:1829399/1‑48 (MQ=11)
                                         ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCg                                                <  1:3592969/47‑1 (MQ=11)
                                         ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                                <  1:3746486/63‑1 (MQ=11)
                                         ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                                <  1:4532308/63‑1 (MQ=11)
                                         ttATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccgg                                <  1:712778/63‑1 (MQ=11)
                                          tATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                                           >  1:2986557/1‑34 (MQ=14)
                                          tATCATGCCTACAAACTTGTGCCTGAACGGTAgga                                                           >  1:3834808/1‑34 (MQ=14)
                                          tATCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGAt                                                      <  1:2194573/40‑1 (MQ=14)
                                           aTCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                                                  <  1:3865459/43‑1 (MQ=14)
                                           aTCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAgg                                                  <  1:1529795/43‑1 (MQ=14)
                                            tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                                              <  1:2434942/46‑1 (MQ=11)
                                            tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                                              <  1:3298669/46‑1 (MQ=11)
                                            tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                                              <  1:2441412/46‑1 (MQ=11)
                                            tCATGCCTACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgca                                     <  1:2503124/55‑1 (MQ=11)
                                               tgCCTACAAACTTGTGCCTGAACGGTAGGACGGATaa                                                    >  1:1375014/2‑37 (MQ=14)
                                               tgCCTACAAACTTGTGCCTGAACGGTAGGACGGATaa                                                    >  1:295680/2‑37 (MQ=14)
                                                gCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCAtt                 >  1:123370/1‑71 (MQ=255)
                                                   tACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGtt                                              <  1:82049/39‑1 (MQ=14)
                                                   tACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggc                               <  1:3542284/54‑1 (MQ=11)
                                                   tACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggc                               <  1:4483043/54‑1 (MQ=11)
                                                   tACAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggc                               <  1:5498161/54‑1 (MQ=11)
                                                   tACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                  <  1:352773/35‑1 (MQ=40)
                                                   tACAAACTTGTGCCTGAACCGTAGGCCGGATAAgg                                                  <  1:667819/35‑1 (MQ=40)
                                                   tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggcatccg                          <  1:3597628/59‑1 (MQ=255)
                                                   tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                     <  1:3936632/48‑1 (MQ=38)
                                                   tACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                     <  1:18786/48‑1 (MQ=38)
                                                   tACAAACTTGTGCCGGATCTGTAGGCCGGATAAGGCGTTTAtgccg                                       <  1:3856010/46‑1 (MQ=11)
                                                   tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGtt                                              <  1:1356805/39‑1 (MQ=14)
                                                   tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTAtgccg                                       <  1:2684260/46‑1 (MQ=11)
                                                   tACAAACTTGTGCCGGATCGGTAGGCCGGATAAGGCGTTTAtgccg                                       <  1:3556127/46‑1 (MQ=11)
                                                     cAAACTTGTGCCTGAACGGTAGGACGGATAAGGCGTTCACGCCgcatccggcag                             >  1:376778/1‑53 (MQ=11)
                                                     cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgcatccggca                              >  1:5661101/1‑53 (MQ=35)
                                                     cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt             <  1:2854365/70‑1 (MQ=255)
                                                     cAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCt             >  1:3725036/1‑70 (MQ=255)
                                                        aCTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATg         <  1:2371740/71‑1 (MQ=255)
                                                          ttGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCgca                                     <  1:4301067/41‑1 (MQ=32)
                                                               ccTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCt  <  1:5197011/71‑1 (MQ=255)
                                                                      |                                                               
CGGCATCCGAGCACCATTGCCTGATGCGACGCTTACCCGTCTTATCAGGCCTACAAACTTGTGCCTGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATCCGAGCACCATTGCCTGATGCGACGCT  >  minE/337081‑337214

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: