Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	850,698	C→A	N83K (AAC→AAA)	ompX →	outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,698	0	C	A	96.2%	259.4 / 2.8	80	N83K (AAC→AAA)	ompX	outer membrane protein X
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base A (44/33);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 5.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

CCGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACC  >  NC_000913/850602‑850794                                                                                                 |                                                                                                 ccGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAg                                                                                               >  2:1218294/1‑100 (MQ=255) ccGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAg                                                                                               >  2:641643/1‑100 (MQ=255)  cGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGt                                                                                              <  1:2273039/100‑1 (MQ=255)    cTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGtac                                                                                            <  1:1558314/100‑1 (MQ=255)    cTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCGTTCACTTACACCGTGAAAAGCGGTGCTGCAAGCACTGGTGACGAAAACAAATAACAGtac                                                                                            >  2:1740207/1‑100 (MQ=255)     tATGAAGAAGCCAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGta                                                                                             <  2:468090/98‑1 (MQ=255)     tATGAAGAAGACACCAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGtact                                                                                           >  1:1738508/1‑100 (MQ=255)     tATGAAGAAGACAACAGCCAGATGGGTGTGATCGGTTCTTTAACTTACACCGAGAAAAGCCGTACTGCTAGCTCTGGAAAAAAAAAAAAAAAACAGtact                                                                                           <  2:1490094/100‑1 (MQ=255)       tGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGtactac                                                                                         <  1:2109103/100‑1 (MQ=255)        gaagaaGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACg                                                                                        >  1:1338416/1‑100 (MQ=255)         aagaagACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACgg                                                                                       <  1:1752088/100‑1 (MQ=255)         aagaagACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACgg                                                                                       <  1:1321671/100‑1 (MQ=255)          agaagaCAACAGCCCGCTGGGTGTTATCGGCTCTTTCACCTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGc                                                                                      <  1:513968/100‑1 (MQ=255)          agaagaCAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGc                                                                                      >  1:1261802/1‑100 (MQ=255)           gaagaCAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCa                                                                                     <  2:1312019/100‑1 (MQ=255)            aagaCAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCAt                                                                                    >  1:1882925/1‑100 (MQ=255)                caacaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCAct                                                                                >  1:1674653/1‑100 (MQ=255)                 aacaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActg                                                                               <  2:87399/100‑1 (MQ=255)                  ncaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgc                                                                              >  1:1264407/2‑100 (MQ=255)                  acaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCCctgc                                                                              >  1:1360210/1‑100 (MQ=255)                  acaGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgc                                                                              >  1:1358797/1‑100 (MQ=255)                   caGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgct                                                                             >  1:1485674/1‑100 (MQ=255)                   caGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCActgct                                                                             <  1:641644/100‑1 (MQ=255)                     gCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGCCTACAACAAAAAACAGTACTATGGCATCACTGCTgg                                                                           >  2:297252/1‑100 (MQ=255)                     gCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTgg                                                                           >  1:1150514/1‑100 (MQ=255)                         gCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTccn                                                                       <  2:1264404/100‑2 (MQ=255)                            gggTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCg                                                                       >  2:225121/1‑97 (MQ=255)                               atctaTCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAc                                                                 <  2:390159/96‑1 (MQ=255)                                gtgATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAcc                                                                <  2:1882918/100‑1 (MQ=255)                                gtgATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAcc                                                                <  2:2187163/100‑1 (MQ=255)                                gtgATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAcc                                                                <  1:168466/100‑1 (MQ=255)                                  gATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGAGTACAACAAAAAACAGTACTACGGCAGCACTGCTTGTCCGGCGTaac                                                                >  2:714294/1‑96 (MQ=255)                                  gATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGc                                                              >  2:1555352/1‑100 (MQ=255)                                  gATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGc                                                              >  1:602161/1‑100 (MQ=255)                                   aTCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAACAAACAGTACTACGGCATCACTGCCAGTCCGGCTTACCGCa                                                             >  2:59365/1‑100 (MQ=255)                                   aTCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAcagat                                                             >  1:390161/1‑96 (MQ=255)                                   aTCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCa                                                             >  1:9116/1‑100 (MQ=255)                                   aTCGGTTCTTTCACTTACACCGAGAAAAGCCATACTGAAAGCACCGGTGCCTACCACAAAAAACAATACCACTGCGTCACTAATGGttaagctt                                                                   >  2:194348/1‑87 (MQ=255)                                     cGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCAtt                                                           <  2:1843978/100‑1 (MQ=255)                                      ggTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGCCCGGCTTACCGCATTa                                                          >  1:260241/1‑100 (MQ=255)                                         tCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACg                                                       <  2:1738503/100‑1 (MQ=255)                                          cTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGa                                                      >  1:390518/1‑100 (MQ=255)                                            ttCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGAc                                                     >  2:880798/1‑99 (MQ=255)                                             tCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTg                                                   <  2:1377651/100‑1 (MQ=255)                                              cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTgg                                                  >  2:2062490/1‑100 (MQ=255)                                                nttACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTCACGACTGGGc                                                >  1:856657/2‑100 (MQ=255)                                                cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGc                                                <  1:1555356/100‑1 (MQ=255)                                                cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGc                                                <  2:260239/100‑1 (MQ=255)                                                cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCACATTAACGACTGGGc                                                >  2:1926419/1‑100 (MQ=255)                                                 ttACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCa                                               <  1:1516739/100‑1 (MQ=255)                                                  tACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCaa                                              >  1:1247423/1‑100 (MQ=255)                                                   acacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAg                                             <  1:2335443/100‑1 (MQ=255)                                                    cacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGc                                            <  1:826734/100‑1 (MQ=255)                                                    cacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGc                                            <  2:1261799/100‑1 (MQ=255)                                                     acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCa                                           >  1:1915025/1‑100 (MQ=255)                                                     acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCa                                           <  1:2062497/100‑1 (MQ=255)                                                       cGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATc                                         >  1:1265261/1‑100 (MQ=255)                                                        gagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCt                                        >  2:522024/1‑100 (MQ=255)                                                         agaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTa                                       <  2:1974826/100‑1 (MQ=255)                                                               gCCGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGgtgt                                 <  2:1666292/100‑1 (MQ=255)                                                                 cGTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAg                               >  2:722791/1‑100 (MQ=255)                                                                  gTACTGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTAACGCATTAACGACTGGGCAAGCATCTACGGTGTAgt                              <  1:714295/100‑1 (MQ=255)                                                                      tGCAAGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggt                          >  2:1924190/1‑100 (MQ=255)                                                                         aaGCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtg                       <  2:1819274/100‑1 (MQ=255)                                                                           gCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                     <  2:2004208/100‑1 (MQ=255)                                                                           gCTCTGGTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtggg                     <  2:9116/100‑1 (MQ=255)                                                                              ctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTa                  >  2:1517524/1‑100 (MQ=255)                                                                                ggTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                >  1:660509/1‑100 (MQ=255)                                                                                ggTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                >  1:1401761/1‑100 (MQ=255)                                                                                ggTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                >  1:1157327/1‑100 (MQ=255)                                                                                ggTGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                >  1:1041030/1‑100 (MQ=255)                                                                                  tGACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGt              <  1:1740212/100‑1 (MQ=255)                                                                                   gACTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTa             <  2:1674648/100‑1 (MQ=255)                                                                                     cTACAACAAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa           >  1:1099793/1‑100 (MQ=255)                                                                                       acaacaAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt         >  1:1781318/1‑100 (MQ=255)                                                                                       acaacaAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt         >  1:1259320/1‑100 (MQ=255)                                                                                       acaacaAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt         >  1:1259113/1‑100 (MQ=255)                                                                                       acaacaAAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt         >  2:923821/1‑100 (MQ=255)                                                                                           caaAAAACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg     <  2:1401758/100‑1 (MQ=255)                                                                                             aaaaACAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc   >  2:2091573/1‑100 (MQ=255)                                                                                              aaaaCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAcc  <  1:1674537/100‑1 (MQ=255)                                                                                                 |                                                                                                 CCGCTATGAAGAAGACAACAGCCCGCTGGGTGTGATCGGTTCTTTCACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACC  >  NC_000913/850602‑850794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑