Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 753,377 T→A 100% S348S (TCA→TCT ycfT ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE753,3770TA92.6% 56.4 / ‑2.0 27S348S (TCA→TCTycfTpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base A (25/0);  total (25/2)
Fisher's exact test for biased strand distribution p-value = 2.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAC  >  minE/753326‑753445
                                                   |                                                                    
cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGTAGCGGGCAGAGAAGGCGGGGAGaac                                                       >  1:780783/1‑67 (MQ=255)
cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGTAGCGGGCAGAGAAGGCGGGGAGaac                                                       >  1:3285640/1‑67 (MQ=255)
cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGTAGCGGGCAGAGAAGGCGGGGAGaac                                                       >  1:2253004/1‑67 (MQ=255)
cGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGTAGCGGGCAGAGAAGGCGGGGAAaac                                                       >  1:1702043/1‑67 (MQ=255)
      ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATc                                               <  1:2935673/69‑1 (MQ=255)
             gCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGa                                          <  1:2452871/67‑1 (MQ=255)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:265309/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:253913/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:272477/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2780070/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2817096/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2832395/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:3080973/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:378444/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:431862/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:529885/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:706690/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2610857/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2546155/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:1104092/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2310699/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2217346/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:2086372/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:1959397/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:1816435/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:1502647/1‑47 (MQ=39)
                                               gCAGAGAAGGCGGGGAGAACAACAGATCGCTGAATGCGCGCTGTGaa                            >  1:1433549/1‑47 (MQ=39)
                                                   tGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc  >  1:2515217/1‑69 (MQ=255)
                                                   tGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc  >  1:1801884/1‑69 (MQ=255)
                                                   tGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc  >  1:1586980/1‑69 (MQ=255)
                                                   tGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc  >  1:1191232/1‑69 (MQ=255)
                                                   |                                                                    
CGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAC  >  minE/753326‑753445

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