Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,172 A→G 43.6% N665N (AAT→AAC htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,1720AG43.6% 79.0 / 233.9 172N665N (AAT→AAChtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (56/41);  new base G (67/8);  total (123/49)
Fisher's exact test for biased strand distribution p-value = 3.41e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 173 total aligned reads displayed.

CAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAT  >  minE/138106‑138224
                                                                  |                                                    
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:2206783/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:1568992/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:2290667/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:2374680/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:1133973/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:2779105/1‑71 (MQ=255)
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                  >  1:28789/1‑71 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                              <  1:1003351/70‑1 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                              <  1:1551453/70‑1 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                              <  1:1543943/70‑1 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGa                                              <  1:2590104/70‑1 (MQ=255)
       ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGGTAt                                                     >  1:2030858/1‑61 (MQ=255)
       ccGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAAt                                           <  1:901224/71‑1 (MQ=255)
        cGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAAt                                           <  1:447126/70‑1 (MQ=255)
          cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                        <  1:2341473/71‑1 (MQ=255)
           tATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                        <  1:2260273/70‑1 (MQ=255)
                   aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCa                                >  1:775973/1‑70 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:856017/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1977668/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2051792/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2178894/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1055168/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:222403/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:936014/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:892599/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2355095/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2473130/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2527564/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:684222/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:667162/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:2837481/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:583445/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:353873/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1627541/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1376026/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1465011/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1471466/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1262145/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1545729/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:167672/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1189599/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:107293/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1084345/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1105746/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGaa                                             >  1:2896679/1‑56 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGaa                                       >  1:944650/1‑62 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGcc                                 >  1:207125/1‑68 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:959964/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2641511/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2700573/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2701597/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1330949/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2860918/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:449742/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:482651/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1870597/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:540301/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:845038/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1128169/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1174955/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:765346/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2079589/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2258532/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1779521/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1615602/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2059062/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2338475/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2009778/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:148548/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1810226/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:24330/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1406149/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:1867078/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:2532795/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTCGCAGAAATTGAACCTGCCAGc                              >  1:1082927/1‑71 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCCCTGTTTGCAGAAATTGAACCTGCCAg                               >  1:2607216/1‑70 (MQ=255)
                    cAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTAATCGCTGTTTGCAGAAATTGAACCTGCCAGc                              >  1:1228971/1‑71 (MQ=255)
                    cAAAACAACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGc                              >  1:3760/1‑71 (MQ=255)
                      aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGt                            >  1:2144523/1‑71 (MQ=255)
                          caccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTccc                        <  1:150523/71‑1 (MQ=255)
                          caccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTccc                        <  1:695308/71‑1 (MQ=255)
                                             aaCTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTccc                        <  1:2553373/52‑1 (MQ=255)
                                               cTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCa   >  1:342905/1‑71 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTTAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:1134972/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGGTGGCAt  <  1:5274/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:1706621/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2016212/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:971081/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:763186/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:753417/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:709993/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2124576/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:1659552/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:463708/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:1209741/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2919417/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:1232932/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2410488/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2660634/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt  <  1:2539557/69‑1 (MQ=255)
                                                                  |                                                    
CAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAT  >  minE/138106‑138224

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