Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,464,200	T→C	9.5%	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.5%	106.2 / 2.4	42	Q473R (CAG→CGG)	algX	alginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (23/15);  new base C (2/2);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

TTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGG  >  NC_002947/1464059‑1464340                                                                                                                                                |                                                                                                                                             ttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAgg                                                                                                                                       >  6:181080/1‑149 (MQ=255)    gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                     <  3:21730/148‑1 (MQ=255)     aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                   <  4:188604/149‑1 (MQ=255)     aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                   <  8:15366/149‑1 (MQ=255)         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGAGGGGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                 >  1:114628/1‑149 (MQ=255)         ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCg                                                                                                                                >  4:264811/1‑148 (MQ=255)           tGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                             >  4:274885/1‑149 (MQ=255)                   aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGACTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCg                                                                                                                     >  6:337840/1‑149 (MQ=255)                   aCCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAAATAGATAGAAAAGCGTCATGG‑‑GTCGCCTCACTGTACTGACTGCCCGGTCTGAGAAGGCACGGCGAACTTTTTTCg                                                                                                                     >  6:255930/1‑149 (MQ=255)                    cccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAATTTGTTACGt                                                                                                                    >  8:185721/1‑149 (MQ=255)                      cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                  <  8:244232/149‑1 (MQ=255)                        gcggcgGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                                                >  1:161051/1‑149 (MQ=255)                               tACCAGCGCGGAGGGGGCGGCACCACCAAACAATCCCAACGCGTGCTGGGGGGGGGTGTATATTTGTTTTAAACGGGTCATGG‑‑GTCGCCTCACTGGAGTGCGTGGGGGGTATGCGAAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                         <  6:328412/149‑1 (MQ=255)                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                        <  7:277587/145‑1 (MQ=255)                                    gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                        >  8:277587/1‑145 (MQ=255)                                                 ggCACCACCAAACACGGCCACCGCGAGCCGGGGGGGGCTGTAAATTCGTTTGAAACGGGTCAGGG‑‑GGCGCCTCACTGGCCGGCCTGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                       <  3:274885/149‑1 (MQ=255)                                                    accaccAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGcc                                                                                    >  3:213617/1‑149 (MQ=255)                                                            aCAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCg                                                                            <  3:264811/149‑1 (MQ=255)                                                                  ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCt                                                                      <  8:222024/149‑1 (MQ=255)                                                                       gcgAGCAGGGCGGGGCTGCAAATGCGTTTGAAAAGGGGCATGG‑‑GGCGCCTCCCTGGCCGGCCTCCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTAcc                                                                 <  5:181080/149‑1 (MQ=255)                                                                         gAGCTGGGCGGGGCTGGTAATTCGTTTGAAAAGGGTCTTGG‑‑GTCGCCTAAGTGGCCTGCGTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                               <  8:82248/149‑1 (MQ=255)                                                                         gAGCCGGGCGGGGCTGGAAATTCGTTTTAAACGGGTCATGG‑‑GGCGCCTCACTGCCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCAACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTc                                                               <  4:253838/149‑1 (MQ=255)                                                                           gCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCg                                                              >  5:125242/1‑148 (MQ=255)                                                                                     gCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGCCCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAg                                                   >  6:253827/1‑149 (MQ=255)                                                                                                gTTTGAAACGGTTAATGG‑‑GTCGCCTCAGTGGCCTGCGTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGctggctgg                                        <  7:185721/149‑1 (MQ=255)                                                                                                 tttGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGAGCGGTATGCGCAGGCAAAGCGCAGTTGTTTAGTTTTCATAGCTTGTCTTAGACAGTAAGGGTGCCGCATGACCGCCCCCATACCTAGAAACCGAGACAGCTCtagctggc                                       >  8:296534/1‑149 (MQ=255)                                                                                                    gAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCa                                                    >  1:79645/1‑133 (MQ=255)                                                                                                    gAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCa                                                    <  2:79645/133‑1 (MQ=255)                                                                                                            tCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCg                             >  4:227926/1‑148 (MQ=255)                                                                                                            tCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCGTCGACCTTCTGGGTGCCGCTTTCCGGCCCCGGTACCTCGAACGCGCGCAGGCGCTGGCTGGCCCCGTCCTCGt                            >  6:304995/1‑149 (MQ=255)                                                                                                              aTGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCCGGCACGGCGAAGTTGTTGCGGTTGCAGAGCTTGGCTTCGACCATCTGGGTGCCGCTTTCCGGCCCCGGTACCTAGAACGCGACCAGGCGCTGGCAGGCCCAGTCCTCGTCt                          >  4:189120/1‑149 (MQ=255)                                                                                                                         cTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATc                 >  7:80919/1‑149 (MQ=255)                                                                                                                             cTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg              >  5:90337/1‑148 (MQ=255)                                                                                                                              tGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTCTACCTCGAGCGCGAGCAGGCGCTGGCTGGCCCAgtcgtc                              >  1:55306/1‑131 (MQ=255)                                                                                                                                gcctgcctGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGa          >  1:187654/1‑149 (MQ=255)                                                                                                                                gcctgcctGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGa          >  3:45763/1‑149 (MQ=255)                                                                                                                                  gtgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc         >  7:126401/2‑148 (MQ=255)                                                                                                                                  ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc         >  3:268005/1‑148 (MQ=255)                                                                                                                                    gctggggcgGTGCGCGCTGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa      <  3:189120/143‑1 (MQ=255)                                                                                                                                    gcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACg        >  4:131959/1‑147 (MQ=255)                                                                                                                                      ctgcGCGGTCTGCGAAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAg     <  4:213617/148‑1 (MQ=255)                                                                                                                                         cgcgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCgg  <  6:125242/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                             TTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGG‑‑GTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGG  >  NC_002947/1464059‑1464340

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