Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,121,633 T→A 47.7% R56R (CGA→CGT btuE ← predicted glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,121,6330TA47.7% 16.5 / 75.4 65R56R (CGA→CGTbtuEpredicted glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (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 = 4.13e-01

CTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACA  >  minE/1121567‑1121692
                                                                  |                                                           
ctctTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATTCCAGCACCATAAAACCtc                                                            >  1:2263581/1‑68 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2216503/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:427827/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2536025/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2571501/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2127998/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:3174739/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2941768/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2278240/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:1824498/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2855923/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:806974/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:798282/66‑2 (MQ=255)
       tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGCGAATCCCAGCACCACAAAACCacgatca                                                       <  1:1839766/66‑2 (MQ=25)
         tcTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                  >  1:2783808/2‑69 (MQ=255)
         gCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                  >  1:2333210/1‑69 (MQ=255)
         gCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                  >  1:3428475/1‑69 (MQ=255)
         gCTGCCCGGTTCTTGTTCCAGAAACTGGTTCCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCa                                                  >  1:631572/1‑69 (MQ=255)
                          ccAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                  >  1:899645/1‑68 (MQ=255)
                          ccAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                  >  1:2003058/1‑68 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:541073/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:206436/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:2362070/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:2272876/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:929009/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:3347877/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:385883/68‑1 (MQ=255)
                          ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:954514/68‑1 (MQ=255)
                                   ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                        <  1:2861827/69‑1 (MQ=255)
                                   ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCt                        <  1:2891951/69‑1 (MQ=255)
                                   ggTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGATTTCTGAATATTCTCCAACTGCt                        <  1:2799681/69‑1 (MQ=255)
                                      tGCACGGGAATCCCAGCACCATAAAACCTCGATCGCCCCAGGCTTTCTGAATATTCTCCAACTGCTCa                      <  1:585735/68‑1 (MQ=255)
                                      tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                     <  1:547664/69‑1 (MQ=255)
                                      tGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCat                     <  1:1020093/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:2022561/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:469848/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2100877/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:630579/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1130251/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3411250/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3395973/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1292487/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3274308/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1296920/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:1321929/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1582942/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:171521/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:630646/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1835752/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2707523/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1859232/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2576707/1‑37 (MQ=21)
                                             gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:628020/1‑37 (MQ=21)
                                                  ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:1735397/39‑1 (MQ=255)
                                                  ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:2347281/39‑1 (MQ=255)
                                                  ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:2551721/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:1759041/1‑69 (MQ=255)
                                                  ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTa         >  1:1370958/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:2598436/63‑1 (MQ=255)
                                                             aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca    <  1:1774602/63‑1 (MQ=255)
                                                             aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca    <  1:1022195/63‑1 (MQ=255)
                                                             aaaCCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCca    <  1:1014727/63‑1 (MQ=255)
                                                                  |                                                           
CTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACA  >  minE/1121567‑1121692

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: