Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,001,010 T→A 16.2% R91R (CGA→CGT gss ← fused glutathionylspermidine amidase and glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,001,0100TA16.2% 173.8 / 23.0 86R91R (CGA→CGTgssfused glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base T (49/23);  new base A (7/7);  total (56/30)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACAC  >  minE/2000954‑2001068
                                                        |                                                          
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTGTCGCAGGAAgcg                                                 >  1:2145754/1‑68 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:1066961/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:1794406/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:1813201/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:2075934/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:1524404/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:2315887/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:2533970/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:1264780/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:3272174/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:334637/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:772987/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:804088/1‑69 (MQ=255)
cgcgGCGAGCCGTTAGGAAATGACTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAgcgc                                                >  1:2254072/1‑69 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:2722276/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:2866742/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:2265844/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:989634/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:1666007/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:1469797/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTGACCACTTCACGCAgg                                                      <  1:121291/54‑1 (MQ=38)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:1549536/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:250950/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:2500718/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:241190/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:1206742/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:273266/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:2797728/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:2084541/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:2084401/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:1797533/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:935057/69‑1 (MQ=255)
         ccGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGaaa                                       <  1:942400/69‑1 (MQ=255)
           gTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAAt                                     <  1:364635/69‑1 (MQ=255)
           gTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAAt                                     <  1:2430722/69‑1 (MQ=255)
           gTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAAt                                     <  1:1182599/69‑1 (MQ=255)
                gAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTcc                                 <  1:928914/68‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:2331084/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:44620/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:2587500/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:681924/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:1161286/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:1530490/69‑1 (MQ=255)
                    tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            <  1:160150/69‑1 (MQ=255)
                               gCAGGATGTTGTCATTTACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                   >  1:3069447/1‑67 (MQ=255)
                               gCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                   >  1:3118337/1‑67 (MQ=255)
                               gCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                   >  1:312684/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:911761/1‑53 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:91481/1‑53 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:3191853/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1679859/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:350201/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:381694/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:608032/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:857385/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1028418/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1641074/1‑67 (MQ=255)
                                    aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1554504/1‑67 (MQ=255)
                                     tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:2960428/1‑66 (MQ=255)
                                     tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1346954/1‑66 (MQ=255)
                                     tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:2563893/1‑66 (MQ=255)
                                     tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:1305184/1‑66 (MQ=255)
                                     tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg              >  1:282534/1‑66 (MQ=255)
                                         gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:670023/1‑48 (MQ=255)
                                         gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:881024/1‑48 (MQ=255)
                                         gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:2648724/1‑48 (MQ=255)
                                         gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:361418/1‑48 (MQ=255)
                                         gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:344468/1‑48 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCATGTCAGTGAAGACtaca   >  1:2939684/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:571920/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:345760/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:639275/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:1632361/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:862634/4‑69 (MQ=255)
                                             ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca   >  1:3115607/4‑69 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAACTCTCCCACGc                            >  1:1447918/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:325304/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:964126/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:215726/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:2235055/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                            >  1:818617/1‑42 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac  >  1:1958547/1‑68 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac  >  1:1772547/1‑68 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac  >  1:1086433/1‑68 (MQ=255)
                                               aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac  >  1:3186526/1‑68 (MQ=255)
                                                  cACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATAcc                      >  1:3257268/1‑45 (MQ=255)
                                                        |                                                          
CGCGGCGAGCCGTTAGGAAATGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACAC  >  minE/2000954‑2001068

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