Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 395,924 T→A 57.0% V425V (GTT→GTA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE395,9240TA57.0% 26.7 / 74.7 72V425V (GTT→GTAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (23/8);  new base A (24/17);  total (47/25)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCGACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTATTTTCG  >  minE/395863‑395988
                                                             |                                                                
aCGTCGACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTgg                                                           <  1:1216574/69‑1 (MQ=255)
aCGTCGACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTgg                                                           <  1:292407/69‑1 (MQ=255)
aCGTCGACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTgg                                                           <  1:34085/69‑1 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1765000/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:2459522/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1787063/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:2158087/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:2329571/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:263725/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:3136478/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:3252835/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:391265/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:5333/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:587102/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:66765/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:883279/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:981766/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:157335/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1411579/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1354717/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1138839/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1331847/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAatat                                                      >  1:1525377/1‑69 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAata                                                       >  1:1611394/1‑68 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAata                                                       >  1:974442/1‑68 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGCAAAGTCGTAGGTCTGGAata                                                       >  1:3418258/1‑68 (MQ=255)
     gACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGTAGAAAAGTCGTAGGTCTGGAata                                                       >  1:2933372/1‑68 (MQ=255)
        attattCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAata                                                       >  1:502304/1‑65 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:2213818/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:1413265/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:3431500/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:2396844/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:1882467/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:1548073/1‑69 (MQ=255)
           attCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGAt                                                >  1:3061261/1‑69 (MQ=255)
              cTGAATGCGCAAACCACGGAAGTTAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCgt                                             >  1:2774064/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgcg                                      >  1:2300311/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgcg                                      >  1:5787/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgcg                                      >  1:3135382/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgcg                                      >  1:1431552/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgcg                                      >  1:298342/1‑69 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAgc                                       >  1:813544/1‑68 (MQ=255)
                     cgcAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTAAgcg                                      >  1:3389951/1‑69 (MQ=255)
                                      aaaGGCGACGGCAGCAAAGTCGTTGGTCTGGATTATCGAGATCGTGTCAGCGGCGATATTcac                           <  1:2133979/63‑1 (MQ=255)
                                      aaaGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTcac                           <  1:1062452/63‑1 (MQ=255)
                                      aaaGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTcac                           <  1:2634941/63‑1 (MQ=255)
                                      aaaGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATAGTGTCAGCGGCGATATTCAAACATCg                    >  1:508857/1‑69 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATct                    <  1:287758/68‑2 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:1493446/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:1758742/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:1873356/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:1103461/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2075539/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:1110828/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2130336/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2180873/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2295345/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:353815/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2406819/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:3340196/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:3311259/68‑1 (MQ=255)
                                        aGGCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:260662/68‑1 (MQ=255)
                                        aGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCg                    <  1:1319405/68‑1 (MQ=255)
                                         ggCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2310670/67‑1 (MQ=255)
                                         ggCGACGGTAGCAAAGTCGTAGGTCTGGAATATCGCGATCGTGTCAGCGGCAATATTCACAACATCg                    <  1:2688288/67‑1 (MQ=255)
                                                  agcaAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGt         >  1:3265869/1‑69 (MQ=255)
                                                  agcaAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGt         <  1:3044397/69‑1 (MQ=255)
                                                     aaaGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTAtt      >  1:2168177/1‑69 (MQ=255)
                                                     aaaGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTCGCCGGTa        >  1:338217/1‑67 (MQ=255)
                                                         tCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTATTTTCg  >  1:1039387/1‑69 (MQ=255)
                                                         tCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTATTTTCg  >  1:1398758/1‑69 (MQ=255)
                                                         tCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTATTTTCg  >  1:3243855/1‑69 (MQ=255)
                                                             |                                                                
ACGTCGACATTATTCTGAATGCGCAAACCACGGAAGTGAAAGGCGACGGCAGCAAAGTCGTTGGTCTGGAATATCGAGATCGTGTCAGCGGCGATATTCACAACATCGAACTGGCCGGTATTTTCG  >  minE/395863‑395988

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: