Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 338,325 A→G 100% A508A (GCT→GCC copA ← copper transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE338,3250AG75.2% 25.4 / 4.2 24A508A (GCT→GCCcopAcopper transporter
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (0/18);  total (3/21)
Fisher's exact test for biased strand distribution p-value = 9.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGCGTCCCGGTTTTATCGAACACTACAGTGTCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAGCCGCCCGCCCGACGCCGGAAATAA  >  minE/338265‑338377
                                                            |                                                    
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGGGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1805894/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1180987/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:847268/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:641406/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3393140/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3390427/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3164335/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3084357/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3017619/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:2882707/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:2241754/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:2218393/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:2157303/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1828141/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1791469/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1605583/67‑1 (MQ=255)
cagCGTCCCGGTTTTATCGAACACTACGGTGTCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:1312627/67‑1 (MQ=255)
                         aCAGTGTCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTcc                                                >  1:1782407/1‑42 (MQ=255)
                         aCAGTGTCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAgccgcc                     >  1:2361376/1‑69 (MQ=255)
                         aCAGTGTCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAgccgcc                     >  1:2566846/1‑69 (MQ=255)
                            gtgtCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAGccgcccg                   <  1:633689/68‑1 (MQ=255)
                               tCGAGCGTACTGGCGCGTTGCAGCGCGTCGGCGTcc                                                <  1:3111119/36‑1 (MQ=25)
                                             gcgTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAGCCGCCCGCCCGACGCCGGAaataa  <  1:1567511/68‑1 (MQ=255)
                                             gcgTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAGCCGCCCGCCCGACGCCGGAaataa  <  1:1086870/68‑1 (MQ=255)
                                                            |                                                    
CAGCGTCCCGGTTTTATCGAACACTACAGTGTCGAGTGTACTGGCGCGTTGCAGCGCGTCAGCGTCCCGCACCAGCACGCCAAACTCAGCCGCCCGCCCGACGCCGGAAATAA  >  minE/338265‑338377

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: