Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	322,140	T→C	24.1%	intergenic (+427/+31)	aguA → / ← PP_0267	agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	322,140	0	T	C	24.1%	97.1 / 30.3	58	intergenic (+427/+31)	aguA/PP_0267	agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (7/7);  total (29/29)
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.63e-01

AGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGC  >  NC_002947/321996‑322284                                                                                                                                                   |                                                                                                                                                 aGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                              <  2:179748/149‑1 (MQ=255)   cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                            <  2:235830/149‑1 (MQ=255)    aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                              >  2:209550/1‑146 (MQ=255)     tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCAGTGTGGCACTGATCGATTGCAGCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                          <  7:9361/149‑1 (MQ=255)        gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGGCTCTTCGCGGGCATGCTCGCGAAGAGGc                                                                                                                                       <  1:209550/149‑1 (MQ=255)        gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGc                                                                                                                                       <  3:352512/149‑1 (MQ=255)        gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGa                                                                                                                                             >  6:177589/1‑143 (MQ=255)        gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGa                                                                                                                                             <  5:177589/143‑1 (MQ=255)                       accCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACt                                                                                                                         >  3:35261/1‑148 (MQ=255)                         ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTa                                                                                                                       >  3:363119/1‑148 (MQ=255)                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                                    <  5:103727/149‑1 (MQ=255)                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                              >  6:366033/1‑123 (MQ=255)                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                              >  6:371043/1‑123 (MQ=255)                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                              <  5:371043/123‑1 (MQ=255)                           acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                              <  5:366033/123‑1 (MQ=255)                                 tCTGCTTGCCGGTCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaa                                                                                                              <  4:35261/149‑1 (MQ=255)                                     cTTGCCGGGCGATACCCCCAATGGTTGATCAGCGAAGCAGAGATGCCTCAGAGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaa                                                                                                              <  1:304256/147‑1 (MQ=38)                                     cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTg                                                                                                          >  2:43994/1‑149 (MQ=255)                                        gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                        <  4:363119/148‑1 (MQ=255)                                             gCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                          >  1:201224/1‑109 (MQ=255)                                             gCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                          <  2:201224/109‑1 (MQ=255)                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                                 >  2:131392/1‑95 (MQ=255)                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                                 <  1:131392/95‑1 (MQ=255)                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaaa                                                                                                             >  5:190054/1‑131 (MQ=255)                                                    ccccAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaaa                                                                                                             <  6:190054/131‑1 (MQ=255)                                                     cccAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                          >  4:32937/1‑149 (MQ=255)                                                           ggTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                            >  6:49452/1‑93 (MQ=255)                                                           ggTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                            <  5:49452/93‑1 (MQ=255)                                                             ttGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                                  >  2:356191/1‑149 (MQ=255)                                                                aTCAGCGAAGCGGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                               >  5:215883/1‑149 (MQ=255)                                                                aTCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                               >  6:167301/1‑149 (MQ=255)                                                                            gagaTGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                                    >  8:61134/1‑148 (MQ=255)                                                                                tGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGCCACCACAAGGCTACGCGGATCGCCCACCTGGa                                                               <  1:43994/149‑1 (MQ=255)                                                                                           tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccg                                                      >  7:307448/1‑149 (MQ=255)                                                                                                          aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGTCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGTATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                       <  3:32937/149‑1 (MQ=255)                                                                                                          aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                       >  3:15962/1‑149 (MQ=255)                                                                                                              aTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTg                                     >  2:191736/1‑147 (MQ=255)                                                                                                                  aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                               >  5:252995/1‑149 (MQ=255)                                                                                                                        ccTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGCGGCAGGGACAGGCAACTTAAACTCAGCAACTGTAGCTCCCCGTCCCCCCACGGCTCCGCGGAGCGCCCCCCTGGCTCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtgg                         >  3:211722/1‑147 (MQ=37)                                                                                                                            tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                      <  6:215883/148‑1 (MQ=255)                                                                                                                            tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                      <  5:167301/148‑1 (MQ=255)                                                                                                                               cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                                <  8:103980/119‑1 (MQ=255)                                                                                                                               cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                                >  7:103980/1‑119 (MQ=255)                                                                                                                               cTGCCGCGTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAAGCAGAAACTGTAGCTACCCGGCCCCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGGCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                  <  4:211722/149‑1 (MQ=255)                                                                                                                                tGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt                 >  5:189411/1‑149 (MQ=255)                                                                                                                                   ctctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt              <  6:189411/149‑1 (MQ=255)                                                                                                                                   ctctTCGCGGGCATGCCCGCGAAGAGGCAGGTACAGGCACCTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCAACCTGGATCTGAGCCGCACTGGTCGCCGAGCTGTTGTAGGGCTTGTCAGCAATGtt              >  2:53314/1‑149 (MQ=255)                                                                                                                                    tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg             >  4:91170/1‑149 (MQ=255)                                                                                                                                    tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                  >  5:92920/1‑144 (MQ=255)                                                                                                                                    tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                  <  6:92920/144‑1 (MQ=255)                                                                                                                                    tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg             >  6:335461/1‑149 (MQ=255)                                                                                                                                    tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg             >  2:197695/1‑149 (MQ=255)                                                                                                                                         gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGt               <  2:200646/142‑1 (MQ=255)                                                                                                                                         gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGt               >  1:200646/1‑142 (MQ=255)                                                                                                                                          cgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt       <  1:356191/149‑1 (MQ=255)                                                                                                                                             gCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgcc    <  8:307448/149‑1 (MQ=255)                                                                                                                                               aTGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgccgc  <  7:61134/149‑1 (MQ=255)                                                                                                                                               aTGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgccgc  <  1:197695/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 AGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGC  >  NC_002947/321996‑322284

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