Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 540,307 C→G 45.4% R236R (CGG→CGC rhtA ← threonine and homoserine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE540,3070CG45.4% 24.1 / 110.1 97R236R (CGG→CGCrhtAthreonine and homoserine efflux system
Reads supporting (aligned to +/- strand):  ref base C (39/14);  new base G (25/19);  total (64/33)
Fisher's exact test for biased strand distribution p-value = 9.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TCTCCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGC  >  minE/540239‑540372
                                                                    |                                                                 
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:3183349/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:2459946/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:2506693/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:2319644/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:2201763/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:2939354/1‑69 (MQ=255)
tctcCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTc                                                                   >  1:181514/1‑69 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:1170412/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:1291620/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:1667450/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:1346343/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:2612658/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:2789573/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:628860/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:494685/66‑1 (MQ=255)
      aGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCgt                                                                <  1:2233307/66‑1 (MQ=255)
                           gCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATc                                        >  1:3149018/1‑69 (MQ=255)
                           gCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATc                                        >  1:2566441/1‑69 (MQ=255)
                              gccagcGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     <  1:1068348/69‑1 (MQ=255)
                              gccagcGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     <  1:2250387/69‑1 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAACCGGGTGAGGGCAATCAtt                                     >  1:2146239/1‑69 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     >  1:1482051/1‑69 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     >  1:1334082/1‑69 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     >  1:2455056/1‑69 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     >  1:1994301/1‑69 (MQ=255)
                              gccagcGCAGGTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCAtt                                     >  1:1234603/1‑69 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:211344/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:2072244/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1987834/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:30845/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1213776/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:295168/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:751314/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1287264/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:654301/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1426988/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:391883/67‑1 (MQ=255)
                                     caggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1885912/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:57135/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:2449199/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1977010/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:882021/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:379833/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:3312824/67‑1 (MQ=255)
                                      aggTTCCATGCTCATTAGCGTACCAAATGTGCGGGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCg                              <  1:1400508/67‑1 (MQ=255)
                                          tCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATaa                         >  1:2027307/1‑69 (MQ=255)
                                          tCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATaa                         >  1:108683/1‑69 (MQ=255)
                                          tCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATaa                         >  1:2063448/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGTGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:718269/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:3289984/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:2802642/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:430828/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:957155/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:3398254/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:124891/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:1798366/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:1775655/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCag                    >  1:1619604/1‑69 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCa                     >  1:1718663/1‑68 (MQ=255)
                                               gCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTACATCGAATAAGGCag                    >  1:3110876/1‑69 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGTTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:3232172/2‑67 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGGAGAGCgg               >  1:1456927/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:492134/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1639178/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:429614/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1686440/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1506725/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:3420905/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1290482/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:3160278/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1808574/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:1007033/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:2528592/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:2893479/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:2834536/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCgg               >  1:2817443/2‑68 (MQ=255)
                                                     taGCGTACCAAATGTGCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGAgg               >  1:2977861/2‑68 (MQ=255)
                                                      aGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTTCAGCGAATAAGGCAGAGCgg               >  1:2206709/1‑67 (MQ=255)
                                                      aGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCa                     >  1:521656/1‑61 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:1840893/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:3159727/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:1414518/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:63605/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:1734252/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:2348202/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:2846025/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:2800753/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGCCAGAGCGGTCgagag        >  1:3151475/1‑68 (MQ=255)
                                                            ccAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATATCCAGCGAATAAGGCAGAGCGGTCgagag        >  1:2217681/1‑68 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGGAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:3403174/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:2628623/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                <  1:2396636/55‑1 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:1942440/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:1903985/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:114249/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:907758/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:923800/1‑55 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                <  1:2409625/55‑1 (MQ=255)
                                                                 tGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGc  >  1:2539127/1‑69 (MQ=255)
                                                                 tGTCCGTGTTGGAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCg                >  1:1766181/1‑54 (MQ=255)
                                                                    |                                                                 
TCTCCGAGGAAAATCATCCCGGAAACGGCAGCCAGCGCCGGTTCCATGCTCATCAGCGTACCAAATGTCCGTGTTGGCAAACGGGTGAGGGCAATCATTTCCAGCGAATAAGGCAGAGCGGTCGAGAGAATAGC  >  minE/540239‑540372

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: