Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,875,977 A→G 100% N262S (AAC→AGC)  ygcG → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,875,9770AG78.3% 40.4 / 7.7 23N262S (AAC→AGC) ygcGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base G (0/18);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 2.97e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGC  >  minE/1875912‑1875989
                                                                 |            
tATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGacaac             >  1:1437501/1‑67 (MQ=255)
tATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGacaac             >  1:616098/1‑67 (MQ=255)
      cTTTACCTGGTTATGTGTGCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1788942/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAATAAGCTTCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1553180/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:788810/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:765318/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:586740/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:384727/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:380666/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:289973/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:2242060/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:2135429/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:2133107/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:2103686/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1920379/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1747079/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1712766/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:166825/66‑1 (MQ=255)
      cTTTACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:1486280/66‑1 (MQ=255)
       tttACCTGGTTATGTGTCCCACAGAAAAAAGCATCAGGTAGAAGTTATCATTCAGACAGcggcgg        <  1:223426/65‑1 (MQ=255)
           ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGc  >  1:479772/1‑67 (MQ=255)
           ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGc  >  1:481823/1‑67 (MQ=255)
           ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGc  >  1:484432/1‑67 (MQ=255)
           ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGc  >  1:1473017/1‑67 (MQ=255)
           ccTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGc  >  1:1239255/1‑67 (MQ=255)
                                                                 |            
TATGGTCTTTACCTGTTTATGTGTACTTCAGAAAAGAGCATCTGGTAGAAGTTATCATTCAGACAACAGCGGTTCCGC  >  minE/1875912‑1875989

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: