Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,759,457	T→G	10.0%	P109P (CCA→CCC)	pvdP ←	pyoverdine biosynthesis‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,759,457	0	T	G	10.0%	139.2 / 7.2	50	P109P (CCA→CCC)	pvdP	pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base G (3/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

GCCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGC  >  NC_002947/4759309‑4759601                                                                                                                                                      |                                                                                                                                                  gccCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGt                                                                                                                                                   <  1:152150/149‑1 (MQ=255) gccCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGACATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGt                                                                                                                                                   <  1:49383/149‑1 (MQ=255)  ccccGCATTGCCAAATACCCCCCACACTTAATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTg                                                                                                                                                  >  3:34560/1‑149 (MQ=255)                      cccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTg‑g                                                                                                                                                >  6:211187/1‑130 (MQ=255)                      cccACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTg‑g                                                                                                                                                <  5:211187/130‑1 (MQ=255)                                attcAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTcggcg                                                                                                                    <  4:34560/145‑1 (MQ=255)                                aTCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTcggcgg                                                                                                                   >  5:204355/1‑149 (MQ=255)                                  ccAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTa                                                                                                                 >  5:162707/1‑149 (MQ=255)                                      cACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTGGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATc                                                                                                               >  2:214475/1‑147 (MQ=255)                                       aCGAGATGGCACTCGTAGCGCGGCATGTGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTcggcg                                                                                                                    <  1:100519/141‑1 (MQ=255)                                       aCGAGATGGCACTCGTAGCGCGGCATGTGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTcggcg                                                                                                                    >  2:100519/1‑141 (MQ=255)                                       aCGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAg                                                                                                            <  6:21384/149‑1 (MQ=255)                                       aCGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAg                                                                                                            <  6:204355/149‑1 (MQ=255)                                         gagaTGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAgg                                                                                                           >  3:184475/1‑148 (MQ=255)                                           gaTGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTaa                                                                                                        >  6:103510/1‑149 (MQ=255)                                              gggaCTCGTAGCGCGGCATTCGGGCGCCATGTGGACTGTCGAACATGCGCCAGCGTTACCGTCCGGGCTTGT‑CGGCAGGCAAGTCGCCGAGCAGCCGTTATCGTG‑GCTCTGCACGGCTGCGCAAGCGGTCTGCGGTATCGAGTTatgc                                                                                                      <  8:84132/145‑4 (MQ=255)                                                     gTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAAt                                                                                                <  4:184475/147‑1 (MQ=255)                                                        gcgcGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAAtgctgc                                                                                           >  4:325082/1‑149 (MQ=255)                                                         cgcgGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCg                                                                                          >  1:232683/1‑149 (MQ=255)                                                          gcgGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCgg                                                                                         <  6:162707/149‑1 (MQ=255)                                                             gCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGccc                                                                                      >  1:181602/1‑149 (MQ=255)                                                                     ccgccATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTgcgc                                                                              <  3:325082/149‑1 (MQ=255)                                                                        ccATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGC‑CACGATATCGGGTGCTCgg                                                                                                                                           >  1:218185/1‑85 (MQ=37)                                                                        ccATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGC‑CACGATATCGGGTGCTCgg                                                                                                                                           <  2:218185/85‑1 (MQ=37)                                                                              gCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAt                                                                     <  6:333665/149‑1 (MQ=255)                                                                              gCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAt                                                                     <  5:103510/149‑1 (MQ=255)                                                                                   gTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGTCCGGCAGGCAAGTCGCCGAGCACCCGATATCGTg                                                                                                                                                  <  1:43028/69‑1 (MQ=255)                                                                                    tCGAACAAGCGCCAGCGTAACCATCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACaa                                                               >  1:244180/1‑149 (MQ=255)                                                                                      gAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAAtt                                                             >  3:112295/1‑149 (MQ=255)                                                                                         cAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTc                                                            >  2:235998/1‑147 (MQ=255)                                                                                            gcgcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCACAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCACAACATCGATGAACAATTCCAGCg                                                       >  7:127555/1‑149 (MQ=255)                                                                                             cgcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGt                                                      <  8:247526/149‑1 (MQ=255)                                                                                              gcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCa                                                                             <  8:225857/125‑1 (MQ=255)                                                                                              gcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCa                                                                             >  7:225857/1‑125 (MQ=255)                                                                                              gcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAt                                                                     <  4:268456/133‑1 (MQ=255)                                                                                              gcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGAt                                                                     >  3:268456/1‑133 (MQ=255)                                                                                              gcCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGtc                                                     <  4:112295/149‑1 (MQ=255)                                                                                                cAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGtctc                                                   >  6:225196/1‑149 (MQ=255)                                                                                                    gTACCCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATCTCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGccc                                               <  8:76740/149‑1 (MQ=255)                                                                                                       aCCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGc                                            <  4:75888/149‑1 (MQ=255)                                                                                                         cGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGcc                                                <  1:121552/143‑1 (MQ=255)                                                                                                         cGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGcc                                                >  2:121552/1‑143 (MQ=255)                                                                                                         cGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCaa                                          <  2:232683/149‑1 (MQ=255)                                                                                                                  ttGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGc                                   >  8:60541/1‑147 (MQ=255)                                                                                                                     t‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGccc                                      >  7:295613/1‑141 (MQ=255)                                                                                                                     t‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGccc                                      <  8:295613/141‑1 (MQ=255)                                                                                                                                  gtgCCGAGC‑CACGATATCGGGTGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGcccc                  >  5:225193/3‑149 (MQ=255)                                                                                                                                   tgCCGAGC‑CACGATATCGGGTGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGccccc                 >  5:138571/2‑149 (MQ=255)                                                                                                                                    gCCGAGCACCCGATATCGTG‑GCTCGGCACAGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGccccat                >  7:54573/1‑147 (MQ=255)                                                                                                                                     ccGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCTCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCtt               <  1:235998/149‑1 (MQ=255)                                                                                                                                            cacGATATCGGGTGCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGa         <  6:138571/147‑1 (MQ=255)                                                                                                                                                     gTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAgcgc  >  2:10483/1‑146 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GCCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGC  >  NC_002947/4759309‑4759601

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