Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,673,550	G→C	11.1%	G108R (GGG→CGG) ‡	PP_4137 →	outer membrane siderophore receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,673,550	0	G	C	11.1%	154.9 / 9.6	54	G108R (GGG→CGG) ‡	PP_4137	outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (26/22);  new base C (2/4);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCT  >  NC_002947/4673409‑4673695                                                                                                                                              |                                                                                                                                                   gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCa                                                                                                                                            >  8:199748/1‑149 (MQ=255) gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCa                                                                                                                                            >  6:376283/1‑149 (MQ=255)    gACATGCAGCGCTTGGGGCAAGGACGAGTGGCCAGCAACCAGGACGTGCTGGACTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGG‑CGCACgg                                                                                                                                         <  4:325340/149‑1 (MQ=255)    gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGG‑CGCACgg                                                                                                                                         <  1:11594/149‑1 (MQ=255)        tGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGG‑CGCACGGttga                                                                                                                                     <  2:54417/149‑3 (MQ=255)           aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGgt                                                                                                                                  <  7:199748/149‑1 (MQ=255)                 tGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTac                                                                                                                            <  1:298804/149‑1 (MQ=255)                    ggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTACCGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTacac                                                                                                                          >  5:297464/1‑148 (MQ=255)                                   ccAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGtt                                                                                                                                       <  1:387740/120‑1 (MQ=255)                                   ccAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGtt                                                                                                                                       >  2:387740/1‑120 (MQ=255)                                    cAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGc                                                                                                         >  1:149771/1‑149 (MQ=255)                                    cAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGG‑CGCACGGtt                                                                                                                                       <  8:206131/119‑1 (MQ=255)                                    cAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGG‑CGCACGGtt                                                                                                                                       >  7:206131/1‑119 (MQ=255)                                         aCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATCCGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTg                                                                                                    <  8:74531/149‑1 (MQ=255)                                         aCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGtt                                                                                                     >  4:300046/1‑148 (MQ=255)                                            aGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGt                                                                                                      >  1:206539/1‑144 (MQ=255)                                                cGTGCTGGCCTACCAGCCAGGGGTGGTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGTCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGc                                                                                              <  1:308500/149‑1 (MQ=255)                                                    cTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCggg                                                                                         >  2:73693/1‑149 (MQ=255)                                                         cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCa                                                                                    >  1:410904/1‑149 (MQ=255)                                                          tACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCAc                                                                                   >  5:262771/1‑149 (MQ=255)                                                           aCCAGCCAGGTGTGTTCGCCCAATCGGCATGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGt                                                                                                                                        <  6:62615/95‑1 (MQ=255)                                                           aCCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGAATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGt                                                                                                                                        >  5:62615/1‑95 (MQ=255)                                                               gccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCt                                                                               >  6:33964/1‑148 (MQ=255)                                                               cccagGGGGGTTCGCCCAAGCGGCAGGTAACGTGGCCATCCAGCTGCCTAAACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCtt                                                                              <  2:149771/148‑1 (MQ=255)                                                                ccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTa                                                                             <  1:73693/149‑1 (MQ=255)                                                                 cagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                            <  3:300046/149‑1 (MQ=255)                                                                       tgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGGCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTg                                                                      <  5:376283/149‑1 (MQ=255)                                                                            gCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCg                                                                  >  3:281023/1‑148 (MQ=255)                                                                              ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                               <  7:1628/149‑1 (MQ=255)                                                                                    ggCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtg                                                          <  7:81293/148‑1 (MQ=255)                                                                                         gTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                            <  3:39719/125‑1 (MQ=255)                                                                                         gTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                            >  4:39719/1‑125 (MQ=255)                                                                                              gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcgg                                                                                      >  5:416183/1‑110 (MQ=255)                                                                                              gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcgg                                                                                      <  6:416183/110‑1 (MQ=255)                                                                                                tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTg                                             >  5:64172/1‑149 (MQ=255)                                                                                                   cATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGcc                                           >  5:155097/1‑148 (MQ=255)                                                                                                    aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGa                                                                 <  8:177214/125‑1 (MQ=255)                                                                                                    aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGa                                                                 >  7:177214/1‑125 (MQ=255)                                                                                                            gTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACCCCATGTTCGCCGGCCTGCCGTTGCCCGGCCCggg                                                                                         >  5:22967/1‑93 (MQ=255)                                                                                                            gTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGCCCGGCCCggg                                                                                         >  5:29008/1‑93 (MQ=255)                                                                                                            gTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                 >  5:271432/1‑149 (MQ=255)                                                                                                               gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGt                               <  5:33964/148‑1 (MQ=255)                                                                                                                aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                         <  2:276402/121‑1 (MQ=255)                                                                                                                aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                         >  1:276402/1‑121 (MQ=255)                                                                                                                 tACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                            <  6:155097/149‑1 (MQ=255)                                                                                                                    gcgGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgcgc                          >  4:9979/1‑148 (MQ=255)                                                                                                                           gggCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGc                   <  6:297464/148‑1 (MQ=255)                                                                                                                             accggaaCCGTGCGCCCCCGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt                >  6:166692/6‑149 (MQ=255)                                                                                                                               cTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGa              >  8:351736/1‑149 (MQ=255)                                                                                                                                   aCCGTGCGCCGGGGG‑CGGACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATTGCTTTGACCaa          <  2:286557/149‑1 (MQ=255)                                                                                                                                      gTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                               >  3:324095/1‑93 (MQ=255)                                                                                                                                      gTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                               <  4:324095/93‑1 (MQ=255)                                                                                                                                           ccGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCt  <  8:329573/149‑1 (MQ=255)                                                                                                                                           ccGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCt  >  3:246652/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                   GTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGG‑CGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCT  >  NC_002947/4673409‑4673695

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