Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,042,223 G→A 100% V52V (GTG→GTA gst → glutathionine S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,042,2230GA78.4% 91.0 / 18.6 69V52V (GTG→GTAgstglutathionine S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (3/12);  new base A (0/54);  total (3/66)
Fisher's exact test for biased strand distribution p-value = 8.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACA  >  minE/1042160‑1042293
                                                               |                                                                      
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:2338845/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:2326767/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:2451472/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:2610013/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:859082/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:3186158/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:1429679/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:1339008/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:525753/68‑1 (MQ=255)
tttAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:1386757/68‑1 (MQ=255)
                     aaaCGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACgg                                                 <  1:992901/66‑1 (MQ=255)
                          gTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTg                                                                    <  1:1001344/42‑1 (MQ=255)
                                  tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                 >  1:3015391/1‑69 (MQ=255)
                                  tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                 >  1:2821295/1‑69 (MQ=255)
                                  tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                 >  1:1868905/1‑69 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTTCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2473944/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2890861/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:790690/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:991137/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2346282/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:988046/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:971577/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2917405/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2960691/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2978980/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:817285/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:3228999/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:380911/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:384991/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:813463/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:391452/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:115267/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:1413902/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:1503419/62‑1 (MQ=255)
                                                    aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:1517631/62‑1 (MQ=255)
                                                         acAGGTACCTGCATTGCTGGTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1023603/51‑1 (MQ=255)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1268230/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1368235/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:728098/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:695569/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1317763/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:562262/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1408040/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:517558/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:513714/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:504524/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1283818/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:432579/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:428223/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2216368/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2077293/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2009973/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:922386/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:3242599/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1518327/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:3219301/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:155020/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:3135936/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1082497/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:3038238/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1552937/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1781814/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1078254/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2927028/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:1947082/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2006018/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                           <  1:2825336/51‑1 (MQ=37)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2315891/56‑1 (MQ=38)
                                                         acAGGTACCTGCATTGCTGCTGGACGATGGTAATTTGCTGACGGAAGGCGTa                           <  1:305514/51‑1 (MQ=37)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:569119/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:498208/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:3386595/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:865048/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:3291959/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:1810903/51‑1 (MQ=255)
                                                               gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                      <  1:2890228/51‑1 (MQ=255)
                                                               gCCTGCATTGCTG‑‑GGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACa  >  1:2255588/1‑69 (MQ=255)
                                                               |                                                                      
TTTAATGAAAAAACGTCTCGAAAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACA  >  minE/1042160‑1042293

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: