Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,001,040 C→T 25.8% A81A (GCG→GCA 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,0400CT25.8% 91.9 / 26.4 77A81A (GCG→GCAgssfused glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (39/18);  new base T (7/13);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 1.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACG  >  minE/2000974‑2001106
                                                                  |                                                                  
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:2331084/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:44620/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:1161286/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:160150/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:681924/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:2587500/69‑1 (MQ=255)
tGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  <  1:1530490/69‑1 (MQ=255)
           gCAGGATGTTGTCATTTACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                                                         >  1:3069447/1‑67 (MQ=255)
           gCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                                                         >  1:312684/1‑67 (MQ=255)
           gCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCAc                                                         >  1:3118337/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:91481/1‑53 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:911761/1‑53 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:857385/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1028418/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:608032/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:3191853/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:350201/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:381694/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1679859/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1641074/1‑67 (MQ=255)
                aTGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1554504/1‑67 (MQ=255)
                 tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1305184/1‑66 (MQ=255)
                 tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:1346954/1‑66 (MQ=255)
                 tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:2960428/1‑66 (MQ=255)
                 tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:282534/1‑66 (MQ=255)
                 tgttgtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAg                                                    >  1:2563893/1‑66 (MQ=255)
                     gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:881024/1‑48 (MQ=255)
                     gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:344468/1‑48 (MQ=255)
                     gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:2648724/1‑48 (MQ=255)
                     gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:670023/1‑48 (MQ=255)
                     gtCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:361418/1‑48 (MQ=255)
                         ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCATGTCAGTGAAGACtaca                                         >  1:2939684/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)
                         ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                         >  1:639275/4‑69 (MQ=255)
                         ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                         >  1:571920/4‑69 (MQ=255)
                         ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                         >  1:345760/4‑69 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAACTCTCCCACGc                                                                  >  1:1447918/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:2235055/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:325304/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:964126/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:818617/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGc                                                                  >  1:215726/1‑42 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                        >  1:1958547/1‑68 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                        >  1:1772547/1‑68 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                        >  1:3186526/1‑68 (MQ=255)
                           aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                        >  1:1086433/1‑68 (MQ=255)
                              cACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATAcc                                                            >  1:3257268/1‑45 (MQ=255)
                                     cGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCa                               <  1:2064612/67‑1 (MQ=255)
                                     cGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCa                               <  1:927569/67‑1 (MQ=255)
                                          gAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAaga                         >  1:2506944/1‑68 (MQ=255)
                                           aaGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAagag                        >  1:1721002/1‑68 (MQ=255)
                                           aaGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAagag                        >  1:13784/1‑68 (MQ=255)
                                           aaGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAagag                        >  1:2108029/1‑68 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:497315/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:687833/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:965294/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1044477/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1092644/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:30902/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1226896/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:2542025/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:2456589/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1278680/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1828098/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:1779516/59‑1 (MQ=255)
                                                      gaAAATCTCCCATGCCATACCCACATCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                      <  1:500366/59‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:2159627/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:1257385/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:1471895/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:2075809/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:365229/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:2591575/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:2972366/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:3101310/68‑1 (MQ=255)
                                                                 aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg  <  1:3204295/68‑1 (MQ=255)
                                                                  |                                                                  
TGCCTGCAATGGCAGGATGTTGTCATTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACG  >  minE/2000974‑2001106

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: