Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,522,583	C→T	Q161Q (CAG→CAA)	yncD ←	putative TonB‑dependent outer membrane receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,522,583	0	C	T	96.3%	83.4 / NA	27	Q161Q (CAG→CAA)	yncD	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/15);  minor base G (0/1);  total (11/16)
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

GCCCACCTGGGTCATCTGCTTTGATATCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGAC  >  NC_000913/1522312‑1522874                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                    gCCCACCTGGGTCATCTGCTTTGATATCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGAcattca                                                                                                                                                                                                                                                                         >  1:330734/1‑300 (MQ=255)    ccccTGGGTCATCTGCTTTGATATCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCNGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTAAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTa                                                                                                                                                                                                                                                                      >  1:513370/3‑300 (MQ=255)    cACCTGGGTCATCTGCTTTGATATCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATANCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTNTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGCAGTAACTACTGACTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTa                                                                                                                                                                                                                                                                      >  2:601233/1‑300 (MQ=255)                           tCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGTTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTTAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTTGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATAc                                                                                                                                                                                                                                               <  2:314590/300‑1 (MQ=255)                           tCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATAc                                                                                                                                                                                                                                               <  1:400618/300‑1 (MQ=255)                           tCCACACTATGGAGAATCAGACTTAATTTGCTGGCTTCACCAATGCGGACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCTTAAAGCGCGTAGTTGAGACGGTGTAATCGACATCGCCAGGCTGGGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATAc                                                                                                                                                                                                                                               <  2:223990/300‑1 (MQ=255)                           tCCACACGGTTGAAAAGCAGACTTAATTTGCTGGCTTCATCAATGCGTACGTCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCTTAAAACGCGTGGTTTAGACGGTGTAATCGACATCGCCAGGCTGGGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCTTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATAc                                                                                                                                                                                                                                               <  2:729994/300‑1 (MQ=255)                             cacaCTATTGAAAATCAGACTTAATTTTCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAg                                                                                                                                                                                                                                             <  2:446152/300‑1 (MQ=255)                             cacaCTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCTCCACTATGGTCACGATAGCCGTTGGTCGTAAAACGCGTGGTTTAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAg                                                                                                                                                                                                                                             <  1:639894/300‑1 (MQ=255)                                acTATTGAAAATCAGACTTGATTTGCTTGCTTCATCAATGCGGACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGc                                                                                                                                                                                                                                          <  2:569811/300‑1 (MQ=255)                                     tGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATNGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTTAGACGGTGTAATCTACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCNGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAgaga                                                                                                                                                                                                                                     <  2:686444/300‑1 (MQ=255)                                           tCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTTGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCNTAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCNGACGCGTTGCCATACAGGGCAGAGAATGGGc                                                                                                                                                                                                                               <  2:522887/300‑1 (MQ=255)                                                           ggCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCcac                                                                                                                                                                                                               >  1:377057/1‑300 (MQ=255)                                                                                            ggcgaaaTTTTTCTTTTCGCCACTATGGTCACGATAGCCGTGGGTCGGAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTNGCGCCCTTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGGTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGt                                                                                                                                                                              <  2:656932/296‑1 (MQ=255)                                                                                                                               atccGTGGTTCGTAAAACGCGTGGTTTATACGGTGTAATCGACATCACCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATAc                                                                                                                                           <  2:355227/298‑1 (MQ=255)                                                                                                                                gCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATAcc                                                                                                                                          >  1:32555/1‑300 (MQ=255)                                                                                                                                                  cgTGGTGGATACGGTGTACTCGACATCGCCAGGCTGTGTGCTGTCTCCCGTTGCNCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTGCTGGCTTCAATGGTTTGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGgcgaat                                                                                                                        <  2:545986/300‑1 (MQ=255)                                                                                                                                                                           tCGCCANGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCNAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCgg                                                                                               <  2:615377/300‑1 (MQ=255)                                                                                                                                                                                          ccGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCNATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGa                                                                                >  2:658610/1‑300 (MQ=255)                                                                                                                                                                                                               gCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATGCAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTgc                                                           >  2:392461/1‑300 (MQ=255)                                                                                                                                                                                                                   tCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCAt                                                       >  1:302416/1‑300 (MQ=255)                                                                                                                                                                                                                             gcgcCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACANTGCTTAAATCGATGTTGGATGTTTGCNCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCg                                             <  2:497641/300‑1 (MQ=255)                                                                                                                                                                                                                                                         ctactGGCTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACANGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCg                 >  1:505654/1‑300 (MQ=255)                                                                                                                                                                                                                                                                cTTCAATGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGGCATTCATTACCCAACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCTCGAATCGAAAGCTGTAAATCTTGCGCATAGTTCTGACGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTg          >  2:381404/1‑300 (MQ=255)                                                                                                                                                                                                                                                                    agtggttggtTGTTGTTTTCCGGTCTGGTTGTTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTc      <  2:290918/298‑1 (MQ=255)                                                                                                                                                                                                                                                                    aaTGGTTGGTGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAAGCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTc      <  2:291400/300‑1 (MQ=255)                                                                                                                                                                                                                                                                        gttggtGGTTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGAc  >  2:377625/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                                    GCCCACCTGGGTCATCTGCTTTGATATCCACACTATTGAAAATCAGACTTAATTTGCTGGCTTCATCAATGCGTACGCCCAGTTTGGCATTGGCTAAATTTTTCTGTGCGCCACTATGGTCACGATAGCCGTGGGTCGTAAAACGCGTGGTTGAGACGGTGTAATCGACATCGCCAGGCTGTGTGCCGTCTCCCGTTGCGCCCGTTGCTTTCAGCCCATAGCGCCAGCTGCCAAAACTGCCGTAGTAACTACTGGCTTCAATGGTTGGTGGCTGTTGTCCGGTCTGGGTGGTGACATTCATTACCCCACCAGACGCGTTGCCATACAGGGCAGAGAAGGGGCCACGCAGCACTTCCACATTTTGCACACTGCTTAAATCGATGTTGGATGTTTGCCCTTGCCCGTCGGGCATGGTGGCGGGAATACCGTCCACATACAGGCGAATACCGCGAATACCGTAAGTGGAGCGGGAGCCAAATCCGCGAATCGACAGCTGTAAATCTTGCGCATAGTTCTGCCGGTTTTGTACCTGCAAACCAGGCACGCCGGTCAGTGATTCGGAC  >  NC_000913/1522312‑1522874

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