Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,687,235 G→A 100% P129S (CCG→TCG)  sseB ← rhodanase‑like enzyme, sulfur transfer from thiosulfate

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,687,2350GA87.3% 52.6 / ‑0.0 31P129S (CCG→TCG) sseBrhodanase‑like enzyme, sulfur transfer from thiosulfate
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (16/11);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCG  >  minE/1687199‑1687290
                                    |                                                       
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTGCCTCTTCACCAATCAACAAACTGAttt                         <  1:1244807/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt                         <  1:3115155/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt                         <  1:3263548/69‑1 (MQ=255)
aTTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGAttt                         <  1:56674/69‑1 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGGGAATTCCCTTCTTCACCAATc                                      >  1:223189/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:3225919/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:3261280/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2946256/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2926496/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2923367/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2648455/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2255422/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2159651/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:2103404/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:1980290/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:1940686/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:185022/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:1686313/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:1531145/1‑44 (MQ=255)
            ccAGTACTTCCTGGCTACTCAGCGAATTCCCTTCTTCACCAATc                                      >  1:81742/1‑44 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:367206/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:3367441/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:860812/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:3212823/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:3170187/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:2907634/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:2636104/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:2633094/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:2450489/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:1418632/60‑1 (MQ=255)
                                aGCGAATTCCCTTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAACTCTTTACCg  <  1:1397517/60‑1 (MQ=255)
                                    |                                                       
ATTCACCGCCTTCCAGGATTTCCTGGCTGCTCAGCGGATTTCCCTCTTCACCAATCAACAAACTGATTTCACGCGGCATAAATTCTTTACCG  >  minE/1687199‑1687290

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