Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,469,811	A→T	11.6%	L111Q (CTG→CAG)	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,811	0	A	T	11.6%	112.3 / 6.1	43	L111Q (CTG→CAG)	algK	alginate biosynthesis protein AlgK
Reads supporting (aligned to +/- strand):  ref base A (26/12);  new base T (2/3);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGC  >  NC_002947/1469672‑1469959                                                                                                                                             |                                                                                                                                                     cTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAggcggc                                                                                                                                             >  2:263179/1‑149 (MQ=255)    ggggAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGcc                                                                                                                                            >  4:293284/1‑147 (MQ=255)         aaCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGgcgc                                                                                                                                     >  7:236267/1‑149 (MQ=255)         aaCTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGgcgc                                                                                                                                     <  2:309332/149‑1 (MQ=255)              tGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGCCTTGGCCGCCTGCTGGTGGCCAAGccggaca                                                                                                                                  <  3:109607/147‑5 (MQ=255)                 gggTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGgc                                                                                                                             >  4:279406/1‑149 (MQ=255)                       cTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTgg                                                                                                                       <  3:293284/149‑1 (MQ=255)                           ggTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGAc                                                                                                                   <  8:246198/149‑1 (MQ=255)                               gagCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTgg                                                                                                               >  7:69869/1‑149 (MQ=255)                                 gCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGTACCGCCTTGGCCGCCTGCTGGTGGCCAAGccg                                                                                                                                      <  1:32332/125‑1 (MQ=255)                                 gCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGCCTTGGCCGCCTGCTGGTGGCCAAGccg                                                                                                                                      >  2:32332/1‑124 (MQ=39)                                      gCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCg                                                                                                         >  1:65769/1‑148 (MQ=255)                                           cGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAgcgc                                                                                                   <  2:59291/149‑1 (MQ=255)                                                tCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCa                                                                                                       >  7:171199/1‑140 (MQ=255)                                                tCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCa                                                                                                       <  8:171199/140‑1 (MQ=255)                                                tCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAg                                                                                              >  1:127392/1‑149 (MQ=255)                                                        ttGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGcc                                                                                                                                            <  5:241957/95‑1 (MQ=255)                                                        ttGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGcc                                                                                                                                            >  6:241957/1‑95 (MQ=255)                                                                gCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTACTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCt                                                                              >  6:148999/1‑149 (MQ=255)                                                                           cgGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAg                                                                   >  6:1356/1‑149 (MQ=255)                                                                                tGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTc                                                              >  1:273842/1‑149 (MQ=255)                                                                                 gCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcg                                                             >  5:307965/1‑149 (MQ=255)                                                                                  cTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTcgc                                                            <  2:65769/149‑1 (MQ=255)                                                                                      aGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGt                                                        >  3:291085/1‑149 (MQ=255)                                                                                       gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                       <  2:127392/149‑1 (MQ=255)                                                                                       gCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTg                                                       >  1:120062/1‑149 (MQ=255)                                                                                        cAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGt                                                        >  1:109126/1‑147 (MQ=255)                                                                                         aGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGc                                                     >  5:26688/1‑149 (MQ=255)                                                                                                       cTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGc                                                                                   <  5:102583/105‑1 (MQ=255)                                                                                                       cTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGc                                                                                   >  6:102583/1‑105 (MQ=255)                                                                                                          aCGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGc                                      >  1:189289/1‑147 (MQ=255)                                                                                                          aCGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGc                                      >  8:104736/1‑147 (MQ=255)                                                                                                             ttCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                                 >  1:186383/1‑149 (MQ=255)                                                                                                             ttCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGc                                 <  3:279406/149‑1 (MQ=255)                                                                                                                ggCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGa                              <  7:104736/149‑1 (MQ=255)                                                                                                                    tGGGTGCTGT‑CCGGCTTGGCCCCCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTg                          >  3:193050/1‑149 (MQ=255)                                                                                                                          cTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTc                    <  8:69869/149‑1 (MQ=255)                                                                                                                            gT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGc                  >  6:49810/1‑149 (MQ=255)                                                                                                                                cGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTg                >  2:36001/1‑148 (MQ=255)                                                                                                                                   cTTGGCCGCCTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGc            >  8:72388/1‑149 (MQ=255)                                                                                                                                     tGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGATGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGcc           >  6:161633/1‑148 (MQ=255)                                                                                                                                       gccgccTGCTGGTGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAg        <  7:72388/149‑1 (MQ=255)                                                                                                                                       gCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAg        >  7:196166/1‑149 (MQ=255)                                                                                                                                           ccaGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCg    <  6:26688/149‑1 (MQ=255)                                                                                                                                             agcagGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGc  >  7:37873/1‑149 (MQ=255)                                                                                                                                             agcagGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGc  >  4:22216/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                     CTTGGGGAAACTCTGCGGGTAGCTGAGGTAGAGCATTGCCAGCGGAATCAGGGTATTGCTTTGGCCCTGTTTCGCGGCCTGCTTGAGCAGGGTTTCGGCTTCCTCACGTTCGGCCTGGGTGCTGT‑CCGGCTTGGCCACCAGCAGGCGGCCAAGGCGCGCTTGGGCGCGTGGCGACGTGGCGGCCGCAGCGCGATAGGTGGCCTCGGCTTCCTTGATTTGCGAGGGGTCGCGGGTGGCCACCTTGATGTCCGCCAGGCCGACTTGCGCTTCGCTGTAGCCCAGGTCGGC  >  NC_002947/1469672‑1469959

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