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
*	minE	2,033,597	0	G	A	8.0%	314.4 / 11.4	99	S12N (AGC→AAC)	tolC	transport 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: ‑