Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,614,521	G→T	7.4%	A132S (GCG→TCG)	pet →	protein Pet18

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,614,521	0	G	T	7.4%	170.3 / 3.5	54	A132S (GCG→TCG)	pet	protein Pet18
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base T (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 = 9.74e-01

ACACTTACTTCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCA  >  NC_002947/3614381‑3614660                                                                                                                                             |                                                                                                                                            acacTTACTTCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGGCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTg                                                                                                                                     <  8:35076/149‑1 (MQ=255)   acttacttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGaa                                                                                                                                   <  6:257243/149‑1 (MQ=255)    cttacttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAAc                                                                                                                                  <  3:86051/149‑1 (MQ=255)     ttacttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACg                                                                                                                                 >  6:180930/1‑149 (MQ=255)        cttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTg                                                                                                                              >  3:207561/1‑149 (MQ=255)         ttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGtt                                                                                                                               <  8:185882/147‑1 (MQ=255)         ttCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGtt                                                                                                                               >  7:185882/1‑147 (MQ=255)           ccTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGc                                                                                                                             <  5:180930/147‑1 (MQ=255)            cTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGaa                                                                                                                                   >  1:75316/1‑140 (MQ=255)            cTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGaa                                                                                                                                   <  2:75316/140‑1 (MQ=255)                gTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGc                                                                                                                      <  2:127933/149‑1 (MQ=255)                 tGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGcc                                                                                                                     <  4:207561/149‑1 (MQ=255)                    cTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTg                                                                                                                   >  7:156018/1‑148 (MQ=255)                     ttCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTa                                                                                                                 >  3:218528/1‑149 (MQ=255)                     ttCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTa                                                                                                                 >  5:257289/1‑149 (MQ=255)                       cGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTAc                                                                                                                <  4:218528/148‑1 (MQ=255)                                  gtgtCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGgcg                                                                                                    >  5:3013/1‑149 (MQ=255)                                     tCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCt                                                                                                  >  6:16456/1‑148 (MQ=255)                                       cGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTc                                                                                                 >  5:263904/1‑147 (MQ=255)                                         cccGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGATTCAAGGCGATCATGCGCGACGCGGCAGCCACCCGAACCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCAGAACCTGGACTCGGCGCTCAAgg                                                                                             >  1:60523/1‑149 (MQ=255)                                          ccGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGgc                                                                                            <  6:263904/149‑1 (MQ=255)                                                  aCGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAAc                                                                                    <  6:34777/149‑1 (MQ=255)                                                     cAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAg                                                                                                                        <  1:228991/110‑1 (MQ=255)                                                     cAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAg                                                                                                                        >  2:228991/1‑110 (MQ=255)                                                        tGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGc                                                                              >  7:279759/1‑149 (MQ=255)                                                         gAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGgc                                                                                                     <  4:296564/125‑1 (MQ=255)                                                         gAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGgc                                                                                                     >  3:296564/1‑125 (MQ=255)                                                              tcatgATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGAc                                                                         <  3:300077/144‑1 (MQ=255)                                                                ccGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAAc                                                                      >  6:45444/1‑149 (MQ=255)                                                                   aTACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTc                                                                   <  4:300840/149‑1 (MQ=255)                                                                      ccGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTc                                                                >  2:130805/1‑149 (MQ=255)                                                                       cgcgCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTc                                                                >  2:281387/1‑148 (MQ=255)                                                                               ccGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAg                                                       <  5:16456/149‑1 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGTCGGTGCTg                                                                                                                                     >  8:256533/1‑70 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGTCGGTGCTg                                                                                                                                     <  1:300549/70‑1 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGTCGGTGCTg                                                                                                                                     <  7:256533/70‑1 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGTCGGTGCTg                                                                                                                                     >  2:300549/1‑70 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAgg                                                                                             <  6:59257/110‑1 (MQ=255)                                                                                cGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAgg                                                                                             >  5:59257/1‑110 (MQ=255)                                                                                          aaGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATg                                                            <  5:162332/133‑1 (MQ=255)                                                                                             gCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGAGAACTTCGTCCATGCCGAGTGGATTACCCTGCa                                         <  8:156018/149‑1 (MQ=255)                                                                                             gCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGt                                                      <  3:160003/136‑1 (MQ=255)                                                                                             gCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGt                                                      >  4:160003/1‑136 (MQ=255)                                                                                                   aTGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACaa                                   <  1:281387/149‑1 (MQ=255)                                                                                                        cgAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCgg                              <  1:130805/149‑1 (MQ=255)                                                                                                         gAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCgg                              >  3:30628/1‑148 (MQ=255)                                                                                                          aaGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCa                            >  4:131363/1‑149 (MQ=255)                                                                                                                cAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAg                      <  8:279759/149‑1 (MQ=255)                                                                                                                   ccACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCa                   <  5:45444/149‑1 (MQ=255)                                                                                                                           tCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGAtt                                                 >  3:298157/1‑111 (MQ=255)                                                                                                                           tCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGAtt                                                 <  4:298157/111‑1 (MQ=255)                                                                                                                             cTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGCCACCCCGGCATTCCAGGCATTCGTCGCtt         >  4:55198/1‑149 (MQ=255)                                                                                                                               aTGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTc       >  1:27788/1‑149 (MQ=255)                                                                                                                                    gCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCa  <  6:3013/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                            ACACTTACTTCCTGCGTGCCTTCGAGGCCTTGGGTGTCCGCCCGGCGCAACGCAGTGAACCGGCCGATACCGCGCCTACCGCAGGTTTCAAGGCGATCATGCGCGAAGCGGCAGCCACCCGTTCCTATGCAGCGGCACTGGCGGTGCTGAACGTTGCCGAAGGCCTGTACCTGGACTGGGCGCTCAAGGCGCCGAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCA  >  NC_002947/3614381‑3614660

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