Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 157,994 G→T 22.7% Q93K (CAA→AAA)  htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606157,9940G→T22.7% 5.7 ‑4.1 0.026277htrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  new base (6/0):  ref base (8/12):  total (14/12)
Fisher's exact test strand distribution p-value = 1.71e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTA  >  REL606/157959‑158028
                                   |                                  
aaaCTGCAATAAATTCTTAAGTGTGATACAGGCTtg                                    <  1:5710862/36‑1 (MQ=255)
 aaCTGCAATAAATTCTTAAGTGTGATACAGGCTtgt                                   <  1:5968288/36‑1 (MQ=255)
  aCTGCAATAAATTCTTAAGTGTGATACAGGCTtgtg                                  <  1:4462818/36‑1 (MQ=255)
     gCAATAAATTCTTAATTGTGATACAGGCTTTTGCGt                               >  1:685664/1‑36 (MQ=37)
     gCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGt                               <  1:4652015/36‑1 (MQ=255)
      cAATAAATTCTTAAGTGTGATACAGGCTTTTTCGtt                              >  1:7223230/1‑36 (MQ=37)
       aaTAAATTCTTAAGTGTGATACAGGCTTTTTCGttt                             >  1:7168403/1‑36 (MQ=37)
       aaTAAATTCTTAAGTGTGATACAGGCTTTTTCGttt                             >  1:3410256/1‑36 (MQ=37)
          aaaTTCTTAAGTGTGATACAGGCTTGTGCGtttttt                          >  1:6934877/1‑36 (MQ=255)
               cTTAAGTGTGATACAGGCTTTTTTGttttttttttt                     >  1:4872139/1‑34 (MQ=21)
               cTTAAGTGTGATACAGGCTTTTGGGttttttttttt                     >  1:2357676/1‑34 (MQ=37)
                ttAAGTGTGATACAGGCTTGTTTGTTTTTTTTTctt                    >  1:6899796/1‑34 (MQ=37)
                ttAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCt                    <  1:5813884/36‑1 (MQ=255)
                 tAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCtt                   <  1:1601672/36‑1 (MQ=255)
                 tAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCtt                   <  1:2061749/36‑1 (MQ=255)
                  aaGTGTGATACAGGCTTGTGCGTTTTTTTTTCCttc                  <  1:3758263/36‑1 (MQ=255)
                    gtgtGATACAGGCTTGTGCGTTTTTTTTTCCttcaa                <  1:1546814/36‑1 (MQ=255)
                      gtgATACAGGCTTGTGCGTTTTTTTTTTTtttaatt              >  1:6646740/1‑36 (MQ=37)
                      gtgATACAGGCTTGTGCGTTTTTTTTTTTttcaatt              >  1:2793193/1‑36 (MQ=37)
                      gtgATACAGGCTTGTGCGTTTTTTTTTCCttcaatt              <  1:5242120/36‑1 (MQ=255)
                      gtgATACAGGCTTGTGCGTTTTTTTTTCCttcaatt              <  1:290363/36‑1 (MQ=255)
                       tgATACAGGCTTGTGCGTTTTTTTTTCCttcaattt             <  1:7025970/36‑1 (MQ=255)
                         aTACAGGCTTGTGCGTTTTTTTTTTCttcaatttca           >  1:435459/1‑36 (MQ=255)
                         aTACAGGCTTGTGCGTTTTTTTTTCCttcaatttca           <  1:5271142/36‑1 (MQ=255)
                         aTACAGGCTTGTGCGTTTTTTTTTCCttcaatttca           <  1:5687276/36‑1 (MQ=255)
                         aTACAGGCTTGTGCGTTTTTTTTTCCttcaatttca           <  1:6250996/36‑1 (MQ=255)
                                cTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATg    >  1:266587/1‑36 (MQ=255)
                                 ttGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGt   >  1:4095229/1‑36 (MQ=255)
                                 ttGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGt   >  1:6008918/1‑36 (MQ=255)
                                  tgtgCGTTTTTTTTTCCTTCAATTTCAATAAATGTa  <  1:3216218/36‑1 (MQ=255)
                                   |                                  
AAACTGCAATAAATTCTTAAGTGTGATACAGGCTTGTGCGTTTTTTTTTCCTTCAATTTCAATAAATGTA  >  REL606/157959‑158028

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 27 ≤  ATCG  < 34 ≤  ATCG  < 37 ≤  ATCG