Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	322,142	T→C	22.2%	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,142	0	T	C	22.2%	78.3 / 20.4	45	intergenic (+427/+31)	aguA/PP_0267	agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (19/16);  new base C (5/5);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

ATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTC  >  NC_002947/322001‑322278                                                                                                                                              |                                                                                                                                         aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  <  2:409705/150‑1 (MQ=255) aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  >  1:512119/1‑150 (MQ=255)   cAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                                <  2:799926/150‑1 (MQ=255)     gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                              <  1:763510/150‑1 (MQ=255)     gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                              <  2:512119/150‑1 (MQ=255)      ggTCAGGACGACCGACCCACGACGTCTGCTTGGCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTGGCGGGCATGCTCGCGAAGGGGCAg                                                                                                                             >  1:469035/1‑150 (MQ=255)            gacgacCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  <  1:11452/139‑1 (MQ=255)            gacgacCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  >  2:11452/1‑139 (MQ=255)               gaccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAgg                                                                                                                                <  2:95326/138‑1 (MQ=255)               gaccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAgg                                                                                                                                >  1:95326/1‑138 (MQ=255)                accgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCaa                                                                                                                   >  2:357786/1‑150 (MQ=255)                 ccgaccCACGACGTCTGATTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  <  2:226485/134‑1 (MQ=255)                 ccgaccCACGACGTCTGATTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                  >  1:226485/1‑134 (MQ=255)                                   ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                        <  1:526287/126‑1 (MQ=255)                                   ttGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACa                                                                                                                        >  2:526287/1‑126 (MQ=255)                                         ggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAg                                                                                                                             <  2:543444/115‑1 (MQ=255)                                         ggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAg                                                                                                                             >  1:543444/1‑115 (MQ=255)                                                  cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcaccac                                                                                 >  2:638687/1‑150 (MQ=255)                                                  cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAgg                                                                                                                                <  1:245533/103‑1 (MQ=255)                                                  cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAgg                                                                                                                                >  2:245533/1‑103 (MQ=255)                                                         gTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                          >  1:793604/1‑150 (MQ=255)                                                         gTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                          <  1:34760/150‑1 (MQ=255)                                                          ttGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAc                                                                         <  2:796399/150‑1 (MQ=255)                                                              tCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                             <  1:352323/94‑1 (MQ=255)                                                              tCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                             >  2:352323/1‑94 (MQ=255)                                                                             tGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGAt                                                      <  1:357786/150‑1 (MQ=255)                                                                                      tGATGGCGGTGTGGCACTGATCTATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccgc                                             >  2:105847/1‑150 (MQ=255)                                                                                         tGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                       <  2:347187/73‑1 (MQ=255)                                                                                         tGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAg                                                                                                                       >  1:347187/1‑73 (MQ=255)                                                                                                          tCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtagt                         >  2:435837/1‑150 (MQ=255)                                                                                                              ttGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                           >  1:11177/1‑80 (MQ=255)                                                                                                              ttGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                           <  2:11177/80‑1 (MQ=255)                                                                                                                 cACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtt                  >  2:376804/1‑150 (MQ=255)                                                                                                                       tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc               >  1:68324/1‑147 (MQ=255)                                                                                                                       tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc               <  2:68324/147‑1 (MQ=255)                                                                                                                       tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGc            >  1:461759/1‑150 (MQ=255)                                                                                                                         ccTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa          >  1:533950/1‑150 (MQ=255)                                                                                                                          cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt         >  2:202351/1‑150 (MQ=255)                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa              <  2:213335/140‑1 (MQ=255)                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa              >  1:213335/1‑140 (MQ=255)                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                       >  1:728324/1‑131 (MQ=255)                                                                                                                               tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                       <  2:728324/131‑1 (MQ=255)                                                                                                                                 ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCGCCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTCGTCAGCAATGTTGCTc  <  1:105847/150‑1 (MQ=255)                                                                                                                                        gCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc    >  1:833738/1‑141 (MQ=255)                                                                                                                                        gCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc    <  2:833738/141‑1 (MQ=255)                                                                                                                                              |                                                                                                                                         ATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTC  >  NC_002947/322001‑322278

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