Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	243,508	T→A	8.8%	*221C (TGA→TGT) ‡	fkl ←	FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	243,508	0	T	A	8.8%	146.3 / 3.6	54	*221C (TGA→TGT) ‡	fkl	FKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base A (3/2);  total (30/24)
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.22e-01

CAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGG  >  NC_002947/243367‑243650                                                                                                                                              |                                                                                                                                                cAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGATGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGc                                                                                                                                          >  5:298742/1‑149 (MQ=255)     gACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAAc                                                                                                                                      >  4:120570/1‑149 (MQ=255)        acCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGc                                                                                                                                   <  4:215931/149‑1 (MQ=255)          agAACGCTTTGACGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATCAGCGCTTCGACCATGAAGACTGCCTGATAGAAGTGCTGCACAACGCACACAAGGAAGAGCGTGCAGTGCAGCACTGAGGTTGGCGTCAGTCGCCAACCGAca                                                                                                                                 >  4:3937/2‑149 (MQ=255)               gCTTTGCCGAGAAGCTGGCCATATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGt                                                                                                                            <  5:185942/149‑1 (MQ=255)                                gCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATAcc                                                                                                             <  3:273430/147‑1 (MQ=255)                                gCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATAcc                                                                                                             >  3:151171/1‑147 (MQ=255)                                gCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATAcc                                                                                                             >  3:146175/1‑147 (MQ=255)                                        gggCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                 >  2:197664/1‑135 (MQ=255)                                        gggCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGaa                                                                                                                 <  1:197664/135‑1 (MQ=255)                                         ggCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGgtgt                                                                                                   <  3:19752/148‑1 (MQ=255)                                           cGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAgg                                                                                                <  4:146175/149‑1 (MQ=255)                                           cGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAgg                                                                                                <  4:151171/149‑1 (MQ=255)                                                   tgCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAt                                                                                         >  7:29984/1‑148 (MQ=255)                                                    gCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAt                                                                                         >  2:315330/1‑147 (MQ=255)                                                     cATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAg                                                                                     <  1:135252/149‑1 (MQ=255)                                                      aTGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAgg                                                                                    >  5:42218/1‑149 (MQ=255)                                                            gcTTCGATCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACTAGGAAGATGGTGCAGTGGAGGCCTGAGGTTTGCGTCAGTCGCCAACCGCCAGGAGTTCGATCTCGAATACCAAGGGTGTGTAGGGTGCAAT‑CAGGTCTccc                                                                              <  3:3937/149‑1 (MQ=255)                                                              ttCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGc                                                                            >  5:209459/1‑149 (MQ=255)                                                                cGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGc                                                                            <  6:291010/147‑1 (MQ=255)                                                                   gCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCt                                                                       >  4:255702/1‑149 (MQ=255)                                                                            cTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCAt                                                              <  8:29984/149‑1 (MQ=255)                                                                                      gAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAATCCAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTgcg                                                     <  1:45962/149‑1 (MQ=255)                                                                                          tgctgcACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGaa                                                <  6:42218/149‑1 (MQ=255)                                                                                              gcACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAg                                                                                     >  4:169898/1‑108 (MQ=255)                                                                                              gcACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGaa                                            >  8:170987/1‑149 (MQ=255)                                                                                                acaACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAAt                                           >  1:157448/1‑148 (MQ=255)                                                                                                  aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAg                                               <  2:91899/142‑1 (MQ=255)                                                                                                  aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAg                                               >  1:91899/1‑142 (MQ=255)                                                                                                  aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTAc                                        >  6:237600/1‑149 (MQ=255)                                                                                                  aaCGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTAc                                        >  6:135478/1‑149 (MQ=255)                                                                                                     gCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGTGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAg                                     <  3:120570/149‑1 (MQ=255)                                                                                                        cacaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGgcg                                  >  8:39765/1‑149 (MQ=255)                                                                                                         acaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGgcgc                                 >  3:52517/1‑149 (MQ=255)                                                                                                         acaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGgcgc                                 >  4:221723/1‑149 (MQ=255)                                                                                                         acaAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGgcgc                                 >  2:238866/1‑149 (MQ=255)                                                                                                          caAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGgcgc                                 >  2:224238/1‑148 (MQ=255)                                                                                                           aaGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGAGCCa                               >  1:160362/1‑149 (MQ=255)                                                                                                               aaGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTg                           >  2:280511/1‑149 (MQ=255)                                                                                                               aaGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTg                           <  1:315330/149‑1 (MQ=255)                                                                                                                    ggTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGcc                      >  8:178874/1‑149 (MQ=255)                                                                                                                        cagtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgt                   <  1:238866/148‑1 (MQ=255)                                                                                                                         agtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGTGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgtgt                 <  6:298742/149‑1 (MQ=255)                                                                                                                         agtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGgtgt                 <  2:157448/149‑1 (MQ=255)                                                                                                                          gtgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGtgtg                >  8:51808/1‑149 (MQ=255)                                                                                                                           tgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                           >  3:317866/1‑105 (MQ=255)                                                                                                                           tgcagGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAgg                                                           <  4:317866/105‑1 (MQ=255)                                                                                                                                   tGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCgg       <  1:224238/149‑1 (MQ=255)                                                                                                                                    gcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTgg                          >  6:64082/3‑129 (MQ=255)                                                                                                                                    gcggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTgg                          <  5:64082/127‑1 (MQ=255)                                                                                                                                      ggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGc      >  2:278510/1‑147 (MQ=255)                                                                                                                                      ggTTGGCGACTGACGCCAACCTCCAGCAGTTCGATCTCGAATACCAAGGGGGTCTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCaa    <  1:278510/149‑1 (MQ=255)                                                                                                                                        ttGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAgg  <  5:237600/149‑1 (MQ=255)                                                                                                                                        ttGGCGACTGACGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGACATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATATGGCGCCGGTGTGCATTTTCGGCaacg  >  7:24427/1‑147 (MQ=255)                                                                                                                                              |                                                                                                                                                CAAGGACACCGAACGCTTTGCCGAGAAGCTGGCCAAATTGGGCGGGTTGTTGCATGAGCGCTTCGAGCTTGAAGACTGCCTGATCGAAGTGCTGCACAACGCACACAAGGAAGAGGGTGCAGTGCAGGCCTGAGGTTGGCGTCAGTCGCCAACCGCCAGCAGTTCGATCTCGAATACCAAGGGGGTGTAGGGTGCAAT‑CAGGTCTCCCGCCCCCTCGGCGCCATAGGCTTGCGCTGAAGGAATTACCAGGCGCCATTTGGCGCCGGTGTGCATTTTCGGCAAGG  >  NC_002947/243367‑243650

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