Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,317,213	C→G	7.4%	P608P (CCG→CCC)	treSA ←	trehalose synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,317,213	0	C	G	7.4%	163.2 / 4.4	54	P608P (CCG→CCC)	treSA	trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base C (24/26);  new base G (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCC  >  NC_002947/3317068‑3317352                                                                                                                                                  |                                                                                                                                             ccGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCggg                                                                                                                                           <  7:182438/149‑1 (MQ=255) ccGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGGGATCTGGATAGCGCGCCCGGCTGGCAGCTCATGCTCCATCACCAGCAGCCCCggg                                                                                                                                           <  4:162226/149‑1 (MQ=255)  cGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCACGTCTGGCAGATCATGCACCATCACCAGCAGCCCCGGgc                                                                                                                                          >  6:289465/1‑149 (MQ=255)       cGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgct                                                                                                                                      >  3:128851/1‑148 (MQ=255)        gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                    <  6:79297/149‑3 (MQ=255)        gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgct                                                                                                                                      <  5:289465/147‑1 (MQ=255)        gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgc                                                                                                                                       <  6:75456/146‑1 (MQ=255)             tCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCt                                                                                                                                 >  5:115963/1‑147 (MQ=255)               aTGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                              <  7:284880/148‑1 (MQ=255)               aTGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                              <  6:115963/148‑1 (MQ=255)                  tCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgc                                                                                                                                       >  3:183760/1‑136 (MQ=255)                  tCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgctgc                                                                                                                                       <  4:183760/136‑1 (MQ=255)                      ccaccGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAg                                                                                                                       >  5:194763/1‑148 (MQ=255)                        accGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAc                                                                                                                    <  1:56651/149‑1 (MQ=255)                        accGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCAc                                                                                                                    >  2:56651/1‑149 (MQ=255)                            gCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCa                                                                                                                                    <  1:321884/129‑1 (MQ=255)                            gCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCa                                                                                                                                    >  2:321884/1‑129 (MQ=255)                             cccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCGCAGGTACCAGCCCCTgcc                                                                                                               >  8:246668/1‑147 (MQ=255)                              ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                    >  4:130446/1‑125 (MQ=255)                              ccAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                    <  3:130446/127‑3 (MQ=255)                                   tgtgAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCtgg                                                                                                         <  4:128851/149‑1 (MQ=255)                                       aaCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGc                                                                                                     >  3:31094/1‑149 (MQ=255)                                        aCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGc                                                                                                     >  5:254454/1‑148 (MQ=255)                                        aCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑g                                                                                                   <  2:154365/149‑1 (MQ=255)                                         cccGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑Ga                                                                                                  >  7:150111/1‑149 (MQ=255)                                          ccGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGGGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GAt                                                                                                 >  6:63566/1‑149 (MQ=255)                                                        ttCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGgc                                                                                   <  7:232830/149‑1 (MQ=255)                                                                     gcgTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCcagca                                                                       <  4:31094/148‑1 (MQ=255)                                                                      cgTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGt                                                                     <  6:254454/149‑1 (MQ=255)                                                                            ggCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                >  7:34445/1‑148 (MQ=255)                                                                             gCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTttcttc                                                               <  7:38692/148‑1 (MQ=255)                                                                             gCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                <  8:150111/147‑1 (MQ=255)                                                                               cGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCAcc                                                            >  5:21544/1‑149 (MQ=255)                                                                                  aGTTCAGCGCAGGGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTg                                                         <  5:165206/149‑1 (MQ=255)                                                                                                   tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                              >  4:307081/1‑64 (MQ=255)                                                                                                   tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                              <  3:307081/64‑1 (MQ=255)                                                                                                   tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCa                                                       >  8:193586/1‑134 (MQ=255)                                                                                                   tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCa                                                       <  7:193586/134‑1 (MQ=255)                                                                                                     gATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCActct                                      >  6:9553/1‑149 (MQ=255)                                                                                                     gATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCActct                                      >  5:17627/1‑149 (MQ=255)                                                                                                       tACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCg                                    <  8:34445/149‑1 (MQ=255)                                                                                                         ccGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCTGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTc                                  <  5:9553/149‑1 (MQ=255)                                                                                                         ccGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGt                                   <  6:17627/148‑1 (MQ=255)                                                                                                          cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                <  1:60176/134‑1 (MQ=255)                                                                                                          cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                >  2:60176/1‑134 (MQ=255)                                                                                                              cccGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCt                             <  6:194763/149‑1 (MQ=255)                                                                                                                cGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGc                           <  6:21544/149‑1 (MQ=255)                                                                                                                             aTGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                            >  6:342504/1‑119 (MQ=255)                                                                                                                             aTGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAAt                                            <  5:342504/119‑1 (MQ=255)                                                                                                                             aTGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGc              >  5:300898/1‑149 (MQ=255)                                                                                                                                 aCCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGGTTGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCg            >  1:271936/1‑148 (MQ=255)                                                                                                                                      cACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAgcg     >  2:286246/1‑149 (MQ=255)                                                                                                                                        tgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGcc   >  4:115027/3‑149 (MQ=255)                                                                                                                                         cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGccc  >  3:221492/1‑149 (MQ=255)                                                                                                                                                 ccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCTGGCAAAt                                                <  7:219929/95‑1 (MQ=255)                                                                                                                                                 ccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAAt                                                >  8:219929/1‑95 (MQ=255)                                                                                                                                                  |                                                                                                                                             CCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGC‑GATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCC  >  NC_002947/3317068‑3317352

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