Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,763,160	C→T	7.4%	S166S (AGC→AGT)	pvdO →	pyoverdine biosynthesis‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,763,160	0	C	T	7.4%	154.3 / 2.5	54	S166S (AGC→AGT)	pvdO	pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base T (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

CCCGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTCGTACACCGCCCCGGTC  >  NC_002947/4763015‑4763307                                                                                                                                                  |                                                                                                                                                     cccGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCg                                                                                                                                                     >  2:18088/1‑147 (MQ=255)  ccGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCg                                                                                                                                                     >  4:89097/1‑146 (MQ=255)  ccGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGaa                                                                                                                                                   <  3:89097/148‑1 (MQ=255)  ccGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAg                                                                                                                                                  >  7:158212/1‑149 (MQ=255)        aGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAgc                                                                                                                                            >  7:47843/1‑149 (MQ=255)          ggCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAgcgc                                                                                                                                          <  7:80946/149‑1 (MQ=255)           gcgcTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAgcgc                                                                                                                                          >  6:116294/1‑148 (MQ=255)           gcgcTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGcgcg                                                                                                                                         <  4:161867/149‑1 (MQ=255)             gcTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGaa                                                                                                                                       <  8:158212/149‑1 (MQ=255)                tCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAatat                                                                                                                                    >  7:57171/1‑149 (MQ=255)                  tGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGc                                                                                                                                  >  3:265730/1‑149 (MQ=255)                   gcgcCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGcc                                                                                                                                 >  5:19844/1‑149 (MQ=255)                        gccAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGTGAAGCCGAGCGCGAATACGCcgcgcg                                                                                                                            >  2:114180/1‑149 (MQ=255)                        gccAGCAAGCCGACCTACAAACAGGGCACGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGTGAAGCCGAGCGCGAATACGCcgcgcg                                                                                                                            >  8:229753/1‑149 (MQ=255)                                   ctacctcAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCaa                                                                                                                <  3:84826/143‑1 (MQ=255)                                    aCCTACAAACAGGGCCCGCTCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAAc                                                                                                               <  5:116294/149‑1 (MQ=255)                                           aaCAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCCATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCg                                                                                                        >  7:249312/1‑149 (MQ=255)                                                    cgcgCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGttccccttc                                                                                               <  4:133265/149‑1 (MQ=255)                                                    cgcgCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGttccccttc                                                                                               <  8:47843/149‑1 (MQ=255)                                                    cgcgCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGttccccttc                                                                                               <  1:232618/149‑1 (MQ=255)                                                    cgcgCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGttccccttc                                                                                               <  8:57171/149‑1 (MQ=255)                                                     gcgcCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCGCCGCCGGTTCAACCGGCCCGttccccttc                                                                                               <  1:18088/149‑1 (MQ=255)                                                       gcCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGttccccttcccc                                                                                            >  8:285298/1‑149 (MQ=255)                                                         cAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTccccttcccctt                                                                                          >  1:266958/1‑149 (MQ=255)                                                          aGCCTGTCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCcccttccccttc                                                                                         <  1:265490/149‑1 (MQ=255)                                                             cTGCCGTGTGCGTGAGCTGGTTCCATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGAt                                                                                      <  8:249312/149‑1 (MQ=255)                                                                        gTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGc                                                                           <  3:224889/149‑1 (MQ=255)                                                                          gAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGa                                                                                                                                        <  6:33071/87‑1 (MQ=255)                                                                          gAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGa                                                                                                                                        >  5:33071/1‑87 (MQ=255)                                                                            gCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTcc                                                                       <  5:105900/149‑1 (MQ=255)                                                                             cTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCa                                                                      <  6:19844/149‑1 (MQ=255)                                                                                gTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCa                                                                      <  2:266958/146‑1 (MQ=255)                                                                                gTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAgg                                                                    >  3:237895/1‑148 (MQ=255)                                                                                          gAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGc                                                         >  5:309574/1‑149 (MQ=255)                                                                                          gAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGc                                                         >  6:268299/1‑149 (MQ=255)                                                                                              aaTACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAg                                                     <  6:309574/149‑1 (MQ=255)                                                                                                tACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCATTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCa                                                   >  1:267419/1‑149 (MQ=255)                                                                                                 aCGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACa                                                                >  6:147086/1‑135 (MQ=255)                                                                                                 aCGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACa                                                                <  5:147086/135‑1 (MQ=255)                                                                                                   gCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCACGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGc                                                >  3:246288/1‑149 (MQ=255)                                                                                                       aGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGTGAAGCCGAGCGCGAATACGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAac                                            <  1:114180/149‑1 (MQ=255)                                                                                                       aGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGTGAAGCCGAGCGCGAATACGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAac                                            <  7:229753/149‑1 (MQ=255)                                                                                                             tGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAaca                                           <  7:267705/144‑1 (MQ=255)                                                                                                             tGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTAc                                      >  2:266856/1‑149 (MQ=255)                                                                                                             tGTCGAAAAAAACCGGCCACAACTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAaca                                           >  8:267705/1‑144 (MQ=255)                                                                                                                 gAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGcc                                  <  4:265730/149‑1 (MQ=255)                                                                                                                                 acCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTc                  <  1:266856/149‑1 (MQ=255)                                                                                                                                 acCTATCGCACGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTc                  <  4:246288/149‑1 (MQ=255)                                                                                                                                   cTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCATCCAGGTACAGCATCGCCCATCACGCAAACACCTACGGCCCCACCCACCGTTAATCg                 >  8:282131/1‑148 (MQ=255)                                                                                                                                            gATAAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTCGTACACCGccc       >  3:302004/1‑149 (MQ=255)                                                                                                                                               cAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAg                                                     >  7:21538/1‑100 (MQ=255)                                                                                                                                               cAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAg                                                     <  8:21538/100‑1 (MQ=255)                                                                                                                                                aGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTCGTACACCGCCCCgg    <  5:268299/148‑1 (MQ=255)                                                                                                                                                 gCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTCGTACACCGCCCCGGTc  >  5:42476/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCCGCCCAGGGCGCTTCTGCGCCGCCAGCAAGCCGACCTACAAACAGGGCCCGCGCCAGCCTGCCGTGTGCGTGAGCTGGTTCGATGTCGAGCAATACGCCAAGTGGCTGTCGAAAAAAACCGGCCACACCTATCGCATGATCAGCGAAGCCGAGCGCGAATATGCCGCGCG‑CGCCGGTTCAACCGGCCCGTTCCCCTTCCCCTTCGATAACGGCGACGCCTCCAGGTACAGCATCGCCCAGCACGCCAACACCTACGGCCCCACCGACGGTTACTCGTACACCGCCCCGGTC  >  NC_002947/4763015‑4763307

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