Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,469,817	C→T	12.1%	R109H (CGC→CAC)	algK ←	alginate biosynthesis protein AlgK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,469,817	0	C	T	12.1%	129.9 / 5.7	49	R109H (CGC→CAC)	algK	alginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base C (25/18);  new base T (3/3);  total (28/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 = 1.00e+00

GGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTC  >  NC_002947/1469677‑1469956                                                                                                                                             |                                                                                                                                             ggAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAgg                                                                                                                                      >  1:21100/1‑149 (MQ=255)  gAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGgc                                                                                                                                     >  2:92828/1‑149 (MQ=255)         tGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTg                                                                                                                             >  1:37889/1‑149 (MQ=255)            gggTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTGGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑Ctgcggc                                                                                                                          >  1:78414/1‑144 (MQ=255)                  cTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGGGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTgg                                                                                                                    <  5:260144/149‑1 (MQ=255)                     aGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGc                                                                                                                   <  4:129683/147‑1 (MQ=255)                      ggTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGAc                                                                                                                >  7:246162/1‑149 (MQ=255)                               ttGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGcc                                                                                                       >  6:278526/1‑149 (MQ=255)                                 gCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                      <  1:244188/148‑1 (MQ=255)                                  ccAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGc                                                                                                     <  2:37889/148‑1 (MQ=255)                                   cAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCACGCGCCCAA‑GCGCGCCTTGGCCGCCTg                                                                                                                     >  1:325119/1‑131 (MQ=255)                                   cAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCACGCGCCCAA‑GCGCGCCTTGGCCGCCTg                                                                                                                     <  2:325119/131‑1 (MQ=255)                                                   ttGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                   <  3:329701/149‑1 (MQ=255)                                                      cTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCgg                                                                                 <  1:150056/148‑1 (MQ=255)                                                      cTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGc                                                                                >  3:244512/1‑149 (MQ=255)                                                         tGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                   <  4:198703/143‑1 (MQ=255)                                                         tGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTc                                                                                   >  3:198703/1‑143 (MQ=255)                                                            cccTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCtt                                                                          >  3:325912/1‑149 (MQ=255)                                                              cTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGa                                                                        >  1:345946/1‑149 (MQ=255)                                                                 tttCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGAttt                                                                     >  4:308250/1‑149 (MQ=255)                                                                     gcgGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGc                                                                   >  8:148082/1‑147 (MQ=255)                                                                           tGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTc                                                           >  3:68180/1‑149 (MQ=255)                                                                                gAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCgg                                                       <  3:320392/148‑1 (MQ=255)                                                                                 aGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGt                                                     <  7:148082/149‑1 (MQ=255)                                                                                  gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                    <  8:246162/149‑1 (MQ=255)                                                                                    aGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGc                                                  >  1:373599/1‑149 (MQ=255)                                                                                           cGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCGCCTGCTGGTGGCCAAgc                                                                                                                                      <  6:123891/59‑2 (MQ=25)                                                                                           cGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCGCCTGCTGGTGGCCAAgc                                                                                                                                      >  5:123891/1‑58 (MQ=25)                                                                                             gCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGAt                                         <  3:308250/149‑1 (MQ=255)                                                                                                  cTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                     <  5:278526/148‑1 (MQ=255)                                                                                                          cGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGcc                             >  8:99475/1‑148 (MQ=255)                                                                                                          cGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGcc                             <  2:373599/148‑1 (MQ=255)                                                                                                            gCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                     <  3:130215/138‑1 (MQ=255)                                                                                                            gCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTcc                                     >  4:130215/1‑138 (MQ=255)                                                                                                            gCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGAc                          >  1:101983/1‑149 (MQ=255)                                                                                                            gCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGAc                          >  4:375430/1‑149 (MQ=255)                                                                                                             ccTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACt                         >  5:325423/1‑149 (MQ=255)                                                                                                                     cTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCtt                  <  2:21100/148‑1 (MQ=255)                                                                                                                        tCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGt                                                     >  2:172749/1‑110 (MQ=255)                                                                                                                        tCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGt                                                     <  1:172749/110‑1 (MQ=255)                                                                                                                          cGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTg             >  4:5622/1‑148 (MQ=255)                                                                                                                           ggCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTa           <  4:244512/149‑1 (MQ=255)                                                                                                                           gccTTGGCCGCCTGCTGGTGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCg                                                                                  >  5:303232/3‑78 (MQ=37)                                                                                                                           gccTTGGCCGCCTGCTGGTGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCg                                                                                  <  6:303232/76‑1 (MQ=37)                                                                                                                             cTTGGCCGCCTGCTGGTGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                              <  2:377157/128‑1 (MQ=255)                                                                                                                             cTTGGCCGCCTGCTGGTGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                              >  1:377157/1‑128 (MQ=255)                                                                                                                                  ccaccaGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAgg    >  6:381826/1‑149 (MQ=255)                                                                                                                                   caccaGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGt   >  4:284562/1‑149 (MQ=255)                                                                                                                                    accaGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTc  >  1:120363/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                             GGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTCCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCG‑CTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTC  >  NC_002947/1469677‑1469956

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