Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,033,597 G→A 8.0% S12N (AGC→AAC)  tolC → transport channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,033,5970GA 8.0% 314.4 / 11.4 99S12N (AGC→AAC) tolCtransport channel
Reads supporting (aligned to +/- strand):  ref base G (41/50);  new base A (3/5);  total (44/55)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAG  >  minE/2033530‑2033666
                                                                   |                                                                     
gcgcTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCt                                                                    <  1:3750285/71‑1 (MQ=255)
   cTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCtt                                                                   <  1:810074/69‑1 (MQ=255)
   cTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCtt                                                                   <  1:3863011/69‑1 (MQ=255)
   cTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTc                                                                 >  1:3123980/1‑71 (MQ=255)
   cTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTc                                                                 >  1:174990/1‑71 (MQ=255)
     aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTg                                                               >  1:962205/1‑71 (MQ=255)
     aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTg                                                               >  1:1500036/1‑71 (MQ=255)
     aaaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCAGAGCCTTTCTg                                                               >  1:3996232/1‑71 (MQ=255)
      aaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTgg                                                              >  1:1647304/1‑71 (MQ=255)
      aaTACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTgg                                                              >  1:1561218/1‑71 (MQ=255)
        tACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGt                                                            <  1:2528857/71‑1 (MQ=255)
          cTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTc                                                          <  1:2513375/71‑1 (MQ=255)
          cTGCTTCACCACAAGGAATGCAAATGAAGAAATTGATCCCCATTCTTATCGGCCTGAGCCTTTCTGGGtt                                                           >  1:3392055/1‑70 (MQ=255)
           tGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTc                                                          <  1:977586/70‑1 (MQ=255)
            gCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAg                                                        >  1:2071620/1‑71 (MQ=255)
             cTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       >  1:187700/1‑71 (MQ=255)
             cTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       >  1:2282443/1‑71 (MQ=255)
             cTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       >  1:639713/1‑71 (MQ=255)
                caccacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTc                                                     <  1:2882429/70‑1 (MQ=255)
                caccacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTc                                                     <  1:2108319/70‑1 (MQ=255)
                 accacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGt                                                   <  1:3215184/71‑1 (MQ=255)
                 accacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGt                                                   >  1:149083/1‑71 (MQ=255)
                 accacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGt                                                   >  1:3488566/1‑71 (MQ=255)
                 accaaAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCg                                                    <  1:2959923/70‑1 (MQ=255)
                  ccacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGtt                                                  <  1:1739973/71‑1 (MQ=255)
                  ccacAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGtt                                                  >  1:1436213/1‑71 (MQ=255)
                     caAGGAATGCAAATGAAGAAGTTGCTCCCCATTCTTATCGGCCTGAGCCttt                                                                  >  1:2031101/1‑52 (MQ=255)
                     caAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCttt                                                                  >  1:2125529/1‑52 (MQ=255)
                     caAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCtt                                                                   <  1:199314/51‑1 (MQ=255)
                     caAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGa                                                <  1:2595890/70‑1 (MQ=255)
                      aatgAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTc                                                     <  1:1715720/61‑1 (MQ=255)
                      aaGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTc                                                     <  1:845086/64‑1 (MQ=255)
                      aaGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTc                                                     <  1:1707194/64‑1 (MQ=255)
                      aaGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGc                                              <  1:1949947/71‑1 (MQ=255)
                       aGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:160308/71‑1 (MQ=255)
                       aGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGc                                              >  1:3188914/1‑70 (MQ=255)
                         gAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCatt                                                       <  1:3879374/59‑3 (MQ=255)
                         gAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       <  1:1635507/59‑1 (MQ=255)
                         gAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       >  1:3505766/1‑59 (MQ=255)
                         gAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       <  1:1659990/59‑1 (MQ=255)
                         gAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGGTCTGTTCg                                                    <  1:4040813/62‑1 (MQ=255)
                          aaTGCAAATGAAGAAATTGCTCCCCGTTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:1591425/68‑1 (MQ=255)
                           aTGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGc                                         >  1:1911287/1‑71 (MQ=255)
                           aTGCAAATGAAGAAATTGCTACCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGc                                         >  1:3190579/1‑71 (MQ=255)
                            tGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAgg                                          <  1:1137945/69‑1 (MQ=255)
                            tGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGcc                                        <  1:909966/71‑1 (MQ=255)
                            tGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGcc                                        <  1:350044/71‑1 (MQ=255)
                               aaaTGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGcc                                        >  1:3719071/1‑68 (MQ=255)
                               aaaTGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGcc                                        >  1:2079744/1‑68 (MQ=255)
                               aaaTGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCga                                      <  1:1363873/70‑1 (MQ=255)
                                  tGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCt                                                                <  1:2652744/41‑1 (MQ=255)
                                   gaagaaATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAAcc                                 >  1:4091438/1‑71 (MQ=255)
                                      taaaTTGCTCGCCATTCTTATCGGCCTGAGCGTTTCTGGGTTCAGTTCGtt                                                  <  1:519122/50‑1 (MQ=255)
                                        aaTTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGAt                             >  1:2909593/1‑70 (MQ=255)
                                        aaTTGCTCCCCATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGaa                                   >  1:3837575/1‑64 (MQ=255)
                                        aaTTGCTCCCCATTCTTATCGGCCTGAACCATTCTg                                                               >  1:3282424/1‑36 (MQ=38)
                                        aaTTGCTCCCCATTCTTAACGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGa                              <  1:784224/69‑1 (MQ=255)
                                         aTTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:3229644/53‑1 (MQ=255)
                                         aTTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:382698/53‑1 (MQ=255)
                                         aTTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:1944349/53‑1 (MQ=255)
                                         aTTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGc                           <  1:2237467/71‑1 (MQ=255)
                                         aTTGCTCCCCACTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATg                            >  1:1636126/1‑70 (MQ=255)
                                          ttGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGt                                                       <  1:1031690/42‑1 (MQ=255)
                                          ttGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:450305/52‑1 (MQ=255)
                                           tGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGc                           >  1:3033060/1‑69 (MQ=255)
                                            gCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGc                                              >  1:1905550/1‑49 (MQ=255)
                                               ccccATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGtt                      >  1:1183371/1‑70 (MQ=255)
                                                cccATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTa                    <  1:4047703/71‑1 (MQ=255)
                                                cccATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGttt                     <  1:3699109/70‑1 (MQ=255)
                                                cccATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGttt                     <  1:3858526/70‑1 (MQ=255)
                                                 ccATTCTTATCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGttt                     <  1:71144/69‑1 (MQ=255)
                                                  cATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTg                                                 >  1:3339272/1‑40 (MQ=255)
                                                    ttCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCg                                                    <  1:1149294/35‑1 (MQ=255)
                                                    ttCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCg                                                    <  1:2456910/35‑1 (MQ=255)
                                                    ttCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTg                                                 >  1:3582819/1‑38 (MQ=255)
                                                    ttCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTg                                                 >  1:1958992/1‑38 (MQ=255)
                                                    ttCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCag                >  1:1031828/1‑71 (MQ=255)
                                                     tCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGtt                                                  <  1:1025482/36‑1 (MQ=255)
                                                      cTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCgaga                                    <  1:2462983/49‑1 (MQ=255)
                                                      cTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagc               <  1:1755818/70‑1 (MQ=255)
                                                      cTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagc               <  1:1337524/70‑1 (MQ=255)
                                                         aTCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATc                  >  1:3000817/1‑64 (MQ=255)
                                                         aTCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagcaag            <  1:1822384/70‑1 (MQ=255)
                                                         aTCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagcaag            <  1:2477198/70‑1 (MQ=255)
                                                         aTCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagcaag            <  1:1399649/70‑1 (MQ=255)
                                                         aTCGGCCTGAACCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAAACTGATGCAAGTTTATCagcaag            >  1:3195874/1‑70 (MQ=255)
                                                          tCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGcc                                             <  1:430236/36‑1 (MQ=255)
                                                             gCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGc       <  1:1860439/71‑1 (MQ=255)
                                                               cTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     <  1:4029215/71‑1 (MQ=255)
                                                               cTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     <  1:610316/71‑1 (MQ=255)
                                                               cTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     >  1:2815117/1‑71 (MQ=255)
                                                               cTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     >  1:1251934/1‑71 (MQ=255)
                                                                tGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTa                    <  1:3956174/55‑1 (MQ=255)
                                                                 gAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     >  1:3429881/1‑69 (MQ=255)
                                                                 gAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCt     >  1:1456043/1‑69 (MQ=255)
                                                                  aGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTa   <  1:1178260/70‑1 (MQ=255)
                                                                  aGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAg  >  1:2647128/1‑71 (MQ=255)
                                                                  aGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAg  >  1:707710/1‑71 (MQ=255)
                                                                  aGCCTTACTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCagcaagc           >  1:527803/1‑62 (MQ=255)
                                                                   |                                                                     
GCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAGCCTTTCTGGGTTCAGTTCGTTGAGCCAGGCCGAGAACCTGATGCAAGTTTATCAGCAAGCACGCCTTAG  >  minE/2033530‑2033666

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: ‑