Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	25,685	G→A	R160K (AGA→AAA)	lspA →	lipoprotein signal peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	25,685	0	G	A	97.1%	98.4 / ‑5.1	34	R160K (AGA→AAA)	lspA	lipoprotein signal peptidase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/17);  total (16/18)
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.72e-01

TGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTAC  >  NC_000913/25448‑25908                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                tGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTCCCTTCTAAAGCGAAAAAACa                                                                                                                                                                                                                     >  2:406541/1‑250 (MQ=255)    ccttggctgggaggatgTATCGCTCGAAGGCCCCGCAGAAGCAAAACAATATCGCTTACGCGCTGTTTAGTGGGGGCGCGCTGGGCACCCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACaata                                                                                                                                                                                                                  <  2:242035/237‑1 (MQ=255)                            cGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGATGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGCAGGGTTTTTACCTTCTAAAGGGAAAAAACAATAATAAACCCTGCCGGATTCGATGCTg                                                                                                                                                                                         >  2:319913/1‑251 (MQ=255)                             gagggCCACGTATAACCTAATCACTATGGCTTACGCGCTTATTATTGGCGGCTCCCTGCTCACCCTTTTCGGCCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCTCTTCGCCACCTTCAACCGTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGa                                                                                                                                                                                        <  2:39271/248‑1 (MQ=255)                                       cAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGACAAAACAATAATAAACCCTGCCGGATGCGATGCTTACGCATCTTAt                                                                                                                                                                              >  1:14212/1‑251 (MQ=255)                                        agaagaTAAACAATATCGCTTACGCGCTGTTTATGGGGGGCGCGCGGGGCACCCTGTCCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATc                                                                                                                                                                             <  1:172820/251‑1 (MQ=255)                                         gAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCCCTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAACCAATAATAAACCCTGCCGGATGCGATGCTTACGCATCTTATcc                                                                                                                                                                            >  2:30149/1‑251 (MQ=255)                                          aaGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGCTGCTGACGCATCTTATCCg                                                                                                                                                                           >  1:314244/1‑251 (MQ=255)                                                  cAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGTTACTGCCATCTGTGTCGGTGCGGCCCTGATTGTGCTGGAAGGTTTTTTACCTTCTTAAGCTAAACAACACTAATAAAACCTGCCGGCTGCGATGCTGAcccctcttttcctgccacc                                                                                                                                                                     >  2:113227/1‑231 (MQ=255)                                                                          gggcggggcgggggtaatccTGTTCGGCCGCCTGTTGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGa                                                                                                                                           <  1:295754/233‑1 (MQ=255)                                                                                cgcgcTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGGGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATTCTGACGCATCTTATCCGGCCTCAAGGTTGCTTCGAAACAGTAGGCCGGATAAGGc                                                                                                                                     >  2:84051/1‑251 (MQ=255)                                                                                                 cggcccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGGCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                    <  2:376832/247‑1 (MQ=255)                                                                                                     ccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAACAAATAATAAACCCTGACGGATGCGATGCTGACGCATCTTATACGGCCTACAGATTACTGCGAAATCATAGGCCGGCTAAAGCGTTTAAGCCTCCTCCGGcca                                                                                                                 >  2:366818/1‑248 (MQ=255)                                                                                                     ccGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaa                                                                                                                <  1:161514/251‑1 (MQ=255)                                                                                                     ccGCCTGTGGCACGGCTTCGTTGGCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaaa                                                                                                                <  2:89875/251‑1 (MQ=255)                                                                                                       gCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGCTGCGATGCTGACGCCTCTTATCCGGCCAACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACCCCGACTCCGGCacaa                                                                                                               >  2:410767/1‑247 (MQ=255)                                                                                                           gtgGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCCCTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGGAAAAACAATAATAAACCCTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGTCGTTTACGCCGCATCCGGCCAAAATCCt                                                                                                          >  2:132337/1‑251 (MQ=255)                                                                                                           gtgGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCTTTGACGCCGCGGCCGGCAAAAATcc                                                                                                           <  1:373976/250‑1 (MQ=255)                                                                                                                  ggCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAatat                                                                                                   <  1:423488/251‑1 (MQ=255)                                                                                                                            gcGATATGATCGACTTCTAAGTCGTCTGCTGGCATTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAAc                                                                                         <  1:84051/250‑1 (MQ=255)                                                                                                                             cGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                                                                                                >  1:105939/1‑243 (MQ=255)                                                                                                                             cGATATGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                                                                                                <  2:105939/243‑1 (MQ=255)                                                                                                                                  tGATCGACTTCTACGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCAt                                                                                   >  1:391053/1‑251 (MQ=255)                                                                                                                                              aCGTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGGGCGGCACTGATTGTGCTGGAAGGTTTTTTCCCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCt                                                                         <  1:232411/249‑1 (MQ=255)                                                                                                                                                gTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCCAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGCATATGTAc                                                                     >  2:139734/1‑251 (MQ=255)                                                                                                                                                gTCGGCGGCTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAc                                                                     <  2:453811/251‑1 (MQ=255)                                                                                                                                                                           aaCCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACt                                          >  2:81496/1‑251 (MQ=255)                                                                                                                                                                              cgtGCCGATACTGCCATATGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGAGGCGATGCTGACGCATCTTATCCGGCCTACAGATTGATGCGAAATCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAGCCTGCATGTCTGAATCTGTACAGTGCAATAGCGCCGTCCTGGTGCACTTCa                                       <  2:242074/249‑1 (MQ=255)                                                                                                                                                                                          gCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcg                           >  1:139488/1‑251 (MQ=255)                                                                                                                                                                                           ccATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcga                          <  1:74196/251‑1 (MQ=255)                                                                                                                                                                                                      ggTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGc                                                                                                                                               <  2:424966/123‑1 (MQ=255)                                                                                                                                                                                                      ggTGCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGc                                                                                                                                               >  1:424966/1‑123 (MQ=255)                                                                                                                                                                                                         gCGGCACTGATTGTGCTGGAAGGTTTTTTACCTTCTAAAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCaccacc             <  1:81496/250‑1 (MQ=255)                                                                                                                                                                                                                    ggggCTGGAAGGTTTTTTACCTTCTAAAGCAAAAAAACAATAAAAACCCCGGCCGAAGGCGATGCGTACGCATCTTATCCGGCCTACAGATTGCTGCAAAATCGTAGGCCGGATAAGGCGTTTCCGCCGCACCCGGAAAAAATCCTAAAAAATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCCCCGTCCTGGTGCACTTCACGCTAAAACCCGACGATGGCACCACCGCCGAGTCTAc  <  2:165093/247‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                TGATCCTGGCAGTGATGATGTATCGCTCGAAGGCCACGCAGAAGCTAAACAATATCGCTTACGCGCTGATTATTGGCGGCGCGCTGGGCAACCTGTTCGACCGCCTGTGGCACGGCTTCGTTGTCGATATGATCGACTTCTACGTCGGCGACTGGCACTTCGCCACCTTCAACCTTGCCGATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTAC  >  NC_000913/25448‑25908

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