Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,121,639 G→A 47.7% V54V (GTC→GTT) ‡ btuE ← predicted glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,121,6390GA47.7% 8.4 / 54.4 65V54V (GTC→GTT) ‡btuEpredicted glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base G (14/20);  new base A (23/8);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

GCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATT  >  minE/1121576‑1121707
                                                               |                                                                    
tcTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                 >  1:2783808/2‑69 (MQ=255)
gCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                 >  1:3428475/1‑69 (MQ=255)
gCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                 >  1:2333210/1‑69 (MQ=255)
gCTGCCCGGTTCTTGTTCCAGAAACTGGTTCCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                                 >  1:631572/1‑69 (MQ=255)
                 ccAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                 >  1:2003058/1‑68 (MQ=255)
                 ccAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                                 >  1:899645/1‑68 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:3347877/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:206436/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:954514/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:929009/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:385883/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:2272876/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:2362070/68‑1 (MQ=255)
                 ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                                 <  1:541073/68‑1 (MQ=255)
                          ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                       <  1:2891951/69‑1 (MQ=255)
                          ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                                       <  1:2861827/69‑1 (MQ=255)
                          ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGATTTCTGAATATTCTCCAACTGCt                                       <  1:2799681/69‑1 (MQ=255)
                             tGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                     <  1:585735/68‑1 (MQ=255)
                             tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                                    <  1:1020093/69‑1 (MQ=255)
                             tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                                    <  1:547664/69‑1 (MQ=255)
                             tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                     <  1:1902983/68‑1 (MQ=255)
                             tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                     <  1:1924596/68‑1 (MQ=255)
                              gCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCa                                     <  1:1784506/67‑1 (MQ=255)
                              gCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCGTTCTGAATATTCTCCAACTGCTCa                                     <  1:1014070/67‑1 (MQ=255)
                                   ggAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTgcg                              >  1:1051120/1‑69 (MQ=255)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGc                                                              >  1:2183274/1‑36 (MQ=19)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2576707/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:3395973/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:3411250/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:3274308/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:3161735/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2980290/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2963872/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1130251/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1292487/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1296920/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1321929/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2707523/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1582942/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2022561/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:171521/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:469848/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:630646/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:630579/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:2100877/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1835752/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:1859232/1‑37 (MQ=21)
                                    gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                                             >  1:628020/1‑37 (MQ=21)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                      <  1:2551721/39‑1 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                      <  1:2347281/39‑1 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                                      <  1:1735397/39‑1 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTcc                                              >  1:432808/1‑47 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCata                                   >  1:2162895/1‑58 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                        >  1:979585/1‑69 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                        >  1:2370488/1‑69 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                        >  1:2322759/1‑69 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                        >  1:1370958/1‑69 (MQ=255)
                                         ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa                        >  1:1759041/1‑69 (MQ=255)
                                                cATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCcaca                 <  1:2055692/69‑1 (MQ=255)
                                                cATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCcaca                 <  1:2035831/69‑1 (MQ=255)
                                                cATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCcaca                 <  1:1562107/69‑1 (MQ=255)
                                                    aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                   <  1:1014727/63‑1 (MQ=255)
                                                    aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                   <  1:1022195/63‑1 (MQ=255)
                                                    aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                   <  1:2598436/63‑1 (MQ=255)
                                                    aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca                   <  1:1774602/63‑1 (MQ=255)
                                                              cgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa    >  1:3353886/1‑68 (MQ=255)
                                                              cgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa    >  1:784032/1‑68 (MQ=255)
                                                              cgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa    >  1:2310166/1‑68 (MQ=255)
                                                              cgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACa    >  1:377734/1‑68 (MQ=255)
                                                              cgaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGATAAGCCACACt               >  1:1629482/1‑57 (MQ=255)
                                                               gaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACAtt  <  1:1928885/69‑1 (MQ=255)
                                                               gaCCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACAtt  <  1:3372894/69‑1 (MQ=255)
                                                               |                                                                    
GCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACACTTTGAGGCGACATT  >  minE/1121576‑1121707

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