Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,042,241 T→C 100% D58D (GAT→GAC gst → glutathionine S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,042,2410TC79.0% 121.9 / 21.8 85D58D (GAT→GACgstglutathionine S‑transferase
Reads supporting (aligned to +/- strand):  ref base T (10/8);  new base C (0/67);  total (10/75)
Fisher's exact test for biased strand distribution p-value = 1.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCC  >  minE/1042181‑1042308
                                                            |                                                                   
aaaCGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACgg                                                                <  1:992901/66‑1 (MQ=255)
             tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                                >  1:2821295/1‑69 (MQ=255)
             tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                                >  1:3015391/1‑69 (MQ=255)
             tACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGaa                                                >  1:1868905/1‑69 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTTCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:2473944/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:115267/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:2346282/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:2890861/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:991137/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:988046/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:391452/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:971577/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:1503419/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:1413902/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:817285/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:813463/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:790690/62‑1 (MQ=255)
                               aaGGGACAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:1517631/62‑1 (MQ=255)
                                    acAGGTACCTGCATTGCTGGTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1023603/51‑1 (MQ=255)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:3135936/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2825336/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:922386/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2927028/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:728098/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:695569/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:3038238/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:3219301/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:562262/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:517558/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:3242599/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:513714/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:428223/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:432579/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:504524/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2216368/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1947082/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2006018/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1317763/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1368235/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2009973/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1781814/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:2077293/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1283818/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1268230/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1552937/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:155020/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1408040/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1518327/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1082497/51‑1 (MQ=37)
                                    acAGGTACCTGCATTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTa                                          <  1:1078254/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:1810903/51‑1 (MQ=255)
                                          gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:3291959/51‑1 (MQ=255)
                                          gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:3386595/51‑1 (MQ=255)
                                          gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:2890228/51‑1 (MQ=255)
                                          gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:865048/51‑1 (MQ=255)
                                          gCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:498208/51‑1 (MQ=255)
                                          gCCTGCATTGCTG‑‑GGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACa                 >  1:2255588/1‑69 (MQ=255)
                                          accTGCTTTGGTGCTGGACGATGGTAGTTTGCTGACGGATGGCGTa                                          <  1:835215/45‑1 (MQ=255)
                                                aTTGCTGCTGGACGATGGTACTTTGCTGACGGAAGGCGTAGCGAt                                     <  1:1306576/45‑1 (MQ=38)
                                                     tgcttGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:1048962/68‑1 (MQ=255)
                                                     tgctgGACGTTGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:2964935/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCTTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:2807387/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:1506533/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:1542816/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:3253106/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:3124591/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:2774130/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:2555962/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:2668502/68‑1 (MQ=255)
                                                     tgctgGACGATGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTcccc         <  1:887652/68‑1 (MQ=255)
                                                         gggtGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:1541795/4‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGccc     >  1:2057397/1‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:2737030/1‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:1080375/1‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:967179/1‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:1263167/1‑68 (MQ=255)
                                                         ggATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGAcc     >  1:2248926/1‑68 (MQ=255)
                                                            tGACGGTACTTTGCTGACGGAAGGCTTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGcc  <  1:1830870/68‑1 (MQ=255)
                                                            tGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGcc  <  1:1193685/68‑1 (MQ=255)
                                                            |                                                                   
AAACGGTGACGATTACTTTGCCGTTAACCCTAAGGGGCAGGTGCCTGCATTGCTGCTGGATGACGGTACTTTGCTGACGGAAGGCGTAGCGATTATGCAGTATCTTGCCGACAGCGTCCCCGACCGCC  >  minE/1042181‑1042308

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: