Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,439,798 A→C 16.7% V168V (GTT→GTG ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,439,7980AC16.7% 237.7 / 57.0 126V168V (GTT→GTGompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base A (63/42);  new base C (8/13);  total (71/55)
Fisher's exact test for biased strand distribution p-value = 9.09e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 131 total aligned reads displayed.

ACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAC  >  minE/1439729‑1439868
                                                                     |                                                                      
aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                       <  1:933613/71‑1 (MQ=255)
aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                       <  1:1318337/71‑1 (MQ=255)
aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                       >  1:3488270/1‑71 (MQ=255)
aCGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGaa                                                                       >  1:1347699/1‑71 (MQ=255)
   aCCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAg                                                                    >  1:3389849/1‑71 (MQ=255)
    ccGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGc                                                                   <  1:3671175/71‑1 (MQ=255)
       ttgttAGTTACGACACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCa                                                                  <  1:188447/69‑1 (MQ=255)
         gttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt                                                              >  1:3033371/1‑71 (MQ=255)
          ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:1015280/71‑1 (MQ=255)
          ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:2604015/71‑1 (MQ=255)
          ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:1725287/71‑1 (MQ=255)
          ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:503659/71‑1 (MQ=255)
          ttagttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:499130/71‑1 (MQ=255)
          ttagttaCGACACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGtt                                                             <  1:571956/71‑1 (MQ=255)
             gttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                           <  1:2935846/70‑1 (MQ=255)
             gttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                          <  1:798875/71‑1 (MQ=255)
             gttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                          <  1:601390/71‑1 (MQ=255)
              ttaCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTGCAgg                                                         <  1:1589190/71‑1 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaaa                                                                >  1:697794/1‑63 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                                        >  1:326855/1‑71 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                                        >  1:2099327/1‑71 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                                        >  1:2938734/1‑71 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGc                                                        >  1:3908289/1‑71 (MQ=255)
               taCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACACTGGTACTGAACAGCAAAGtt                                                             >  1:2765155/1‑66 (MQ=255)
                aCGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGcc                                                       >  1:3794144/1‑71 (MQ=255)
                    cACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTc                                                    >  1:2167676/1‑70 (MQ=255)
                      cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAAc                                                                      >  1:2335284/1‑50 (MQ=255)
                      cTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                 >  1:744548/1‑71 (MQ=255)
                       tAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAAc                                                 >  1:3632776/1‑70 (MQ=255)
                           aaaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGt                                                              <  1:1760896/53‑1 (MQ=255)
                           aaaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                          >  1:1087946/1‑57 (MQ=255)
                            aaGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCaaa                                                                >  1:1402010/1‑50 (MQ=255)
                              gCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCga                                         <  1:3173951/71‑1 (MQ=255)
                              gCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCg                                          >  1:2038521/1‑70 (MQ=255)
                              gCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCg                                          >  1:2685902/1‑70 (MQ=255)
                               ccTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCg                                          >  1:1339562/1‑69 (MQ=255)
                                   cacaAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:3147482/1‑71 (MQ=255)
                                   cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:2077659/1‑71 (MQ=255)
                                   cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:3946223/1‑71 (MQ=255)
                                   cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:1986726/1‑71 (MQ=255)
                                   cACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGaagaag                                    >  1:2249199/1‑71 (MQ=255)
                                     ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                           >  1:2549583/1‑46 (MQ=255)
                                     ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                          >  1:924170/1‑47 (MQ=255)
                                     ccAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAg                                                          <  1:1781687/47‑1 (MQ=255)
                                      cAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTc                                                            >  1:1971971/1‑44 (MQ=255)
                                      cAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                 <  1:3281883/71‑1 (MQ=255)
                                      cAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCa                                 <  1:1939185/71‑1 (MQ=255)
                                          tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                             <  1:134798/71‑1 (MQ=255)
                                          tGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                             <  1:1767127/71‑1 (MQ=255)
                                             gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTa                           >  1:3070716/1‑70 (MQ=255)
                                             gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                          >  1:2273931/1‑71 (MQ=255)
                                             gTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTAc                          >  1:3023141/1‑71 (MQ=255)
                                              ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCa                                                           >  1:3601290/1‑37 (MQ=255)
                                              ttGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCa                                                   <  1:532337/45‑1 (MQ=255)
                                                gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGa                                             >  1:3145691/1‑49 (MQ=255)
                                                gCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACgg                        >  1:2088991/1‑70 (MQ=255)
                                                 tcGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACga                        >  1:3865664/2‑68 (MQ=2)
                                                  cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACgat                       >  1:2134175/1‑67 (MQ=255)
                                                  cGTTGTTACCTTGATACTGCACCGCGAAGTTCAGACAATTTACCAGACCGAAGAAATCAGTATTACgat                       >  1:3357153/1‑67 (MQ=255)
                                                   gTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAgg                    >  1:2133154/1‑71 (MQ=255)
                                                    tttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGt                   >  1:1662705/1‑71 (MQ=255)
                                                     ttttACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAgtgt                             >  1:2405916/1‑60 (MQ=255)
                                                      gttACCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTc                  <  1:2924173/69‑1 (MQ=255)
                                                        tACCCTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgc                <  1:2593206/70‑1 (MQ=255)
                                                         aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg               <  1:1926335/70‑1 (MQ=255)
                                                         aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg               <  1:3891873/70‑1 (MQ=255)
                                                         aCCTTGATACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTcgcg               <  1:1496269/70‑1 (MQ=255)
                                                         aCCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                       <  1:1806140/62‑1 (MQ=255)
                                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg               >  1:1247060/1‑69 (MQ=255)
                                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg               >  1:202365/1‑69 (MQ=255)
                                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTcgcg               <  1:2765595/69‑1 (MQ=255)
                                                          cccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGa              <  1:1979738/70‑1 (MQ=255)
                                                          cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaac             <  1:2843971/71‑4 (MQ=255)
                                                          cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaac             <  1:1605507/71‑4 (MQ=255)
                                                          cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa              <  1:3432106/70‑3 (MQ=255)
                                                          cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgcaa              <  1:212156/70‑3 (MQ=255)
                                                          cccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGTtgca               <  1:3448962/69‑2 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTATGGTAGGTCGCGAAg            <  1:541593/71‑1 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                       >  1:320875/1‑60 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGt                       >  1:1132752/1‑60 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg            <  1:3569742/71‑1 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg            <  1:3639299/71‑1 (MQ=255)
                                                           ccTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg            <  1:2413742/71‑1 (MQ=255)
                                                           ccTGATATTGCAAAGCAAAGTTCAAGCCATCGACCTAACCAAAGAAGCCATTGTTACGGTAGGt                   <  1:2829234/64‑1 (MQ=255)
                                                             tgatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGc           >  1:2967278/4‑70 (MQ=9)
                                                             tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc          >  1:1207240/1‑71 (MQ=255)
                                                             tGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGcc          >  1:3453501/1‑71 (MQ=255)
                                                              gatACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca         >  1:2334103/3‑70 (MQ=255)
                                                               gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAg            >  1:2475107/1‑67 (MQ=255)
                                                               gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:277325/1‑71 (MQ=255)
                                                               gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:3323222/1‑71 (MQ=255)
                                                               gTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGt        >  1:2406348/1‑71 (MQ=255)
                                                               atACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca         <  1:956235/69‑2 (MQ=255)
                                                               atACTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGcca         <  1:2362503/69‑2 (MQ=255)
                                                                tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt       <  1:165070/71‑1 (MQ=255)
                                                                tACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGtt       <  1:294004/71‑1 (MQ=255)
                                                                 aCTGCACCGCGAAGTTCAGACCATTTACCAGACCGAAGAAATCAGTATTACGATAAGTCGCGACGCCAttg      >  1:1124359/1‑70 (MQ=255)
                                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCgaag                                       <  1:3823167/37‑1 (MQ=255)
                                                                  cTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTAc     <  1:1444877/71‑1 (MQ=255)
                                                                     aaCAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAc  >  1:3771385/1‑71 (MQ=255)
                                                                     |                                                                      
ACGACCGTTGTTAGTTACGCCACTAGTAAAGCCTTCACCAGATGGGTTGCCGTTTTTACCCTGGTACTGAACAGCAAAGTTCAGGCCGTCAACCAGACCGAAGAAGTCAGTGTTACGGTAGGTCGCGAAGCCGTTACCAC  >  minE/1439729‑1439868

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: