Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,813,385	A→G	N257S (AAC→AGC)	emrB →	multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,813,385	0	A	G	95.6%	129.1 / ‑3.1	45	N257S (AAC→AGC)	emrB	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (21/22);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CTGCTGGTTATTGGTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGG  >  NC_000913/2813237‑2813531                                                                                                                                                     |                                                                                                                                                   ctgctgGTTATTGGTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGcc                                                                                                                                                  >  2:530344/1‑151 (MQ=255)   gctgGTTATTGGTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGcc                                                                                                                                                  <  1:655932/149‑1 (MQ=255)       gTTATTGGTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAg                                                                                                                                            <  2:127865/151‑1 (MQ=255)          aTTGGGTTCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTGCCGTGGTGGCGGTGGTGGCGATCTGATTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCg                                                                                                                                         <  1:436343/151‑1 (MQ=255)             ggTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATc                                                                                                                                      <  2:164448/151‑1 (MQ=255)              gTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCt                                                                                                                                     >  2:456312/1‑151 (MQ=255)                aTCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGt                                                                                                                                   >  1:511235/1‑151 (MQ=255)                  cGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCg                                                                                                                                 >  2:375587/1‑151 (MQ=255)                                ttATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTcgcg                                                                                                                   >  2:531835/1‑151 (MQ=255)                                  aTGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCa                                                                                                                 <  2:320616/151‑1 (MQ=255)                                        gACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGATGAACGACGATAGCCCGATAGTCGATCTGGCGTTGTTTAAGTCGCGCAACTTc                                                                                                            >  2:301030/1‑150 (MQ=255)                                            gcgGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCAt                                                                                                       >  1:303835/1‑151 (MQ=255)                                             cgGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATc                                                                                                      <  2:649500/151‑1 (MQ=255)                                                 aaaGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGgct                                                                                                  <  2:177011/151‑1 (MQ=255)                                                  aaGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGgctg                                                                                                 <  1:92282/151‑1 (MQ=255)                                                      gCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCtt                                                                                              >  2:705718/1‑150 (MQ=255)                                                           aCTGGTTTTCATCACGGGAAATTATCAGCCTTACCGTGATGGCGGTGGTGGCTATCTGCTTCCTGATTTTCTGGGAGCCGACCGACGATAGCTCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGCCTGCTTGTGTa                                                                                         <  1:199817/150‑1 (MQ=255)                                                                 tttCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCt                                                                                  >  2:501733/1‑151 (MQ=255)                                                                     atcaCAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTcgcg                                                                              >  2:634157/1‑151 (MQ=255)                                                                     atcaCAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTcgcg                                                                              <  2:341653/151‑1 (MQ=255)                                                                       cacaGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGta                                                                            <  1:513703/151‑1 (MQ=255)                                                                           ggAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATg                                                                        <  2:299281/151‑1 (MQ=255)                                                                              aaTTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGctc                                                                     >  2:736596/1‑151 (MQ=255)                                                                               aTTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGctct                                                                    >  2:577782/1‑151 (MQ=255)                                                                                    catcCTTCCCGTGGTGGCGGTGGTGGCTATCTGTTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCTATCTATCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTc                                                               <  1:548656/151‑1 (MQ=255)                                                                                    catcCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTc                                                               >  2:30317/1‑151 (MQ=255)                                                                                       ccTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGgc                                                            >  2:464944/1‑151 (MQ=255)                                                                                          tACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCt                                                         >  2:286972/1‑151 (MQ=255)                                                                                              gtggtgGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTg                                                     <  2:736344/151‑1 (MQ=255)                                                                                                      ggtggtGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTctgctg                                             <  2:43617/151‑1 (MQ=255)                                                                                                       gtggtgGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCtgctgc                                            <  1:356479/151‑1 (MQ=255)                                                                                                             gCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCa                                        <  1:678648/149‑1 (MQ=255)                                                                                                                   tGCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGc                                >  1:424036/1‑151 (MQ=255)                                                                                                                    gCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGttgttgaa                               <  2:432804/151‑3 (MQ=255)                                                                                                                    gCTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGc                                <  2:511287/150‑1 (MQ=255)                                                                                                                     cTTCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGc                                >  1:35112/1‑149 (MQ=255)                                                                                                                      ttCCTGATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCagg                             >  2:446187/1‑151 (MQ=255)                                                                                                                           gATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAg                                       <  1:398566/136‑1 (MQ=255)                                                                                                                           gATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAg                                       >  2:398615/1‑136 (MQ=255)                                                                                                                           gATTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTGTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCGCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCaggagg                          <  1:43371/149‑1 (MQ=255)                                                                                                                            aTTGTCTGGGAGCTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCt                       <  2:721628/151‑1 (MQ=255)                                                                                                                                        cTGACCGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACgg           >  2:480235/1‑151 (MQ=255)                                                                                                                                            ccGACGATAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGAc       >  1:773783/1‑151 (MQ=255)                                                                                                                                                 gaTAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTggg  >  2:561074/1‑151 (MQ=255)                                                                                                                                                 gaTAGCCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTggg  <  2:293904/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CTGCTGGTTATTGGTATCGGCAGCCTGCAGATTATGCTCGACCGCGGTAAAGAGCTGGACTGGTTTTCATCACAGGAAATTATCATCCTTACCGTGGTGGCGGTGGTGGCTATCTGCTTCCTGATTGTCTGGGAGCTGACCGACGATAACCCGATAGTCGATCTGTCGTTGTTTAAGTCGCGCAACTTCACCATCGGCTGCTTGTGTATCAGCCTCGCGTATATGCTCTACTTCGGCGCTATTGTTCTGCTGCCGCAGTTGTTGCAGGAGGTCTACGGTTACACGGCGACCTGGG  >  NC_000913/2813237‑2813531

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