Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,065,068	C→T	8.3%	A259A (GCG→GCA)	PP_0923 ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,065,068	0	C	T	8.3%	139.9 / 2.4	48	A259A (GCG→GCA)	PP_0923	acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/21);  new base T (2/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GCGGATGTCTGACAGTTAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACAC  >  NC_002947/1064925‑1065202                                                                                                                                                |                                                                                                                                       gCGGATGTCTGACAGTTAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAG‑CCGCGCCTGc                                                                                                                                  <  2:12334/149‑1 (MQ=255)          tGACAGTTAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTgg                                                                                                                          >  7:263626/1‑149 (MQ=255)              aGTTAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGgcgc                                                                                                                       >  4:260135/1‑148 (MQ=255)                ttAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGTGCCTGCCGACGTGGCGCCAt                                                                                                                    >  1:89264/1‑149 (MQ=255)                ttAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGTGCCTGCCGACGTGGCGCCAt                                                                                                                    >  1:94107/1‑149 (MQ=255)                   gagGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCAtctc                                                                                                                 >  4:81717/1‑149 (MQ=255)                              gTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCaa                                                                                                      <  3:260135/149‑1 (MQ=255)                                 aTACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGTGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAAcc                                                                                                    <  2:89264/148‑1 (MQ=255)                                 aTACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGTGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAAcc                                                                                                    <  2:94107/148‑1 (MQ=255)                                     tGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGaa                                                                                               >  1:5386/1‑149 (MQ=255)                                     tGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGaa                                                                                               >  8:152155/1‑149 (MQ=255)                                     tGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGaa                                                                                               >  4:210983/1‑149 (MQ=255)                                       gCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACa                                                                                             >  4:151876/1‑149 (MQ=255)                                                cgGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGt                                                                                    <  2:53674/149‑1 (MQ=255)                                                       aCCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGc                                                                             <  3:151876/149‑1 (MQ=255)                                                         ccGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGAGAAAGGGCCAGAACAGGCAGTGAAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCt                                                                           >  2:319184/1‑149 (MQ=255)                                                          cGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTg                                                                          <  3:81717/149‑1 (MQ=255)                                                                tACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCAt                                                                                                                    >  5:370595/1‑101 (MQ=255)                                                                tACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCAt                                                                                                                    <  6:370595/101‑1 (MQ=255)                                                                     gCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAAc                                                                                                     >  5:151277/1‑111 (MQ=255)                                                                     gCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAAc                                                                                                     <  6:151277/111‑1 (MQ=255)                                                                     gCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGc                                                                >  5:130189/1‑148 (MQ=255)                                                                       aaTGCTCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAAc                                                             <  2:348294/149‑1 (MQ=255)                                                                                 gtgtAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCg                                                   >  5:69671/1‑149 (MQ=255)                                                                                      ggAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGt                                              >  1:293006/1‑149 (MQ=255)                                                                                           gggTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCtt                                         <  2:5386/149‑1 (MQ=255)                                                                                              ttCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGcc                                      <  7:152155/149‑1 (MQ=255)                                                                                                 aCCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCAc                                   >  3:332023/1‑149 (MQ=255)                                                                                                 aCCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCAc                                   >  7:60240/1‑149 (MQ=255)                                                                                                 aCCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCAc                                   <  3:210983/149‑1 (MQ=255)                                                                                                   ccGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGc                                 >  4:305559/1‑149 (MQ=255)                                                                                                   ccGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGc                                 <  2:355799/149‑1 (MQ=255)                                                                                                       gAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGc                                                                             <  4:114668/101‑1 (MQ=255)                                                                                                       gAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGc                                                                             >  3:114668/1‑101 (MQ=255)                                                                                                           gggCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTAtt                         <  8:60240/149‑1 (MQ=255)                                                                                                           gggCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTAtt                         <  8:263626/149‑1 (MQ=255)                                                                                                             gCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGt                                                                         >  6:275681/1‑99 (MQ=255)                                                                                                             gCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGt                                                                         <  5:275681/99‑1 (MQ=255)                                                                                                                     aGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTg                                                                 >  7:219692/1‑99 (MQ=255)                                                                                                                     aGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTg                                                                 <  8:225567/99‑1 (MQ=255)                                                                                                                     aGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTg                                                                 <  8:219692/99‑1 (MQ=255)                                                                                                                     aGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTg                                                                 >  7:225567/1‑99 (MQ=255)                                                                                                                             cAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGAt        <  6:69671/148‑1 (MQ=255)                                                                                                                              aGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATg       >  3:175703/1‑148 (MQ=255)                                                                                                                              aGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATg       >  5:291634/1‑148 (MQ=255)                                                                                                                               ggCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGca     <  6:130189/149‑1 (MQ=255)                                                                                                                                       ggCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCacac  >  6:319281/1‑144 (MQ=255)                                                                                                                                       ggCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACCTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCacac  <  5:319281/144‑1 (MQ=255)                                                                                                                                                |                                                                                                                                       GCGGATGTCTGACAGTTAGAGGTGTGTGGGTGATACTGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGCAATGCGCTCGGTGTAGGAGCGGGTTCACCCGCGAAAGGGCCAGAACAGGCAGTGCAGGCATTCAGGCAGCCCGCGCCTGCCGACGTGGCGCCATCTCCACGTCTTCAACCCCGAACAGCGCCAGGTGAATCGCCTGTTGCGCCTGCAACGCCAGCACCGCACGTTCCTTGCCCACGCTACCTATTGGCTGCAGCAGGTGGATGCACAC  >  NC_002947/1064925‑1065202

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