Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,322,829 C→A 48.4% T81T (ACG→ACT sbp ← sulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,322,8290CA48.4% 11.5 / 65.2 62T81T (ACG→ACTsbpsulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (24/8);  new base A (23/7);  total (47/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TGATCCACTCTTTATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGTTTAC  >  minE/2322765‑2322878
                                                                |                                                 
tGATCCACTCTTTATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTg                                                >  1:798946/1‑68 (MQ=255)
tGATCCACTCTTTATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTg                                                >  1:3267302/1‑68 (MQ=255)
            tATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGc                                     <  1:2818690/67‑1 (MQ=255)
            tATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGc                                     <  1:1435965/67‑1 (MQ=255)
            tATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGc                                     <  1:3089080/67‑1 (MQ=255)
            tATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGc                                     <  1:2019675/67‑1 (MQ=255)
            tATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGc                                     <  1:326606/67‑1 (MQ=255)
              tCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTc                                  <  1:1032863/68‑1 (MQ=255)
              tCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTc                                  <  1:2772216/68‑1 (MQ=255)
              tCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTc                                  <  1:3399911/68‑1 (MQ=255)
                   ccgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATAcc                            >  1:580037/1‑69 (MQ=255)
                   ccgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATAcc                            >  1:710228/1‑69 (MQ=255)
                   ccgcccacGTTCCGCAATTGCGTCCACGTGATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATa                              >  1:441602/1‑67 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGac                                              >  1:1833063/1‑50 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATAcc                            >  1:1913093/1‑68 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATAcc                            >  1:2712079/1‑68 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:1715756/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:1790417/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:3247566/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:2503655/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:2445433/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:827308/1‑69 (MQ=255)
                    cgcccgcGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:2818122/1‑69 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:2671243/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:1185929/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:2901048/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:3014054/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:3209375/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:3399135/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:249691/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:160008/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:2184289/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:1745341/1‑68 (MQ=255)
                        cgcgTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTg                        >  1:2579097/1‑68 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCATCGTGACAACATCAGCTTCAATACCGTTGATTACCGAct              >  1:2507548/1‑64 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATa                              >  1:1519090/1‑49 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCg                           >  1:1691849/1‑52 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:1986280/1‑69 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:824187/1‑69 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:820226/1‑69 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:550429/1‑69 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:424264/1‑69 (MQ=255)
                                     tGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGc          >  1:301499/1‑69 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATAcc                            >  1:1835929/1‑48 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:3173552/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:1374696/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:2336157/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:2784554/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:2798814/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:2297965/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:163649/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:1904863/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:3212681/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:2008795/1‑49 (MQ=255)
                                        gTCCACGTCATAGGCCAGAGCCAGAGTGACAACATCAGCTTCAATACCg                           >  1:3102831/1‑49 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:3083311/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:2297567/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:2023646/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:732093/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:1428951/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:861335/62‑3 (MQ=255)
                                                    ggCCAGAGCCAGAGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:947590/62‑3 (MQ=255)
                                                     gccagagccagaGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGtttgc  <  1:2168572/61‑3 (MQ=255)
                                                                |                                                 
TGATCCACTCTTTATCAATCCGCCCGCGTTCCGCAATTGCGTCCACGTCATAGGCCAGAGCCAGCGTGACAACATCAGCTTCAATACCGTTGATTACCGACGTCGCTTGTTTAC  >  minE/2322765‑2322878

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: