Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,472,712	G→C	10.2%	Q39E (CAG→GAG)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,472,712	0	G	C	10.2%	136.2 / 7.8	49	Q39E (CAG→GAG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base G (17/27);  new base C (3/2);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-02

GGTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTG  >  NC_002947/1472566‑1472850                                                                                                                                                   |                                                                                                                                            ggTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGc                                                                                                                                            <  8:222240/148‑1 (MQ=255)  gTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGc                                                                                                                                            >  2:220444/1‑147 (MQ=255)          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                                                                                               >  4:274057/1‑149 (MQ=255)                                            tACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAg                                                                                               >  3:40088/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)                                                                                                                                                   |                                                                                                                                            GGTGGCATTACCCATTTTCTGCACCTTGCGGCGCAGGTACGGGTACACGCCGTAGAGGAAGATCATGCCGCGGATGAAATGGGTGGCGCCCATGGAGTAGCGCCAGATGCCGACTGCGCCGACCAGGAAGATGAAGTGCTTCGACTGCGAGT‑CGAAGATATCAGCTGGCAGGGCCAGGGCGATCAGCATGAGCAGGCTCATGTAGAGCAGCCACCCGGCGCACTGCAACAGCACTGTCTGGAGCCTTTGCATGTTCAGCATCCGTCGAGAATTCGGGAAAGGGTG  >  NC_002947/1472566‑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: ‑