Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,634,516:1	+C	7.7%	coding (978/1872 nt)	PP_2304 →	peptidyl‑prolyl cis‑trans isomerase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,634,516	1	.	C	7.7%	198.2 / 10.9	52	coding (978/1872 nt)	PP_2304	peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base . (23/25);  new base C (2/2);  total (25/27)
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.15e-01

TCACCGACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGAC  >  NC_002947/2634371‑2634663                                                                                                                                                   |                                                                                                                                                    tCACCGACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑Gta                                                                                                                                                  <  2:200306/149‑1 (MQ=255)      gACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGa                                                                                                                                             <  2:143520/149‑1 (MQ=255)          cccAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCg                                                                                                                                         <  3:25396/149‑1 (MQ=255)          cccAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGCGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCg                                                                                                                                         <  4:105645/149‑1 (MQ=255)              ggCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTCGTATAC‑ACCCgggcc                                                                                                                                     <  7:111736/149‑4 (MQ=255)                   aGGCCCGCGCCGAAGATATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGaa                                                                                                                                <  7:149470/149‑1 (MQ=255)                   aGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGaa                                                                                                                                <  2:72542/149‑1 (MQ=255)                   aGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGaa                                                                                                                                <  3:177673/149‑1 (MQ=255)                     gCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAaga                                                                                                                              <  6:293477/149‑1 (MQ=255)                       ccGCGCCGAAGAGATCGAGCAGCTTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAgg                                                                                                                            <  7:219313/149‑1 (MQ=255)                              gAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTa                                                                                                                     >  3:100602/1‑149 (MQ=255)                               aaGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGAcc                                                                                                                                           >  8:40177/1‑126 (MQ=255)                               aaGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGAcc                                                                                                                                           <  7:40177/126‑1 (MQ=255)                               aaGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTa                                                                                                                     >  8:434/1‑148 (MQ=255)                                               tCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACCTGCTGCAGGATGGGc                                                                                                    >  8:179871/1‑149 (MQ=255)                                                 tGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAg                                                                                                  >  2:329158/1‑149 (MQ=255)                                                   gCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAgg                                                                                                 >  1:166658/1‑148 (MQ=255)                                                           tGAGGACTTTGCCGCGTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCAc                                                                                        >  6:140946/1‑149 (MQ=255)                                                                              gCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTgg                                                                     <  2:294139/149‑1 (MQ=255)                                                                               ccAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGc                                                                    >  2:183226/1‑149 (MQ=255)                                                                                cAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCt                                                                   <  2:45175/149‑1 (MQ=255)                                                                                cAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGCATGGGCAGGTATCTGAACCGGTGCGCACCGAGTTTgcc                                                                    >  7:172165/1‑146 (MQ=255)                                                                                    gAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCa                                                               <  4:100602/149‑1 (MQ=255)                                                                                         ctctgATGACCCAGGCTCCGCCAACTCGGGTGGTGACCTCGGCTTTGCCGGCCCGTGT‑GTATAGGAGCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATc                                                             <  6:279712/146‑1 (MQ=255)                                                                                         ctctcAGGACCCAGGCTCCGCCAACACCGGTTGTGAACTCGGCTTTGCCGGCCCGGGT‑GTAAACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCCGTGCGCACCGAGTTTGGCTACCATc                                                             >  5:279712/1‑146 (MQ=255)                                                                                            tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                        >  2:229117/1‑148 (MQ=255)                                                                                            tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                        >  4:280863/1‑148 (MQ=255)                                                                                            tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                        >  1:160907/1‑148 (MQ=255)                                                                                              aGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                        <  1:183226/146‑1 (MQ=255)                                                                                              aGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCa                                                       <  7:434/147‑1 (MQ=255)                                                                                                gACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAgct                                                   >  3:30218/1‑149 (MQ=255)                                                                                                   ccAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAgctgct                                                >  6:32148/1‑149 (MQ=255)                                                                                                   ccAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAgctgct                                                >  6:16146/1‑149 (MQ=255)                                                                                                   ccAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAgctgct                                                >  8:94830/1‑149 (MQ=255)                                                                                                     aGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTgg                                              <  3:71224/149‑1 (MQ=255)                                                                                                       gCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGc                                            <  5:32148/149‑1 (MQ=255)                                                                                                          ccGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGt                                                                                    <  7:174092/106‑1 (MQ=255)                                                                                                          ccGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGt                                                                                    >  8:174092/1‑106 (MQ=255)                                                                                                            gCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCa                                       <  1:229117/149‑1 (MQ=255)                                                                                                               aaCACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGgc                                    <  5:16146/149‑1 (MQ=255)                                                                                                               aaCACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAACCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATc                                                             >  2:42388/1‑124 (MQ=255)                                                                                                                    cGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGgcgc                                  >  1:136373/1‑146 (MQ=255)                                                                                                                              cTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAg                     >  5:212044/1‑149 (MQ=255)                                                                                                                                 ggCTTTGCCGGCCCGGGT‑GTATTCGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCTGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGc                    <  6:316359/147‑1 (MQ=255)                                                                                                                                      tGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGcc              <  8:167653/148‑1 (MQ=255)                                                                                                                                        gaaggccGGGTCGTATAC‑ACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGc                                            <  4:33493/111‑1 (MQ=255)                                                                                                                                        gaaggccGGGTCGTATAC‑ACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGc                                            >  3:33493/6‑116 (MQ=255)                                                                                                                                        gaaggccGGGTCGTATAC‑ACCCGGCCTACGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGc           >  2:312452/6‑149 (MQ=255)                                                                                                                                        ccGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTg                                                                      <  5:124187/90‑1 (MQ=255)                                                                                                                                        ccGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTg                                                                      >  6:124187/1‑90 (MQ=255)                                                                                                                                        ccGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGc           >  5:159080/1‑149 (MQ=255)                                                                                                                                              cGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAg     <  4:30218/149‑1 (MQ=255)                                                                                                                                                 gt‑gtATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGAc  <  2:160907/149‑1 (MQ=255)                                                                                                                                                 gt‑gtATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGAc  <  1:329158/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TCACCGACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGAC  >  NC_002947/2634371‑2634663

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