Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,228,129	C→T	7.2%	S310S (TCG→TCA)	PP_4606 ←	ferric siderophore receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,228,129	0	C	T	7.2%	158.7 / 2.0	55	S310S (TCG→TCA)	PP_4606	ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base T (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

AGGCGGCGCGTGGCAGGGATGTTCAGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCATAGGTCGAGCC  >  NC_002947/5227983‑5228275                                                                                                                                                   |                                                                                                                                                    aGGCGGCGCGTGGCAGGGATGTTCAGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGa                                                                                                                                                   >  8:301444/1‑149 (MQ=255)     ggcgCGTGGCAGGGATGTTCAGCGGGTGGCTGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑cg                                                                                                                                              >  2:86617/1‑149 (MQ=255)      gcgcGTGGCAGGGATGTTCAGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑cgc                                                                                                                                             >  8:85789/1‑149 (MQ=255)                    tGTTCAGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAg                                                                                                                               <  5:34988/149‑1 (MQ=255)                         aGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGACAGGCCAGCGACGGCGCCACCAGCGCTT‑CGCGGTTCACGCGGTAGTTCCg                                                                                                                          >  4:47484/1‑149 (MQ=255)                         aGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGACCATCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGAACGCCGTAGTTACg                                                                                                                          >  1:37490/1‑149 (MQ=255)                          gCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGc                                                                                                                         <  3:195021/149‑1 (MQ=255)                            gggTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGc                                                                                                                         >  5:110722/1‑147 (MQ=255)                                            cGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCg                                                                                                       <  7:301444/149‑1 (MQ=255)                                             gAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGt                                                                                                      >  3:6952/1‑149 (MQ=255)                                              caCGCGGTGCCACGGGCGAACGGATACCCAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCCTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTg                                                                                                     <  8:205239/148‑1 (MQ=255)                                              aaCGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTc                                                                                                           <  8:70002/143‑1 (MQ=255)                                              aaCGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTc                                                                                                           >  7:70002/1‑143 (MQ=255)                                                   ggTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGTGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTcga                                                                                                >  8:149224/1‑149 (MQ=255)                                                      gCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGAt                                                                                                                                                  <  4:65573/97‑1 (MQ=255)                                                      gCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGAt                                                                                                                                                  >  3:65573/1‑97 (MQ=255)                                                       ccACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTc                                                                                                  >  7:208316/1‑143 (MQ=255)                                                       ccACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTc                                                                                                  <  8:208316/143‑1 (MQ=255)                                                         aCGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATTCCGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGAcc                                                                                           <  1:86617/149‑1 (MQ=255)                                                             tCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCg                                                                                      >  7:169582/1‑149 (MQ=255)                                                               gAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCgg                                                                                     >  6:76012/1‑148 (MQ=255)                                                               gAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGt                                                                                    <  2:223820/149‑1 (MQ=255)                                                                aaCGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCgg                                                                                     >  8:36910/1‑147 (MQ=255)                                                                     aTAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCg                                                                               >  4:196644/1‑148 (MQ=255)                                                                         agaAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCagta                                                                                                                    <  6:59862/107‑1 (MQ=255)                                                                         agaAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCagta                                                                                                                    >  5:59862/1‑107 (MQ=255)                                                                         agaAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGtt                                                                          >  1:145030/1‑149 (MQ=255)                                                                           aaaTTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGc                                                                        <  4:250865/149‑1 (MQ=255)                                                                            aaTTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCt                                                                       <  3:196644/149‑1 (MQ=255)                                                                            aaTTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCt                                                                       <  3:138103/149‑1 (MQ=255)                                                                             aTTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTg                                                                      <  6:110722/149‑1 (MQ=255)                                                                             aTTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTg                                                                      >  2:212785/1‑149 (MQ=255)                                                                                 aCGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGTGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCAc                                                                  <  7:149224/149‑1 (MQ=255)                                                                                  cGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCAcc                                                                 <  4:6952/149‑1 (MQ=255)                                                                                  cGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCAcc                                                                 <  7:85789/149‑1 (MQ=255)                                                                                    gTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCAcc                                                                 >  3:15999/1‑147 (MQ=255)                                                                                          gTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGAAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTgc                                                           <  5:313015/147‑1 (MQ=255)                                                                                             ggCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGccc                                                       <  4:15999/148‑1 (MQ=255)                                                                                               ccAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTg                                                    >  6:256016/1‑149 (MQ=255)                                                                                                 aGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCt                                                  <  7:36910/149‑1 (MQ=255)                                                                                                   caccacCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTg                                                 >  4:204517/1‑148 (MQ=255)                                                                                                   caccacCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTg                                                 <  3:204517/148‑1 (MQ=255)                                                                                                      caccTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGTGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGccc                                                       <  8:207746/139‑1 (MQ=255)                                                                                                      caccTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGTGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGccc                                                       >  7:207746/1‑139 (MQ=255)                                                                                                          tGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGccc                                         >  4:294597/1‑149 (MQ=255)                                                                                                                    cGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCt                                <  5:256016/148‑1 (MQ=255)                                                                                                                     gCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGc                              <  2:145030/149‑1 (MQ=255)                                                                                                                       ccagccagGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCg                            >  8:157069/1‑149 (MQ=255)                                                                                                                            caggccagCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCtt                       <  3:294597/149‑1 (MQ=255)                                                                                                                               gccagCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTAGGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGcc                    >  1:231174/1‑149 (MQ=255)                                                                                                                               gccagCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGcc                    >  7:53246/1‑149 (MQ=255)                                                                                                                                       cGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCAt            <  2:231174/149‑1 (MQ=255)                                                                                                                                        ggCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTTGTCGACGACCAGGCTGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTTTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCATa           <  2:37490/149‑1 (MQ=255)                                                                                                                                               ccaGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCATAGGTCGAg    >  1:264729/1‑149 (MQ=255)                                                                                                                                                 aGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCATAGGTCGAGcc  >  5:128060/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGGCGGCGCGTGGCAGGGATGTTCAGCGGGTGGCCGTTGTTGCCGAACGCGGTGCCACGGTCGAACGGATAGAGAAATTCACGGTGCTCGTAGGCCAGCACCACCTGGGTGTCTTCGCCCAGCCAGGCCAGCGACGGCGCCACCAGCGATT‑CGCGGTGCACGCCGTAGTTACGCCAGTAGTCCTCGTCTTCGTGGTCGACGACCAGGCGGTAGGCGAAGTTGCTGTCACCCAGTGCGCCCGTGCTGTCGAGCCCGCCACCGCTGCCGTTCTTGCCGCTGCCATAGGTCGAGCC  >  NC_002947/5227983‑5228275

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