Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 1,592,974 C→T S300N (AGC→AAC)  cysA ← sulfate/thiosulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,592,9740CT98.5% 187.0 / ‑6.5 65S300N (AGC→AAC) cysAsulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (0/64);  minor base G (0/1);  total (0/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGG  >  minE/1592959‑1593021
               |                                               
tAGTGACCTTTCGTGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3004152/62‑1 (MQ=255)
tAGTGACCTTTCGGTTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3224245/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2453688/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:1153206/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2559932/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2743490/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2778623/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2843490/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2999545/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3033781/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3061786/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3149944/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3239094/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3454695/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3459867/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:3547743/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:438400/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:502393/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:740112/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:77102/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:873108/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:952411/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2430721/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:1180186/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:1597059/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:1774304/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:1871948/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:204986/62‑1 (MQ=255)
tAGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2386280/62‑1 (MQ=255)
tAGTGACCTTTCGGGGTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2450599/62‑1 (MQ=255)
 aGTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcg   <  1:2311331/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1571257/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:325125/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3401086/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1693146/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1685914/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3534490/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3539960/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2520140/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3588685/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1137670/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1349794/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:641945/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1226952/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:791520/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1178999/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:891462/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2370815/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2151/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2334283/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2822176/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2882313/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2883194/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2883327/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2961078/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2986403/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2291706/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:1766326/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2594738/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2022526/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3066901/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3076932/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:2590498/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:3195874/61‑1 (MQ=255)
  gTGACCTTTCGGGTTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTAcgg  <  1:322314/61‑1 (MQ=255)
               |                                               
TAGTGACCTTTCGGGCTGGCTTCCAGTACCTGTACCGGCAGCGGCGAATCGAGGCTGGTACGG  >  minE/1592959‑1593021

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: ‑