Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,634,523	Δ1 bp	7.1%	coding (985/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,523	0	G	.	7.1%	178.6 / 6.7	56	coding (985/1872 nt)	PP_2304	peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base . (2/2);  total (26/30)
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

GACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACAAACT  >  NC_002947/2634376‑2634668                                                                                                                                                     |                                                                                                                                                  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                                                             >  4:280863/1‑148 (MQ=255)                                                                                       tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                             >  2:229117/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                                                                                         >  8:174092/1‑106 (MQ=255)                                                                                                     ccGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGt                                                                                         <  7:174092/106‑1 (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)                                                                                                                                                tataCGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACaaa    <  3:280863/149‑1 (MQ=255)                                                                                                                                                tataCGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACaa     <  5:140946/148‑1 (MQ=255)                                                                                                                                                  taCGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACAAACt  >  4:98682/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGT‑GTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACAAACT  >  NC_002947/2634376‑2634668

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