Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 139,705 T→C 16.7% Q154Q (CAA→CAG htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE139,7050TC16.7% 224.0 / 28.0 136Q154Q (CAA→CAGhtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base T (69/44);  new base C (10/13);  total (79/57)
Fisher's exact test for biased strand distribution p-value = 1.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 136 total aligned reads displayed.

TTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTT  >  minE/139643‑139772
                                                              |                                                                   
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:2946808/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:1706373/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:1905314/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               <  1:2070033/69‑1 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:2232853/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:268126/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:2823959/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:1281586/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:3298814/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               <  1:100206/69‑1 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:36701/1‑69 (MQ=255)
ttAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGa                                                               >  1:1180834/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTg                                                                   >  1:2939106/1‑63 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGaa                                                              >  1:3014777/1‑68 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:1071774/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:90514/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:1627833/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:654654/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:2285136/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:2376289/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:2443457/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:2611467/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:62521/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:1364024/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:126543/1‑69 (MQ=255)
  aaTGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAAc                                                             >  1:3224452/1‑69 (MQ=255)
     gCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTCGAGGAACGTc                                                          <  1:1053444/69‑1 (MQ=255)
      ccGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTc                                                          >  1:933159/1‑68 (MQ=255)
      ccGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTc                                                          >  1:1623578/1‑68 (MQ=255)
                 tAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              <  1:1648354/69‑1 (MQ=255)
                 tAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              <  1:2359837/69‑1 (MQ=255)
                 tAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              <  1:1644478/69‑1 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCa                                                  >  1:1352630/1‑63 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:2285845/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:2900362/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:1346478/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:3416418/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:1340056/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:1603957/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:2221525/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:3056430/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:2052282/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:1947245/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGtt                                            >  1:3239482/1‑69 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                                             >  1:2103966/1‑68 (MQ=255)
                   aCGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                                             >  1:1158495/1‑68 (MQ=255)
                   aCGATGGATAAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGt                                             >  1:2943867/1‑68 (MQ=255)
                     gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:2735954/68‑1 (MQ=255)
                     gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:1932299/68‑1 (MQ=255)
                     gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:1818288/68‑1 (MQ=255)
                     gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:2649183/68‑1 (MQ=255)
                     gATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:1512616/68‑1 (MQ=255)
                      aTGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:1424009/67‑1 (MQ=255)
                      aTGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:1819848/67‑1 (MQ=255)
                        ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:3301159/65‑1 (MQ=255)
                        ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:2216036/65‑1 (MQ=255)
                        ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:661862/65‑1 (MQ=255)
                        ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTg                                           <  1:2245822/65‑1 (MQ=255)
                        ggATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATc                                        >  1:1038818/1‑68 (MQ=255)
                           tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              >  1:1737039/1‑59 (MQ=255)
                           tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              >  1:1881248/1‑59 (MQ=255)
                           tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACg                                              >  1:778831/1‑59 (MQ=255)
                           tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGtt                                     <  1:78475/68‑1 (MQ=255)
                           tCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGtt                                     <  1:811234/68‑1 (MQ=255)
                                aTAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACg                                 >  1:756191/1‑67 (MQ=255)
                                  aGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTATCGTCAt                             <  1:442784/69‑1 (MQ=255)
                                  aGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCAt                             <  1:2103968/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGTGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:749256/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:1542301/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:1157484/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:2716367/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:1444888/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:3377166/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:205781/69‑1 (MQ=255)
                                   gTTTTGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:1022105/69‑1 (MQ=255)
                                   gTTTTGATAATTCATCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:2322446/69‑1 (MQ=255)
                                    ttttGATAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTa                            <  1:2030272/68‑1 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:558154/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:3187997/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:3082853/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:3069268/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:2896663/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:2474370/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:2365083/2‑68 (MQ=255)
                                         atAATTCTTCATAATCCATGCCTGGGGGACATCGATATCCAGGCGTTGATCGTTAACATCGTAACGaa                       >  1:866857/2‑68 (MQ=255)
                                              ttttCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGaa                  <  1:1530033/68‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:559956/69‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:1014071/69‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:1226880/69‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:1256283/69‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:1861666/69‑1 (MQ=255)
                                                      aCCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAggg         <  1:1713788/69‑1 (MQ=255)
                                                           ggCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATa      >  1:2675347/1‑67 (MQ=255)
                                                           ggCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAt    >  1:1825668/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:2977560/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:274251/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:2601714/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:2467767/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:1611192/1‑69 (MQ=255)
                                                             cTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAAttt  >  1:2161066/1‑69 (MQ=255)
                                                              |                                                                   
TTAATGCCGTTTTCCCATAACGATGGATCAACATAGTTTTGGTAATTTTTCATTACCCAGGCTTGAGGAACGTCTATATCCAGACGTTGATCGTTAACGTCATAACGAACAGAAGCCTGAGGGATAATTT  >  minE/139643‑139772

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: