Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,284,114 T→C 48.0% S1965P (TCA→CCA)  yeeJ → adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,284,1140TC48.0% 2.7 / 114.0 125S1965P (TCA→CCA) yeeJadhesin
Reads supporting (aligned to +/- strand):  ref base T (39/26);  new base C (14/46);  total (53/72)
Fisher's exact test for biased strand distribution p-value = 5.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 125 total aligned reads displayed.

GTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284048‑1284163
                                                                  |                                                 
gTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCa                                                 <  1:1168627/69‑1 (MQ=255)
gTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCa                                                 <  1:1590152/69‑1 (MQ=255)
           tttCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      >  1:1241583/1‑69 (MQ=255)
           gttCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      >  1:3119208/2‑69 (MQ=255)
               cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:2397105/65‑1 (MQ=255)
               cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:2792111/65‑1 (MQ=255)
               cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:1937503/65‑1 (MQ=255)
               cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:3295064/65‑1 (MQ=255)
               cATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:3397578/65‑1 (MQ=255)
                  aaCGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCa                                      <  1:1229057/62‑1 (MQ=255)
                   aCGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCt                               >  1:1815947/1‑68 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1090379/1‑69 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:281650/1‑69 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1889943/1‑69 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:1705178/1‑69 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:426715/1‑69 (MQ=255)
                     ggCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGc                            >  1:841766/1‑69 (MQ=255)
                        gtcaCAATACAGACACAGATAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:2204535/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3195167/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1193529/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3004057/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:2788779/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:852052/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1269158/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1867221/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1834730/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3157958/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1690874/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:3271160/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1615067/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:601545/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:1387002/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:523499/66‑1 (MQ=255)
                        gtcaCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCttt                         <  1:244441/66‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:980217/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2587025/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2451561/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:95081/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:477059/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:3066497/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:309375/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:81608/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:788609/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:3308477/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:522026/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:45563/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1942884/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:2261118/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1515444/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1196447/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1291701/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1993847/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1575628/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATc                      <  1:1652571/69‑1 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:2296253/1‑69 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:1924680/1‑69 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:786642/1‑69 (MQ=255)
                           aCAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATc                      >  1:771484/1‑69 (MQ=255)
                                acagacacagacAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGa                 >  1:5151/1‑69 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTTGCAACTCAAGCACCGCCCATGTTCTAGCttt                         <  1:3044072/59‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:1705196/60‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:3383530/60‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:380671/60‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:1304132/60‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:3099334/60‑1 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:1556537/1‑67 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:502280/1‑67 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:82770/1‑67 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:280502/1‑67 (MQ=255)
                                 cagacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTg                  >  1:2410764/1‑67 (MQ=255)
                                   gacacagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATcc              <  1:801326/69‑1 (MQ=255)
                                       cagacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCttt                         <  1:2203479/54‑1 (MQ=255)
                                         gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:1734703/1‑69 (MQ=255)
                                         gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:959548/1‑69 (MQ=255)
                                         gacaACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGAc        >  1:1432185/1‑69 (MQ=255)
                                         gacaACCGTCAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGa         >  1:2198153/1‑68 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGccc                                       >  1:1993294/1‑36 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:317396/1‑69 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:2048239/1‑69 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:1398825/1‑69 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:2735281/1‑69 (MQ=255)
                                           caACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTa      >  1:2238324/1‑69 (MQ=255)
                                             aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:1026425/69‑1 (MQ=255)
                                             aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2478479/69‑1 (MQ=255)
                                             aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2485235/69‑1 (MQ=255)
                                             aCCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2695576/69‑1 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:790081/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1171693/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:201990/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1301504/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3004898/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1447427/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3096995/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:3171438/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1595098/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATcc              >  1:1386683/1‑58 (MQ=255)
                                              ccGTGAAGGTTACTGGCAACCCAAGCACCGCACATGTTGCTAGCTTTATCGCTGATc               >  1:3370054/1‑57 (MQ=255)
                                              ccGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:2623362/68‑1 (MQ=255)
                                              ccGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATc    <  1:549640/68‑1 (MQ=255)
                                               cGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATcgc  >  1:1545416/1‑69 (MQ=255)
                                                                  |                                                 
GTGACTGCATCTTTCCATAACGGCGTAACAATACAGACACAGACAACCGTGAAAGTCACTGGCAACTCAAGCACCGCCCATGTTGCTAGCTTTATCGCTGATCCATCGACTATCGC  >  minE/1284048‑1284163

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: