Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,471,405	A→C	13.8%	S477A (TCG→GCG)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,471,405	0	A	C	13.8%	53.4 / 5.9	29	S477A (TCG→GCG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base A (12/13);  new base C (2/2);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.55e-01

GGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCA  >  NC_002947/1471261‑1471548                                                                                                                                                 |                                                                                                                                                ggCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGAcc                                                                                                                                            >  2:8194/1‑150 (MQ=255)   cTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                          <  2:265183/150‑1 (MQ=255)   cTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTGCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                          <  1:167297/150‑1 (MQ=255)         aCTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACCGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGCACCAGGTGTTg                                                                                                                                    <  2:6660/150‑1 (MQ=255)             aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgacc                                                                                                                                      <  2:14197/144‑6 (MQ=255)             aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatgacc                                                                                                                                      >  1:14197/1‑139 (MQ=255)                   gACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttgtt                                                                                                                          >  1:296957/1‑150 (MQ=255)                                                    ttGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGACAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGTAAGCTGGCGAGGTCACGCTTGAGGGCAGTgg                                                                                         <  1:157127/150‑1 (MQ=255)                                                     tGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCACCGTGGCGCCTGGTTCGCGAGGCCCCGGTGTTGAACCACTGTTGAAGGCTGGCGAGGTCACGCTTGAGGGCAGTggg                                                                                        <  1:8194/150‑1 (MQ=255)                                                           tGTTAGATCCGATCCGTGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGACAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                  >  2:100464/1‑150 (MQ=255)                                                           tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                  >  2:227315/1‑150 (MQ=255)                                                             ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                          >  2:295854/1‑92 (MQ=255)                                                             ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                          <  1:295854/92‑1 (MQ=255)                                                                     gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttg                                                                                                                            >  1:372299/1‑98 (MQ=255)                                                                     gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttg                                                                                                                            <  2:372299/98‑1 (MQ=255)                                                                          gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTACACGACTGCtt                                                                   >  2:278143/1‑150 (MQ=255)                                                                                         gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTc                                                                                                          <  1:420649/96‑1 (MQ=255)                                                                                         gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTc                                                                                                          >  2:420649/1‑96 (MQ=255)                                                                                         gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaagaa                                                    <  1:306143/150‑1 (MQ=255)                                                                                              tGAACAGCACAGCAACGAAGATGCTGACAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGt                                               <  1:100464/150‑1 (MQ=255)                                                                                                        agcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcat                                     >  2:20643/1‑150 (MQ=255)                                                                                                                   tGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCga                          <  2:296957/150‑1 (MQ=255)                                                                                                                            ccGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGt                 >  2:213303/1‑150 (MQ=255)                                                                                                                               agaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgcg                              >  2:152041/1‑134 (MQ=255)                                                                                                                               agaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgcg                              <  1:152041/134‑1 (MQ=255)                                                                                                                                     tCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga        <  2:326378/150‑1 (MQ=255)                                                                                                                                     acctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTCGTCGAGGCGGAAGAACACGTAGATCTtc                                       >  2:246262/4‑119 (MQ=255)                                                                                                                                     acctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTCGTCGAGGCGGAAGAACACGTAGATCTtc                                       <  1:246262/116‑1 (MQ=255)                                                                                                                                           tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCa  >  1:216406/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                GGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCA  >  NC_002947/1471261‑1471548

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑