Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,845,199 G→T 100% P266T (CCG→ACG)  cysD ← sulfate adenylyltransferase, subunit 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,845,1990GT95.2% 211.8 / 3.3 63P266T (CCG→ACG) cysDsulfate adenylyltransferase, subunit 2
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base T (60/0);  total (63/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCG  >  minE/1845162‑1845237
                                     |                                      
cTGGTGGTGGAAACCATCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:589806/1‑54 (MQ=38)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:2024884/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1666009/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1667005/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1736514/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:174443/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1882187/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:192077/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:2020557/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:2023605/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1014465/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:2131218/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:316567/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:468326/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:818667/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:861330/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:980419/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:987230/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1294934/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1213254/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1430220/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1240365/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1151302/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1364613/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1438212/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1129721/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:1116426/1‑54 (MQ=255)
cTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCAttt                        >  1:127980/1‑54 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCa                           >  1:1511431/1‑44 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:501115/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:230384/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:2146495/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:259846/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:2143633/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:347862/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:473618/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1003293/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:55191/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:704519/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:742144/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1145305/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:826213/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:91343/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1644342/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1808929/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1652435/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1553866/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1704891/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1746538/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1764059/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1766708/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:2129690/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:186527/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1907753/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1413874/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1981456/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1985689/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:131907/1‑67 (MQ=255)
       tggAAACCAGCATCTCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:2117218/1‑67 (MQ=255)
       tggAAACCAGCATATCTTCGATGATTTCCGTCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAg    >  1:1871061/1‑67 (MQ=255)
           aaCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCg  >  1:1767756/1‑65 (MQ=255)
           aaCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCg  >  1:1690303/1‑65 (MQ=255)
           aaCCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCg  >  1:1545673/1‑65 (MQ=255)
                                     |                                      
CTGGTGGTGGAAACCAGCATCTCTTCGATGATTTCCGGCAGTGTTTGTGCATTTGACTCCACCGCACCGGTCAGCG  >  minE/1845162‑1845237

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