Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,086,684 C→A 57.7% I153I (ATC→ATA yqjG → predicted S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,086,6840CA57.7% 19.8 / 80.9 78I153I (ATC→ATAyqjGpredicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (22/11);  new base A (16/29);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTAC  >  minE/2086619‑2086735
                                                                 |                                                   
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATcc                                                    >  1:3041651/1‑67 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:907513/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:3321021/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:3062880/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:2845762/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:2528050/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:2504710/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:1609274/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:2155104/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:2000943/1‑68 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:2514760/1‑69 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:885652/1‑69 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:3238642/1‑69 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:3070591/1‑69 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:2944467/1‑69 (MQ=255)
aCTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGc                                                  >  1:2270915/1‑69 (MQ=255)
aCTGTTCCAGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCg                                                   >  1:1328167/1‑68 (MQ=255)
      cccGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCa                                                 >  1:693793/1‑64 (MQ=255)
      cccGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            >  1:2363825/1‑69 (MQ=255)
      cccGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            >  1:3324156/1‑69 (MQ=255)
            cTGTGGGACAAAAAGAACCTCACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:862260/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:1005849/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:2691028/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:1171316/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:3094330/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:1555967/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:3169412/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:2214862/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:2449140/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:2243470/63‑1 (MQ=255)
            cTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTc                                            <  1:2303712/63‑1 (MQ=255)
                 ggACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGt                                  >  1:2904749/1‑68 (MQ=255)
                 ggACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGt                                  >  1:3415036/1‑68 (MQ=255)
                 ggACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCAACCGCATg                                               >  1:2225748/1‑55 (MQ=255)
                                   cATCGTCAGCAACGAATCATCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1617730/69‑1 (MQ=39)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTTGGTGCGaa               <  1:3262766/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1056881/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:3332760/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:3281553/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:3107966/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:526156/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:542721/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:593182/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:767026/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1309470/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1032548/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2822952/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1876440/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2073722/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2139288/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1664808/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1646979/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:222357/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1836846/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2486261/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1606973/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1005549/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:1927085/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2634181/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2758427/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:27937/69‑1 (MQ=255)
                                   cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCCCGTTTGATGCGCTGGGTGCGaa               <  1:2772323/69‑1 (MQ=255)
                                    aTCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa               <  1:2317907/68‑1 (MQ=255)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:1958383/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:872820/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:752555/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:71843/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:1878711/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:619579/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:2228339/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:2093880/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:429786/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:3295588/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:2449748/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:1124460/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:2573056/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:3092583/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:1469635/1‑64 (MQ=38)
                                                     agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac  >  1:153232/1‑64 (MQ=38)
                                                                 |                                                   
ACTGTTCCCGTGCTGTGGGACAAAAAGAACCACACCATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTAC  >  minE/2086619‑2086735

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: