Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,139,906 T→C 100% L428L (TTA→CTA)  ydjN → predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,139,9060TC76.0% 33.2 / 5.7 29L428L (TTA→CTA) ydjNpredicted transporter
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base C (22/0);  total (22/7)
Fisher's exact test for biased strand distribution p-value = 6.41e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGCCAGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGTTAGTGGCTCGATGACAGCTGGCACGCTGACCAGCCAGTGGCTGAAGCAAA  >  minE/1139840‑1139963
                                                                  |                                                         
cTGCCAGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTa                                                         <  1:960776/69‑1 (MQ=255)
cTGCCAGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTa                                                         <  1:3204059/69‑1 (MQ=255)
     aGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGt                                                    <  1:1294904/69‑1 (MQ=255)
     aGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGt                                                    <  1:180550/69‑1 (MQ=255)
     aGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGt                                                    <  1:2302603/69‑1 (MQ=255)
     aGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGt                                                    <  1:2466174/69‑1 (MQ=255)
           ccTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTa                                                         <  1:2271445/58‑1 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGa                                         >  1:1404525/1‑63 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:1094070/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:1264685/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:333586/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:1671642/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:3072063/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:2978059/1‑64 (MQ=255)
                      tGTTAATCTCCGTTGAACCGCTTATCGACATGGGTCGTACGGCGCTAAACGTTAGTGGCTCGAt                                        >  1:2932763/1‑64 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGAccag                   >  1:123257/3‑52 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:2715353/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:269490/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:2489512/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:3011513/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:1917842/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:3104079/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:3105497/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:1502920/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:1469723/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:3352660/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaaa  >  1:890815/3‑69 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCaa   >  1:2985455/3‑68 (MQ=255)
                                                       gtcGTACGGCGCTAAACGTTAGTGGCTCGATGACCGCTGGCACGCTGACCAGCCAGTGGCTGAAGCa    >  1:1461016/3‑67 (MQ=255)
                                                                  |                                                         
CTGCCAGTAACCCTGGTGGCGCTGTTAATCTCCGTTGAACCGCTTATCGACATGGGCCGTACGGCGTTAAACGTTAGTGGCTCGATGACAGCTGGCACGCTGACCAGCCAGTGGCTGAAGCAAA  >  minE/1139840‑1139963

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: