Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,759,458:1	+T	10.9%	coding (326/1623 nt)	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,458	1	.	T	10.9%	162.0 / 15.7	46	coding (326/1623 nt)	pvdP	pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base . (21/20);  new base T (3/2);  total (24/22)
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.19e-01

CCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGC  >  NC_002947/4759310‑4759601                                                                                                                                                       |                                                                                                                                                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                                                                                                                                           <  2:218185/85‑1 (MQ=37)                                                                       ccATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGC‑CACGATATCGGGTGCTCgg                                                                                                                                           >  1:218185/1‑85 (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                                                >  2:121552/1‑143 (MQ=255)                                                                                                        cGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGcc                                                <  1:121552/143‑1 (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)                                                                                                                                                       |                                                                                                                                                CCCCGCATTGCCAAATACCCCCCACACTTCATCCAGCACGAGATGGCACTCGTAGCGCGGCATGCGGCCGCCATCGGCACTGTCGAACAAGCGCCAGCGTAACCGTCCCGGCTTGT‑CGGCAGGCAAGTCGCCGAGCACCCGATATCGTG‑GCTCGGCACCGCTGCGCAAGCGCTCGGCGGTATCGAGGTAACCGCGAATGCTGCGGCCCTTCTGCGCAACATCGATGAACAATTCCAGCGTCTCGCCCTGCAAGCCCGGCAAGCCGGCCTGCGCCCCCTTGAACGACCAGCGC  >  NC_002947/4759310‑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: ‑