Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,464,200	T→C	9.7%	Q473R (CAG→CGG)	algX ←	alginate biosynthesis protein AlgX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,464,200	0	T	C	9.7%	130.3 / 5.2	53	Q473R (CAG→CGG)	algX	alginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (30/16);  new base C (2/3);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

CTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464056‑1464346                                                                                                                                                   |                                                                                                                                                   cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTgcgc                                                                                                                                                <  7:129098/149‑1 (MQ=255) cTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTgcg                                                                                                                                                 >  8:129098/1‑148 (MQ=255)    ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                              <  6:81444/148‑1 (MQ=255)     tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCAGACCGCGCAgg                                                                                                                                             <  1:347325/148‑1 (MQ=255)      cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCa                                                                                                                                               >  3:13256/1‑145 (MQ=255)      cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCa                                                                                                                                               <  4:13256/145‑1 (MQ=255)      cGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCa                                                                                                                                           >  7:204552/1‑149 (MQ=255)       gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          <  3:91652/149‑2 (MQ=255)         tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCCCCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTCAAAAGGGTCAGGGCTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                        >  5:277315/1‑149 (MQ=255)         tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACgg                                                                                                                                        >  8:190642/1‑149 (MQ=255)             cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    <  3:231430/149‑1 (MQ=255)                gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCCGGCACGGCGAAgtt                                                                                                                                 >  1:348764/1‑149 (MQ=255)                gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgtt                                                                                                                                 >  8:198812/1‑149 (MQ=255)                gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                           >  1:151490/1‑139 (MQ=255)                gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                           <  2:151490/139‑1 (MQ=255)                        ccTGCGGCGGTACCCGCGGGGGGTGCGGGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGTAATTTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    <  3:238845/138‑1 (MQ=255)                         cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACATGATGCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGt                                                                                                                          >  1:123845/1‑149 (MQ=255)                          ggcggggGCACCTGGGGGGCGGGCGGGCAACCACAAAACACGCCCACCGCGCGCGGGGGGGGGCTGTAATTTCGTTTCAAACGGGTCTGGGGCCGCCTCAC‑GG‑GCCTGCCTGGGGGGTCGGCGCAGCCACGGCGAAGTTGTTGCGTTTg                                                                                                                       <  8:204552/148‑1 (MQ=255)                            cgggggAACCTGCGCGGCGTGCGCGGCACCACCAAACAAGGCCACCGCGACCAGGGCGGGGCTGGAAATTCGTTTTAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                       <  7:327698/147‑1 (MQ=255)                                    ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                                             >  3:290879/1‑149 (MQ=255)                                    ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                                             >  3:285772/1‑149 (MQ=255)                                        cgcgGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGCGTCGCCTCAC‑TG‑GCCTGCCCGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                         <  8:142005/149‑1 (MQ=255)                                            gcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTgc                                                                                                                                                  >  4:114819/1‑104 (MQ=255)                                            gcgTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTgc                                                                                                                                                  <  3:114819/104‑1 (MQ=255)                                                    ggCACCACCAACCACGCCCAACGCTCGCTGGGGGGGGCTGGAAATTCGTTTCAAACGGGCCTTGGGGCGCCTCAG‑GG‑GCCGGCCGGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                             <  7:198812/149‑1 (MQ=255)                                                         caccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATGAGTTGGAAAAGGGTCATGGGTCGCCTCAA‑TG‑GCATGCCTGGGCGGTCGGCGAGGGCACGGCGATGTTGTTGCGTTTGCAGCGCTTGGCTTGGCCGTTCTGGTTGCCGc                                                                                        >  6:206568/1‑149 (MQ=255)                                                          accaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCt                                                                                       >  3:53051/1‑149 (MQ=255)                                                           ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCTTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                              <  2:61272/109‑1 (MQ=255)                                                            caAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                       >  3:208563/1‑115 (MQ=255)                                                             aaaCAAGGCCAACGCGACCTGGGCGGTGCTGTAATTTGGTTTGAAAAGGGCCATGGGTGGCCTCAC‑TT‑GCGTGCCTGCGGGGTTTGCGATGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTc                                                                                    <  6:64519/149‑1 (MQ=255)                                                                 aaGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGGTTGAAAAGGGTCAAGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGCCCTTCTGGGTGCCGCTTTCCGGc                                                                                >  5:270134/1‑149 (MQ=255)                                                                           cgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               >  7:94100/1‑124 (MQ=255)                                                                                        gCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             >  8:264863/1‑145 (MQ=255)                                                                                        gCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             <  7:264863/145‑1 (MQ=255)                                                                                            gAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               >  7:243263/1‑107 (MQ=255)                                                                                                 tcttttGAAACGGGCCTTGGGGCCCCTCAG‑TG‑GCCTGCCTGCGGGGTCTGCGCTGCCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggc                                                 <  4:53051/145‑1 (MQ=255)                                                                                                         aaaGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGt                                        >  6:340060/1‑149 (MQ=255)                                                                                                           aGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                                      <  7:190642/149‑1 (MQ=255)                                                                                                            gggTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                        >  7:104298/1‑130 (MQ=255)                                                                                                               tCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                          >  5:300532/1‑141 (MQ=255)                                                                                                                 aTGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCt                                >  7:53289/1‑149 (MQ=255)                                                                                                                     gTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                            <  5:340060/149‑1 (MQ=255)                                                                                                                          cacAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTAGCGCATc                       >  3:41665/3‑149 (MQ=255)                                                                                                                          cTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTccagatgc                       >  8:181769/1‑142 (MQ=255)                                                                                                                             aC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTg                     >  6:264534/1‑148 (MQ=255)                                                                                                                              c‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg                    >  7:287710/1‑148 (MQ=255)                                                                                                                                    cgtgcctgcGCAGACCGAGCACGCA‑GGCCAAGTTGTTGAGTTTGGAGAGCTTGGATGCGAACTTCTGGGTGCCGCTGTACGGCCCCTGTACCGCGACCGCCCGCCGGCGCTGGCAGGCCCAGTCCTCGTCTTCGCGCATCTGGatgcc               >  8:264074/3‑144 (MQ=255)                                                                                                                                     ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCAGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc               >  5:204243/1‑148 (MQ=255)                                                                                                                                      tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACg              >  3:107997/1‑148 (MQ=255)                                                                                                                                            cgcgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGc       >  4:283549/1‑149 (MQ=255)                                                                                                                                              cgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCg     <  8:287710/149‑1 (MQ=255)                                                                                                                                                 tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  >  7:216603/1‑149 (MQ=255)                                                                                                                                                 tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  >  5:29350/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCAC‑TG‑GCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464056‑1464346

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