Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,121,626 T→C 58.0% M59V (ATG→GTG)  btuE ← predicted glutathione peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,121,6260TC58.0% 21.9 / 54.2 76M59V (ATG→GTG) btuEpredicted glutathione peroxidase
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (23/21);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACA  >  minE/1121562‑1121692
                                                                |                                                                  
ttAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATaaaa                                                                >  1:1028110/1‑69 (MQ=255)
     ctctTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATTCCAGCACCATAAAACCtc                                                            >  1:2263581/1‑68 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:798282/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2278240/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:427827/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2216503/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2127998/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:806974/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2536025/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:3174739/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2941768/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2571501/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:2855923/66‑2 (MQ=255)
            tCGCTACCCGGTTCTTGCTCCAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCacgatca                                                       <  1:1824498/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:2003058/1‑68 (MQ=255)
                               ccAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATtc                                  >  1:899645/1‑68 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:929009/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:2272876/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:541073/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:385883/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:2362070/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:206436/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:3347877/68‑1 (MQ=255)
                               ccAGAAACTGGTTACACGGGAATCCCAGCACCACAAAACCACGATCAGCCCAGGCTTTCTGAATATtc                                  <  1:954514/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:1130251/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:469848/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1292487/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:628020/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1296920/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1321929/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:630579/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3411250/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3395973/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:630646/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3274308/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:3161735/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1582942/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2980290/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2963872/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:171521/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1835752/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:1859232/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2707523/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2576707/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2022561/1‑37 (MQ=21)
                                                  gAATCCCAGCACCACAAAACCACGATCAGCCCAGGCt                                              >  1:2100877/1‑37 (MQ=21)
                                                       ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:1735397/39‑1 (MQ=255)
                                                       ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:2551721/39‑1 (MQ=255)
                                                       ccAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGaa                                       <  1:2347281/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)
                                                                |                                                                  
TTAATCTCTTCATCGCTGCCCGGTTCTTGTTCCAGAAACTGGTTGCACGGGAATCCCAGCACCATAAAACCTCGATCGACCCAGGCTTTCTGAATATTCTCCAACTGCTCATATTGCGGCGTTAAGCCACA  >  minE/1121562‑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: