Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,000,557 A→G 15.9% G242G (GGT→GGC gss ← fused glutathionylspermidine amidase and glutathionylspermidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,000,5570AG15.9% 86.7 / 7.8 56G242G (GGT→GGCgssfused glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/21);  new base G (2/7);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

TCCTGCTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTG  >  minE/2000498‑2000622
                                                           |                                                                 
tCCTGCTCGGCACTCTCGGTAATGGTGTAGTAATGATACGGATCCTGATTTATCACCTGGCCGTTGGcc                                                          >  1:2062418/1‑69 (MQ=255)
tCCTGCTCGGCACTCTCGGTAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGcc                                                          >  1:3044948/1‑69 (MQ=255)
    gCTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGca                                                      <  1:1080208/69‑1 (MQ=255)
    gCTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGca                                                      <  1:1246745/69‑1 (MQ=255)
     cTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGca                                                      >  1:2196367/1‑68 (MQ=255)
     cTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGca                                                      >  1:1179058/1‑68 (MQ=255)
     cTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGATGGCCTGc                                                       >  1:489217/1‑67 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:1038797/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:2546533/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:2852490/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:288834/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:312384/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTTGGCCTGTACATAGGCGtttt                                          <  1:975532/66‑1 (MQ=255)
                   tAATGGTGTAGTAATGATACGGATCCTGATTGATCACCTGGCCGTAGGCCTGTACATAGGCGtttt                                          <  1:2862550/66‑1 (MQ=255)
                             gtaATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCttt                             >  1:1825523/1‑69 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGca                                                      >  1:2677683/1‑41 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:1351246/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:906145/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:813748/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:1822045/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:1575217/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCg                                              >  1:1282506/1‑49 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTc                          >  1:2852803/1‑69 (MQ=255)
                                aTGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTc                          >  1:1446215/1‑69 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCtttttaa                         >  1:683500/1‑65 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:732081/1‑67 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:3001311/1‑67 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:2513511/1‑67 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:2509803/1‑67 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:1227029/1‑67 (MQ=255)
                                   aTAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCa                         >  1:1595469/1‑67 (MQ=255)
                                      agtaTCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:3396053/66‑1 (MQ=255)
                                      aGGATCGTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:363397/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:3154529/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:1075226/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:111915/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:3407300/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:2277631/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:3028231/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:225442/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:1881469/69‑1 (MQ=255)
                                      aGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAg                    <  1:2754017/69‑1 (MQ=255)
                                        gATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGcc                  <  1:1107491/69‑1 (MQ=255)
                                            cTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCAttt              <  1:1687092/69‑1 (MQ=255)
                                            cTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCAttt              <  1:810731/69‑1 (MQ=255)
                                            cTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCAttt              <  1:1223149/69‑1 (MQ=255)
                                                    tCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCa       <  1:2199646/68‑1 (MQ=255)
                                                    tCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCa       <  1:2821665/68‑1 (MQ=255)
                                                    tCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCa       <  1:780320/68‑1 (MQ=255)
                                                    tCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCa       <  1:58950/68‑1 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACt   >  1:134082/1‑68 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTg  >  1:164765/1‑69 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTg  >  1:2019688/1‑69 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTg  >  1:3221998/1‑69 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTg  >  1:1519579/1‑69 (MQ=255)
                                                        cTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTg  >  1:1526892/1‑69 (MQ=255)
                                                           |                                                                 
TCCTGCTCGGCACTCTCGGTAATGGTGTAGTAATGATAAGGATCCTGATTGATCACCTGACCGTTGGCCTGCACATAGGCGTTTTGCAGCGGATCTTTTTCATCCAGCCATTTACCGTCAAACTG  >  minE/2000498‑2000622

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