Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,472,714	G→C	10.6%	S38W (TCG→TGG) ‡	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,472,714	0	G	C	10.6%	132.1 / 7.9	47	S38W (TCG→TGG) ‡	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base G (16/26);  new base C (3/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.07e-01

ACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTG  >  NC_002947/1472575‑1472850                                                                                                                                            |                                                                                                                                          aCCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGT‑CGAAg                                                                                                                                  <  2:87762/148‑1 (MQ=255)   ccATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGata                                                                                                                               >  6:325613/1‑149 (MQ=255)   ccATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGata                                                                                                                               >  4:150090/1‑149 (MQ=255)    cATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGatat                                                                                                                              <  3:105974/149‑1 (MQ=255)             cACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCAGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTgg                                                                                                                       >  6:203037/1‑147 (MQ=255)              aCCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGA‑CGAAGATATCAGCTGGcag                                                                                                                    >  1:308060/1‑149 (MQ=255)               ccTTGCGGCGCAGGTACGGGTCCACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGcagg                                                                                                                   <  3:253532/149‑1 (MQ=255)                cTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGcaggg                                                                                                                  >  4:291950/1‑149 (MQ=255)                     ggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGcagggccagg                                                                                                             <  1:134085/149‑1 (MQ=255)                     ggcgCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGcagggccagg                                                                                                             <  4:269072/149‑1 (MQ=255)                         cAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGa                                                                                                         >  4:92720/1‑149 (MQ=255)                                gggTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCACATGGAGTAGCGCCAGATGCCGACTGCGACGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCAt                                                                                                  >  8:45555/1‑149 (MQ=255)                                  gTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCAt                                                                                                  >  3:205852/1‑147 (MQ=255)                                   tACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAg                                                                                               >  3:40088/1‑149 (MQ=255)                                   tACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAg                                                                                               >  4:274057/1‑149 (MQ=255)                                              gagGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATg                                                                                     <  5:23019/148‑1 (MQ=255)                                              gagGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATg                                                                                     <  4:62642/148‑1 (MQ=255)                                               agGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATg                                                                                     >  2:7043/1‑147 (MQ=255)                                                 gAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGa                                                                                                                                          <  1:42849/93‑1 (MQ=255)                                                 gAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGa                                                                                                                                          >  2:42849/1‑93 (MQ=255)                                                 gAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTa                                                                                   >  3:15814/1‑147 (MQ=255)                                                 gAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGT‑CGAAg                                                                                                                                  <  5:326639/100‑1 (MQ=255)                                                 gAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTCGCAGT‑CGAAg                                                                                                                                  >  6:326639/1‑100 (MQ=255)                                                        aTGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCATCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGcc                                                                           <  5:79668/148‑1 (MQ=255)                                                              cgGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGgc                                                                    <  6:147643/149‑1 (MQ=255)                                                                 aTGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCa                                                                 <  1:7043/149‑1 (MQ=255)                                                                    aaaTGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACt                                                               <  3:92720/148‑1 (MQ=255)                                                                     aaTGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGc                                                             <  4:40088/149‑1 (MQ=255)                                                                                   aTGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTg                                               <  3:291950/149‑1 (MQ=255)                                                                                      gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCTGCGCACTGCAACAGCACTGTCTGGa                                             <  5:203037/148‑1 (MQ=255)                                                                                      gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGa                                             <  2:308060/148‑1 (MQ=255)                                                                                      gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGa                                             <  6:26833/148‑1 (MQ=255)                                                                                      gAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAACATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGa                                             <  5:325613/148‑1 (MQ=255)                                                                                             gcCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGc                                     <  2:118254/149‑1 (MQ=255)                                                                                                aGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATg                                  >  6:321900/1‑149 (MQ=255)                                                                                                  aTGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGTGCGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACTCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGt                                 >  1:75093/1‑149 (MQ=255)                                                                                                   tGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTc                               <  4:205852/149‑1 (MQ=255)                                                                                                     ccGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCa                              <  3:274057/148‑1 (MQ=255)                                                                                                      cGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGc                            <  3:150090/149‑1 (MQ=255)                                                                                                       gACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCa                           <  5:34422/149‑1 (MQ=255)                                                                                                                gACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCgag                  <  1:220444/149‑1 (MQ=255)                                                                                                                       aaGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCTATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCg            <  7:45555/148‑1 (MQ=255)                                                                                                                       aaGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCgg           <  4:15814/149‑1 (MQ=255)                                                                                                                        aGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCgg           >  7:277363/1‑148 (MQ=255)                                                                                                                           ctgatatcTTCGACTCGCAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGa         >  5:91599/8‑147 (MQ=255)                                                                                                                                atcTTCGACTCGCAGC‑CGAAGCTATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCATCAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTg  >  3:193484/3‑149 (MQ=255)                                                                                                                                       aCTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCgg           <  2:257502/133‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          ACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTG  >  NC_002947/1472575‑1472850

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