Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,845,787	T→G	10.0%	S26R (AGT→AGG)	pfeS‑II →	histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,845,787	0	T	G	10.0%	137.8 / 5.8	48	S26R (AGT→AGG)	pfeS‑II	histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (18/25);  new base G (2/3);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGC  >  NC_002947/1845640‑1845932                                                                                                                                                    |                                                                                                                                                  tCACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTACCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGtt                                                                                                                                                  <  5:54965/149‑1 (MQ=255)   aCGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTtat                                                                                                                                                <  6:318120/149‑1 (MQ=255)   aCGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATATTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTta                                                                                                                                                 >  7:273961/1‑148 (MQ=255)    cGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTtata                                                                                                                                               <  5:45789/149‑1 (MQ=255)      aaCGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATAcc                                                                                                                                             >  7:135328/1‑149 (MQ=255)            aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                       <  1:24241/149‑1 (MQ=255)            aGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                       >  2:24241/1‑149 (MQ=255)            aGATCCGCACCGTGGGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGGGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCCCAGGTATAACTGAGGc                                                                                                                                       <  7:85118/149‑1 (MQ=255)                     ccGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                              <  2:260187/149‑1 (MQ=255)                      cGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATg                                                                                                                              >  1:38203/1‑148 (MQ=255)                          tgGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAAc                                                                                                                         <  2:38203/149‑1 (MQ=255)                             ggCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCa                                                                                                                      <  5:142740/149‑1 (MQ=255)                               cAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAg                                                                                                                     >  8:40837/1‑148 (MQ=255)                                   ggTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGc                                                                                                                <  8:41321/149‑1 (MQ=255)                                    gTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGcc                                                                                                               <  6:151738/149‑1 (MQ=255)                                        cGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTc                                                                                                            <  7:40837/148‑1 (MQ=255)                                        cGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTc                                                                                                            <  1:43275/148‑1 (MQ=255)                                               aGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGa                                                                                                    <  5:330845/149‑1 (MQ=255)                                                gTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                               >  2:134752/1‑121 (MQ=255)                                                gTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACAt                                                                                                                               <  1:134752/121‑1 (MQ=255)                                                     aGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCgg                                                                                               >  3:30923/1‑148 (MQ=255)                                                              gAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGAc                                                                                     >  7:152799/1‑149 (MQ=255)                                                                   aaCCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGc                                                                                >  3:127128/1‑149 (MQ=255)                                                                      cATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGCTCGGCCTCAGGTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCggg                                                                             <  8:306839/149‑1 (MQ=255)                                                                                  ccATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATgccgcc                                                                  >  2:204573/1‑148 (MQ=255)                                                                                   cATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACGCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGa                                                                <  8:273961/149‑1 (MQ=255)                                                                                     ttCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAg                                                              >  5:336459/1‑149 (MQ=255)                                                                                            ttCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACa                                                       <  7:207384/149‑1 (MQ=255)                                                                                              cTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAgg                                                     >  5:313522/1‑149 (MQ=255)                                                                                               tGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                    >  2:208484/1‑149 (MQ=255)                                                                                                 gAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                    <  6:336459/147‑1 (MQ=255)                                                                                                 gAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGc                                                    <  3:96964/147‑1 (MQ=255)                                                                                                   aGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTgg                                                <  6:313522/149‑1 (MQ=255)                                                                                                      tGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGgcg                                              >  1:124363/1‑148 (MQ=255)                                                                                                        gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                           >  1:125806/1‑149 (MQ=255)                                                                                                        gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                           <  8:152799/149‑1 (MQ=255)                                                                                                        gCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAg                                           >  7:161049/1‑149 (MQ=255)                                                                                                         ccATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGt                                          <  1:208484/149‑1 (MQ=255)                                                                                                                 cTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCt                                                  >  3:285312/1‑133 (MQ=255)                                                                                                                 cTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCt                                                  <  4:285312/133‑1 (MQ=255)                                                                                                                               cctgcTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGAc                     <  4:30923/148‑1 (MQ=255)                                                                                                                                   ctgATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAATCTGAACAGGCCTGGCGCAGTGGCGGGCAGGCAGGGCTGGACCAgt                 <  2:124363/148‑1 (MQ=255)                                                                                                                                    tgatgaTCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                        <  8:147851/92‑1 (MQ=255)                                                                                                                                    tgatgaTCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTAt                                                                        >  7:147851/1‑92 (MQ=255)                                                                                                                                        tgccGGCCCCAGGTATAACTGAGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGTAGGGCTGGACCAgtgggtgg           >  7:111610/4‑149 (MQ=255)                                                                                                                                             gCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGcc      <  4:127128/149‑1 (MQ=255)                                                                                                                                               cTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCAt    <  7:45341/149‑1 (MQ=255)                                                                                                                                                tCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTCGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATg   >  6:245970/1‑149 (MQ=255)                                                                                                                                                  aGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTGGGAGCCGGCACGCCAGAGCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGTGCTGGACCAGTGGGTGGCGGCCTTGc  <  5:245970/148‑1 (MQ=255)                                                                                                                                                    ttATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggc                                       <  3:74441/109‑1 (MQ=255)                                                                                                                                                    ttATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTggc                                       >  4:74441/1‑109 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCACGAACGGCAGATCCGCACCGTGTGGGGCAAGGGTTACGTGCTCAGTGCCAGCGAGGCGGAGTGAACCATGCTCGACCGCCATTCGCTGTTCTGGAAGCTGGCCATCTTGCTGGTGGGCTTCTGCCTGCTGATGATCGGCCTCAGTTATACCTGGGGCCGGCACATGGAAACCCAGAATGCCTTCCTCTCGGAGCCGGCACGCCAGACCTTGCGGGGTTATGCCGCCGAAGCTGAACAGGCCTGGCGCAGTGGCGGGCGGGCAGGGCTGGACCAGTGGGTGGCGGCCATGC  >  NC_002947/1845640‑1845932

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