Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,376,393	A→T	11.2%	L391Q (CTG→CAG)	dnaK ←	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,376,393	0	A	T	11.2%	133.0 / 6.2	52	L391Q (CTG→CAG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (27/19);  new base T (2/4);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTCTG  >  NC_002947/5376250‑5376538                                                                                                                                                 |                                                                                                                                                  gAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGGCAGCGGGCTGACGTCCAGCAGCCgg                                                                                                                                              >  2:41307/1‑149 (MQ=255)    tGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACg                                                                                                                                           >  4:94543/1‑149 (MQ=255)       tCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCt                                                                                                                                        >  8:196324/1‑149 (MQ=255)         aCGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCtt                                                                                                                                       >  3:11780/1‑148 (MQ=255)          cGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagcc                                                                                                                                     <  6:199748/149‑5 (MQ=255)          cGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCTGCTGCTGGACGTcagcc                                                                                                                                     <  6:196107/149‑5 (MQ=255)                    gTTGTCGTCGGAAGTCGAGAAAACCTGCGACTTCTTGGTCGGTATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGGTTCGATACCCAGGTTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCg                                                                                                                            <  7:251830/148‑1 (MQ=255)                        tcgtcgTCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAg                                                                                                                       <  7:235307/149‑1 (MQ=255)                          gtcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTa                                                                                                                     <  6:272604/149‑1 (MQ=255)                           tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAc                                                                                                                    <  7:196324/149‑1 (MQ=255)                           tcgGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTAc                                                                                                                    >  2:51962/1‑149 (MQ=255)                                    gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCAcc                                                                                                              >  7:150835/1‑146 (MQ=255)                                    gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCAcc                                                                                                              <  8:150835/146‑1 (MQ=255)                                    gagaACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTg                                                                                                           >  5:46777/1‑149 (MQ=255)                                           ccTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGgca                                                                                                    <  8:13381/149‑1 (MQ=255)                                             tGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                                                               <  7:3942/104‑1 (MQ=255)                                             tGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCcagcagcag                                                                                                                                               >  8:3942/1‑104 (MQ=255)                                                   ttcttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCAt                                                                                            >  3:245181/1‑149 (MQ=255)                                                     cttGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATg                                                                                           <  3:198418/148‑1 (MQ=255)                                                          tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCt                                                                                                            >  2:49539/1‑126 (MQ=255)                                                          tCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCt                                                                                                            <  1:49539/126‑1 (MQ=255)                                                               taggtggtGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTc                                                                                <  7:105360/147‑1 (MQ=255)                                                               aTGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTc                                                                                <  4:11780/149‑1 (MQ=255)                                                                   tggtgTTCTTCTCGATCAGAGCGGTCATCACGCCGCCCAT‑GGTTTCGATACACAGGGTCAGCGGGCTGACGTCCATCAGCAGGACGTCTTTCACGTCACCGGCAAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTcc                                                                            >  1:10140/1‑149 (MQ=255)                                                                      tgtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCATAGGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTc                                                                                >  1:48810/1‑143 (MQ=255)                                                                      tgtTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTc                                                                                <  2:48810/142‑1 (MQ=255)                                                                              tcGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTc                                                                 >  2:246353/1‑149 (MQ=255)                                                                                gATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCt                                                               >  8:224444/1‑149 (MQ=255)                                                                                gATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCt                                                               >  1:61048/1‑149 (MQ=255)                                                                                             ctcccgccgcCCAT‑GGTTTCGAAACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCttt                                                  <  1:41307/145‑1 (MQ=255)                                                                                             atcaCGCCGCCCAT‑GGTTTCGATACACAGGGTCAGCGGGCTGACGTCCAGCAGCAGGAAGTCTTTCACGTCACCGGCCAGAACCGCAACCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTAGCGCGcttctt                                                   >  8:96218/1‑148 (MQ=255)                                                                                               caCGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTAc                                                >  1:43045/1‑149 (MQ=255)                                                                                                 cgccgcCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCg                                              <  6:46777/149‑1 (MQ=255)                                                                                                  gccgccCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGACAGTACCGCACACTGGAAGGCAGCAGCCATGGCAACGGCTTCGTACGGGTTGACGACCTTGCTCGCTTCTTTACCg                                              >  6:164665/1‑148 (MQ=255)                                                                                                      cccAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCgaagaa                                         <  4:245181/149‑1 (MQ=255)                                                                                                             gTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGc                                   >  5:253965/1‑149 (MQ=255)                                                                                                              tttCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTTGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGc                                   >  2:66994/1‑148 (MQ=255)                                                                                                               ttCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCa                                 >  1:240251/1‑149 (MQ=255)                                                                                                                 cGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCAc                                <  7:224444/148‑1 (MQ=255)                                                                                                                       ccAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACttctttc                         >  8:25433/1‑149 (MQ=255)                                                                                                                         aGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCTACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTg                       >  6:223462/1‑149 (MQ=255)                                                                                                                         aGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTTCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTg                       >  3:43508/1‑149 (MQ=255)                                                                                                                         aGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTg                       <  1:246353/149‑1 (MQ=255)                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCTACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGc                <  5:223462/149‑1 (MQ=255)                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgcgcgc                                                        >  4:136908/1‑109 (MQ=255)                                                                                                                                gCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTgcgcgc                                                        <  3:136908/109‑1 (MQ=255)                                                                                                                                  gggCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGTGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCgg              >  8:241396/1‑149 (MQ=255)                                                                                                                                   gaaagACGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCgg              >  6:105397/5‑148 (MQ=255)                                                                                                                                        aCGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGcttctt                                                   <  7:8692/106‑1 (MQ=255)                                                                                                                                        aCGTCCTGCTGCTGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGcttctt                                                   >  8:8692/1‑106 (MQ=255)                                                                                                                                          gTCCAGCTGCAGGTCGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTTCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCgg       <  7:40492/148‑1 (MQ=255)                                                                                                                                            ccAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCa                                 <  8:179257/120‑1 (MQ=255)                                                                                                                                            ccAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCa                                 >  7:179257/1‑120 (MQ=255)                                                                                                                                            ccAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTc    >  6:55008/1‑149 (MQ=255)                                                                                                                                             cagcagcagGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCACTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTc    <  3:94543/148‑1 (MQ=255)                                                                                                                                              agcagcagGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTCt   >  6:18509/1‑148 (MQ=255)                                                                                                                                              agcagcagGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTCTg  >  6:22323/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GAATGGTCACGGCGCTCTGGTTGTCGTCGGCAGTCGAGAACACCTGCGACTTCTTGGTCGGGATGGTGGTGTTCTTCTCGATCAGCGCGGTCATCACGCCGCCCAT‑GGTTTCGATACCCAGGGTCAGCGGGCTGACGTCCAGCAGCAGGACGTCTTTCACGTCACCGGCCAGTACCGCACCCTGGATGGCAGCACCCATGGCAACGGCTTCGTCCGGGTTGACGTCCTTGCGCGCTTCTTTACCGAAGAAGTCGGCCACTTCTTTCTGCACCAGCGGCATGCGGGTCTG  >  NC_002947/5376250‑5376538

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: ‑