Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,001,058 A→G 44.4% F75F (TTT→TTC 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,0580AG44.4% 4.5 / 34.3 45F75F (TTT→TTCgssfused glutathionylspermidine amidase and glutathionylspermidine synthetase
Reads supporting (aligned to +/- strand):  ref base A (14/11);  new base G (7/13);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

TTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACGCATTGCCAC  >  minE/2000999‑2001115
                                                           |                                                         
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCATGTCAGTGAAGACtaca                                                  >  1:2939684/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:862634/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:639275/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:571920/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:1632361/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:3115607/4‑69 (MQ=255)
ttgaCCACTTCACGCAGGAAGCGCAGTGAGAAAATCTCCCATGCCATACCCACGTCAGTGAAGACtaca                                                  >  1:345760/4‑69 (MQ=255)
  aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                                 >  1:1086433/1‑68 (MQ=255)
  aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                                 >  1:1958547/1‑68 (MQ=255)
  aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                                 >  1:3186526/1‑68 (MQ=255)
  aaCCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCacac                                                 >  1:1772547/1‑68 (MQ=255)
            cGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCa                                        <  1:927569/67‑1 (MQ=255)
            cGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCa                                        <  1:2064612/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:30902/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1044477/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:2542025/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:965294/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1779516/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1092644/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1226896/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1278680/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:1828098/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACGTCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:2456589/59‑1 (MQ=255)
                             gaAAATCTCCCATGCCATACCCACATCAGTGAAGACTACACCGTAGTTCAGAAAGAGaa                               <  1:500366/59‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:2075809/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:2591575/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:1257385/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:1471895/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:365229/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:3204295/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:2159627/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:3101310/68‑1 (MQ=255)
                                        aCGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           <  1:2972366/68‑1 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGATCAAATTCAACg           >  1:2479479/1‑60 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCa                    >  1:2717520/1‑51 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCa                    >  1:1101270/1‑51 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           >  1:1683298/1‑60 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACg           >  1:6277/1‑60 (MQ=255)
                                                cccACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACGCATTGCCAc  >  1:1796196/1‑69 (MQ=255)
                                                           |                                                         
TTAACCACTTCTCGCAGGAAGCGCAGCGAGAAAATCTCCCACGCCATACCCACGTCAGTAAAGACCACACCGTAATTCAGAAAGAGAAAACGGCGAGCAAATTCAACGCATTGCCAC  >  minE/2000999‑2001115

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: