Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,700,761 C→T 27.6% T123T (ACC→ACT hofQ → predicted fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,700,7610CT27.6% 70.6 / 24.1 65T123T (ACC→ACThofQpredicted fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base C (26/21);  new base T (12/6);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

CGCCAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAG  >  minE/2700699‑2700826
                                                              |                                                                 
cgccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             <  1:2561824/69‑1 (MQ=255)
cgccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             <  1:952244/69‑1 (MQ=255)
cgccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             <  1:2373819/69‑1 (MQ=255)
cgccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             <  1:2837307/69‑1 (MQ=255)
cgccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             >  1:2911970/1‑69 (MQ=255)
 gccAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCaa                                                             <  1:489280/68‑1 (MQ=255)
                              gCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCg                               >  1:1282041/1‑69 (MQ=255)
                              gCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCg                               >  1:1994876/1‑69 (MQ=255)
                               cAAATCTGCCGCTGGTAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:3038090/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:458254/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:536871/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:2104429/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:1859171/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGcc                                                       >  1:1360949/1‑44 (MQ=255)
                               cAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGc                                                        >  1:853098/1‑43 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:2738980/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:460513/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:264604/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:2642503/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:954930/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:952115/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:539940/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:2366843/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:1909212/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:215338/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:2099768/1‑66 (MQ=255)
                                  aTCTGCCGCTGGAAAATCGCAATATTACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCgg                              >  1:1943000/1‑66 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:2297232/68‑1 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:611579/68‑1 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:1712715/68‑1 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:509261/68‑1 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:1749646/68‑1 (MQ=255)
                                         gCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGc                     <  1:677295/68‑1 (MQ=255)
                                            ggAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACt                 <  1:2148270/69‑1 (MQ=255)
                                               aaaTCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTg                <  1:351391/67‑1 (MQ=255)
                                               aaaTCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTg                <  1:2725159/67‑1 (MQ=255)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:1062280/51‑1 (MQ=37)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:229869/51‑1 (MQ=37)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:616095/51‑1 (MQ=37)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:825886/51‑1 (MQ=37)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:3204842/51‑1 (MQ=37)
                                                        tattACTCTGCAATACGCCGACGCCGGAGAGCTGGCGAAAGCGGGGGAGaa                       <  1:967038/51‑1 (MQ=37)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGa               >  1:2754609/1‑57 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:458485/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:870260/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:394649/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2565841/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2286667/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:1064896/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:1112295/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:1414822/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:1822943/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2087989/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2293210/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2440034/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2703990/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:2857918/1‑68 (MQ=255)
                                                          taACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAggg    >  1:3305846/1‑68 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:2469009/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:3354765/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:2500563/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:1998475/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:2517737/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:1164431/68‑1 (MQ=255)
                                                            aCCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAg  <  1:1123861/68‑1 (MQ=255)
                                                              |                                                                 
CGCCAGGAGGCGGAGCAGGCGCGGGCGCAGGCAAATCTGCCGCTGGAAAATCGCAGTATAACCCTGCAATACGCCGACGCGGGAGAACTGGCGAAAGCGGGGGAGAAGCTACTGAGTGCCAAAGGGAG  >  minE/2700699‑2700826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: