Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE585,7550GT36.0% 78.0 / 44.1 50A47E (GCG→GAG) cydCfused cysteine transporter subunits and membrane component and ATP‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (0/32);  new base T (18/0);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 5.54e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGTCCGGCAACCCCCGCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTCA  >  minE/585734‑585815
                     |                                                            
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGCACCAGCCGGAAAGTGTCAACAGAcc                      >  1:207829/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGa                        >  1:108139/1‑60 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:795036/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:553211/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:41995/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:660088/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:352129/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:350673/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:257895/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:217210/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:174722/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:793408/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:847616/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:101520/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:1025716/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAcc                      >  1:1048941/1‑62 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAc                       >  1:857509/1‑61 (MQ=255)
aGTCCGGCAACCCCCGCAACCTCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGAc                       >  1:1045184/1‑61 (MQ=255)
               gCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGagcag       <  1:846875/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1000482/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:635455/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:637368/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:683503/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:955450/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:724962/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:776321/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:948663/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:933612/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:893960/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:963429/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:795075/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:81447/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:158263/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1013805/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1034000/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1089719/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1126619/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1140101/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1158165/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:1178452/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:129512/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:546705/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:159740/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:173974/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:263945/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:348550/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:394690/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:465601/62‑1 (MQ=255)
                    cGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACAGATACTGGCGAGCAGCGTca  <  1:666590/62‑1 (MQ=255)
                    cGCTGAGGCCGAAAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:3680/62‑1 (MQ=255)
                     gCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTca  <  1:13636/61‑1 (MQ=255)
                     |                                                            
AGTCCGGCAACCCCCGCAACCGCTGAGGCCGAGAGGAACCAGCCGGAAAGTGTCAACAGACCGATACTGGCGAGCAGCGTCA  >  minE/585734‑585815

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: