Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,464,200	T→C	9.2%	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.2%	113.2 / 2.7	45	Q473R (CAG→CGG)	algX	alginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (26/13);  new base C (2/2);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

GCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464055‑1464346                                                                                                                                                  |                                                                                                                                                   gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                                 >  5:167824/1‑149 (MQ=255) gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCGgcg                                                                                                                                                 >  5:172191/1‑149 (MQ=255) gCTTTTCGATCCCCTGGTAGTAACCATGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATGCGTTTGAAAAGGGTCATGGGTCGCCTCCGTGGCCTGCCTGCGCGGTCTgcg                                                                                                                                                 >  4:59242/1‑149 (MQ=255)        gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGAGCGGCACCACCAAACAAGTCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCa                                                                                                                                           <  2:52208/148‑1 (MQ=255)         aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          <  6:74364/148‑2 (MQ=255)         aTCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          >  5:74364/1‑147 (MQ=255)           ccccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGc                                                                                                                                       >  3:137698/1‑149 (MQ=255)             ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCg                                                                                                                                      >  6:273269/1‑148 (MQ=255)                 gtagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgtt                                                                                                                                 >  2:28245/1‑149 (MQ=255)                     taACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGttgttg                                                                                                                             >  3:138251/1‑149 (MQ=255)                         ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAg                                                                                                                                   >  7:281784/1‑139 (MQ=255)                             gggcggTACCAGGGGGGGGTGCGCGCAACCACAAAACATGCCCAACGCTAGCAGGGCGGGGGTGTAAATTCGTTTCAAAAGGGTAATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCa                                                                                                                     <  8:134114/148‑1 (MQ=255)                              ggcgtTACCCGCGCGGCTGGCGCGCAACAACCAAACACGGCCAACGCGAGCAGGGGGGGGGTTTAAATTCTTTTTAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCag                                                                                                                    <  6:167824/149‑1 (MQ=255)                              ggcggAACCCGGGCGGCGTGCGCGGAACCACAAAACAAGCCCACCGCGCGCCGGGGGGGGCTGTAAATTCGTTTTAAACGGGTAATGGGGCGCCTCACGGGCGTGCCTGCGCGGTCTGCGCTGGCACGGCGAAGTTGTTGCGTTTGCag                                                                                                                    <  6:172191/149‑1 (MQ=255)                                      cAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                                                            >  4:127303/1‑149 (MQ=255)                                            ggcgGGGGCGGAACCACAAAACATGGCCACCCGGTGGTGGGCGGGGGTGTAATTTGTTTTTAAATGGGTTTTGGGGGGCCTCACGGGCCTGCGGGGGGGGTGGGGCCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGAc                                                                                                      <  4:137698/149‑1 (MQ=255)                                                       caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgtt                                                                                                                              >  5:326633/1‑114 (MQ=255)                                                            ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCCCTGGTCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGCCCTTCTGGGTGACGCtt                                                                                      >  7:149259/1‑149 (MQ=255)                                                            ccaAACAAGGCCAACGCGAGCAGGCCGGGGCTGGAAATTCGTTTGAAAAGCGTCATGGCTCGCCTCACTGGCCTTCCTCCGCGGTCTGCGAAGGCACGGCGAAGTTGATCAGTTTGCAGAGCTTGGCTTCGACCTTCTGCGTGCCGCtt                                                                                      >  7:91922/1‑149 (MQ=255)                                                                  aaGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGc                                                                                >  8:18082/1‑149 (MQ=255)                                                                    ggCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTATGCGCAGGCACGGCGAAGTTTTTTGGTTTTCAGAGCTTGGCTTCGCCCTTCTGGGTGACGCTTTCCGGacc                                                                              >  8:281161/1‑149 (MQ=255)                                                                                                 ttCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               >  3:121830/1‑103 (MQ=255)                                                                                                       tGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTTTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                            <  1:28245/148‑1 (MQ=255)                                                                                                         aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGTGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAg                                         >  7:156777/1‑149 (MQ=255)                                                                                                         aaaaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCa                                          <  7:18082/148‑1 (MQ=255)                                                                                                                    gggTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTc                              >  2:67026/1‑149 (MQ=255)                                                                                                                      gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                       <  5:227226/138‑1 (MQ=255)                                                                                                                      gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTc                                       >  6:227226/1‑138 (MQ=255)                                                                                                                      gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                            >  7:280038/1‑149 (MQ=255)                                                                                                                          ccTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                        >  4:50725/1‑149 (MQ=255)                                                                                                                           cTCACGGCCCGGCGTGCGCGGGGTCCCCCGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATc                       <  6:204932/149‑1 (MQ=255)                                                                                                                            tCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGc                                                           >  1:215194/1‑112 (MQ=255)                                                                                                                            tCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGc                                                           <  2:215194/112‑1 (MQ=255)                                                                                                                              aCTGGCCTGCCTGCGCGGGCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTgg                    <  3:59242/149‑1 (MQ=255)                                                                                                                                    ctgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTACTGGGTGCCGCTTTACTGCACATGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGAc               >  7:200351/1‑148 (MQ=255)                                                                                                                                     tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGaaacctc             >  8:190585/1‑144 (MQ=255)                                                                                                                                     tgcctgcGCAGACCGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCa                                                          >  3:116693/1‑104 (MQ=255)                                                                                                                                       cctgcGCTGTGTGCGCTGGCACGGCGAAGTTGTTGTGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa            <  8:156777/148‑1 (MQ=255)                                                                                                                                         tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCg         >  1:155940/1‑149 (MQ=255)                                                                                                                                          gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                        >  3:160889/1‑101 (MQ=255)                                                                                                                                          gcgcGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAgg                                                        <  4:160889/101‑1 (MQ=255)                                                                                                                                              ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCgg                                                                                 <  3:168701/72‑1 (MQ=255)                                                                                                                                              ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCgg                                                                                 >  4:168701/1‑72 (MQ=255)                                                                                                                                               gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc   <  3:127303/149‑1 (MQ=255)                                                                                                                                                tCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  >  1:181925/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464055‑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: ‑