Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	850,742	A→G	D98G (GAC→GGC)	ompX →	outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,742	0	A	G	98.4%	218.1 / ‑4.5	63	D98G (GAC→GGC)	ompX	outer membrane protein X
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/32);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

TACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTT  >  NC_000913/850651‑850837                                                                                            |                                                                                                cgatctgagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCaa                                                                                         <  2:875587/94‑1 (MQ=255)  acacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAg                                                                                        <  1:889498/100‑1 (MQ=255)  acacCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAg                                                                                        <  1:1296225/100‑1 (MQ=255)  acacCGAGAAAAGCCGTACTACAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAg                                                                                        <  2:762123/100‑1 (MQ=255)   cacCGAGAAAAGCCGTACTGCAAGCTATGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTACGGATTACCGCATTAACGCCTGGGaaa                                                                                         >  2:163411/1‑98 (MQ=255)    acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCa                                                                                      <  2:12051/100‑1 (MQ=255)    acCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCa                                                                                      <  2:460981/100‑1 (MQ=255)       gagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCaaagatcg                                                                                   >  1:875594/1‑94 (MQ=255)          aaaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACg                                                                                <  2:264099/100‑1 (MQ=255)             aGCCGTACTGCAAGCTCTGGTGACTACACCAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgtg                                                                             >  2:1286489/1‑100 (MQ=255)             aGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgtg                                                                             <  2:798811/100‑1 (MQ=255)                  tACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtg                                                                        <  1:1451634/100‑1 (MQ=255)                    cTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtggg                                                                      <  1:1353893/100‑1 (MQ=255)                      gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtg                                                                    <  1:1118124/100‑1 (MQ=255)                      gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtg                                                                    >  1:998482/1‑100 (MQ=255)                       cAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgt                                                                   >  2:200915/1‑100 (MQ=255)                          gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCCAGCATCTACGGTGTAgtgggtgtggg                                                                >  2:411693/1‑100 (MQ=255)                          gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgtgg                                                                 >  2:341104/1‑99 (MQ=255)                           ctctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgtgggt                                                               >  2:995376/1‑100 (MQ=255)                              tGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTAt                                                            <  2:1039923/100‑1 (MQ=255)                               ggTGACTACAACAAAAACCAGTACTACGGCATGACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGGGGGTTATg                                                           >  1:1492026/1‑100 (MQ=255)                               ggTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGAATCTACGGTGTAGTGGGTGTGGGTTATg                                                           >  1:678770/1‑100 (MQ=255)                                    cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                      >  1:1099353/1‑100 (MQ=255)                                    cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                      >  2:257818/1‑100 (MQ=255)                                     tACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAt                                                     <  1:1432972/100‑1 (MQ=255)                                      acaacaAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCTAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt                                                    <  2:966161/100‑1 (MQ=255)                                       caacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTc                                                   <  1:38224/100‑1 (MQ=255)                                       caacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTc                                                   <  2:1348289/100‑1 (MQ=255)                                       caacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGGTATGGTAAAtt                                                    >  2:27448/1‑99 (MQ=255)                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                 <  2:146586/100‑1 (MQ=255)                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                 >  1:1323375/1‑100 (MQ=255)                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                 >  2:1283951/1‑100 (MQ=255)                                         acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                 >  2:1283965/1‑100 (MQ=255)                                          caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                >  2:40190/1‑100 (MQ=255)                                            aaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAc                                              <  1:1390641/100‑1 (MQ=255)                                               aCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCAc                                           >  2:13006/1‑100 (MQ=255)                                                 ccgatctacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  1:1256535/95‑1 (MQ=255)                                                 cAGTACTACGGCATCGCTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  2:214205/100‑1 (MQ=255)                                                 cAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  2:91236/100‑1 (MQ=255)                                                  aGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGa                                        <  1:1306255/100‑1 (MQ=255)                                                  aGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGa                                        >  2:498897/1‑100 (MQ=255)                                                   gTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGaa                                       >  2:155677/1‑100 (MQ=255)                                                     actacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAAt                                      <  1:8380/99‑1 (MQ=255)                                                     actacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATa                                     >  1:345032/1‑100 (MQ=255)                                                      ctacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAc                                    >  2:530108/1‑100 (MQ=255)                                                       tacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAcc                                   <  2:76144/100‑1 (MQ=255)                                                        acGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGagatcgg                                  >  2:1256527/1‑94 (MQ=255)                                                        acGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATAcca                                  >  1:76144/1‑99 (MQ=255)                                                         cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 >  2:1362028/1‑100 (MQ=255)                                                         cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 <  1:1334259/100‑1 (MQ=255)                                                            cATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAcc                              <  1:40190/100‑1 (MQ=255)                                                              tCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTa                            <  1:1495725/100‑1 (MQ=255)                                                               cACTGCTGGTCCGGCTTACAGCATTAACGGCTGGGCAAGCATTTACGGTGTAGTGGGGGTGGGGTATGGTAAATTTCAGACCACTAAATAccc                                  >  2:1288515/1‑93 (MQ=255)                                                                 ctgctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACaa                         >  1:636545/1‑100 (MQ=255)                                                                   gctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAac                       <  1:970458/100‑1 (MQ=255)                                                                   gctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAac                       >  1:140308/1‑100 (MQ=255)                                                                   gctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAac                       >  1:139872/1‑100 (MQ=255)                                                                      ggTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    <  1:163411/100‑1 (MQ=255)                                                                       gTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    <  1:554489/99‑1 (MQ=255)                                                                           ggCTTACCGCATTAACGGCTGGGCAAGCATCTACGGAGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACAcc               >  1:1111327/1‑100 (MQ=255)                                                                            gCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCa              <  2:948164/100‑1 (MQ=255)                                                                                       tAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGt   >  2:176453/1‑100 (MQ=255)                                                                                        aaCGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGtt  <  2:345031/100‑1 (MQ=255)                                                                                            |                                                                                                TACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTT  >  NC_000913/850651‑850837

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