Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,317,256	C→G	13.0%	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,256	0	C	G	13.0%	83.8 / 2.9	29	P608P (CCG→CCC)	treSA	trehalose synthase A
Reads supporting (aligned to +/- strand):  ref base C (14/11);  new base G (2/2);  total (16/13)
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.00e+00

CCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGC  >  NC_002947/3317111‑3317403                                                                                                                                                  |                                                                                                                                                    ccGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGgc                                                                                                                                                 <  1:101706/150‑1 (MQ=255)       cGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGgctgctgg                                                                                                                                           <  1:118384/150‑3 (MQ=255)        gTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCAc                                                                                                                                          >  2:14439/1‑150 (MQ=255)              cATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGAAAACTAACCTAGACACAACACCCTACGGGCAGTTCAAGCGCGATCACCAGCAGCCCCGGGCTGCTCACCTCAgg                                                                                                                                    <  1:14439/150‑1 (MQ=255)                tGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTa                                                                                                                                  >  1:167820/1‑150 (MQ=255)                            gCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGAc                                                                                                                      <  1:520026/150‑1 (MQ=255)                                        aCCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGa                                                                                                          <  2:85169/150‑1 (MQ=255)                                                            tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                      <  1:471935/150‑1 (MQ=255)                                                            tCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGcaca                                                                                      >  1:530540/1‑150 (MQ=255)                                                                 aaTCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                     <  1:345111/98‑1 (MQ=255)                                                                 aaTCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAg                                                                                                                                     >  2:345111/1‑98 (MQ=255)                                                                   tCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCcagca                                                                               >  2:23496/1‑150 (MQ=255)                                                                        tCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttc                                                                          <  2:615687/150‑1 (MQ=255)                                                                          cTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTttctt                                                                        >  1:287982/1‑150 (MQ=255)                                                                                 aaGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTg                                                                 <  1:23496/150‑1 (MQ=255)                                                                                       aGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCa                                                           >  1:61716/1‑150 (MQ=255)                                                                                           cAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCATCACCTGACGGCt                                                                                                                  >  2:786675/1‑91 (MQ=255)                                                                                           cAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCt                                                                                                                  <  1:786675/91‑1 (MQ=255)                                                                                            aGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGa                                                      >  1:339249/1‑150 (MQ=255)                                                                                                aTCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCa                                                  <  2:530540/150‑1 (MQ=255)                                                                                                   tGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCActc                                               >  1:584105/1‑150 (MQ=255)                                                                                                          cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagcca                                                                                  >  2:374303/1‑108 (MQ=255)                                                                                                          cgcgCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCAcagcca                                                                                  <  1:374303/108‑1 (MQ=255)                                                                                                                                   cATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGtataa               >  2:522731/1‑149 (MQ=255)                                                                                                                                     tCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAgcg             >  1:568925/1‑150 (MQ=255)                                                                                                                                        tgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGc                      <  1:698070/136‑1 (MQ=255)                                                                                                                                        tgagcagcCCGGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGc                      >  2:698070/3‑138 (MQ=255)                                                                                                                                        ccAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCCAACCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGccc          >  1:696834/1‑150 (MQ=255)                                                                                                                                         cagcagCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCg         >  1:768213/1‑150 (MQ=255)                                                                                                                                                cccGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTATCGTCCAGCTGCCTGTCTAGGCTGCCGTATATGGCGCGGGCAggc  >  1:311426/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CCGTCTCGTTGATCATGTCGACCACCGGCCCAGGTGTGAACCCGGTCAACAGCAGTTCCTCGGCAATCGCGTCCTGGCCGAAGTTCAGCGCAGTGATCTGGATACCGCGCCCGGCTGGCAGCTCATGCACCATCACCAGCAGCCCCGGGCTGCTCACCTCAGGTACCAGCACCTGACGGCTGGTGGCGATGCCGTAGGCCTGGCGCACAGCCAGCAGTTTCTTCACCTTGCAGGCAAATGAATCACTCTCGTCCAGCTGCCTGTCCAGGCTGCCGTATAGCGCCCGGGCACGC  >  NC_002947/3317111‑3317403

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