Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,344,302	A→T	7.3%	intergenic (+33/‑123)	PP_3817 → / → PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,344,302	0	A	T	7.3%	142.3 / 2.3	55	intergenic (+33/‑123)	PP_3817/PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (26/24);  new base T (2/2);  total (28/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 = 4.91e-01

GGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGC  >  NC_002947/4344155‑4344449                                                                                                                                                     |                                                                                                                                                    gggCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAg                                                                                                                                                    <  4:129897/149‑1 (MQ=255)      tGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCAccg                                                                                                                                               >  8:135426/1‑149 (MQ=255)        ccGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCAccgc                                                                                                                                              <  8:257282/148‑1 (MQ=255)        ccGATGCAGCCGGGCGACAGCGTGAAAGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCAccgcc                                                                                                                                             >  1:10923/1‑149 (MQ=255)         cGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATAGAGCAccgccc                                                                                                                                            >  5:53957/1‑149 (MQ=255)               aGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCAccgc                                                                                                                                              <  6:247538/141‑1 (MQ=255)               aGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCAccgc                                                                                                                                              >  5:247538/1‑141 (MQ=255)                       gACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGAt                                                                                                                              <  2:10923/149‑1 (MQ=255)                         cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                            <  5:264926/149‑1 (MQ=255)                         cAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCt                                                                                                                            >  2:270700/1‑149 (MQ=255)                             gTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCa                                                                                                                        <  7:96899/149‑1 (MQ=37)                                aaGGTGCGGTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCTCATTAACGCTCTTATACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                      <  3:240340/148‑1 (MQ=255)                                      cGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCa                                                                                                               >  5:267808/1‑149 (MQ=255)                                            gCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTg                                                                                                         >  1:300877/1‑149 (MQ=255)                                             cAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTgg                                                                                                        <  1:96728/149‑1 (MQ=255)                                                  gCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGc                                                                                                    <  6:267808/148‑1 (MQ=255)                                                          ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                           >  1:249086/1‑149 (MQ=255)                                                          ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                           >  2:44403/1‑149 (MQ=255)                                                          ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGAGCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                           >  4:126473/1‑149 (MQ=255)                                                           tgaatatatcCCGAATGGGCCGAGCGCGACCTGACCAAGCCGGCGGGGGCGAATTAACGCTCTTGAAGGTTGCAGCGCGG‑GCCTTATCGCGCACCGCCCGCGCCCCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                           <  2:300877/142‑1 (MQ=255)                                                            gCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                      >  3:240771/1‑120 (MQ=255)                                                            gCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                      <  4:240771/120‑1 (MQ=255)                                                                      cGAATGGCCCGCCAGCCACCTGAGCACGGCGGCCGGGGCGCATTAACGCTCTTGCACGTTGCCCCGCCT‑GCCACATCGAGCACCGCCCGCGCGCCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatca                                                                               <  7:135426/149‑1 (MQ=255)                                                                        aaTGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTatcatc                                                                             <  1:44403/149‑1 (MQ=255)                                                                            ggCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGc                                                                         <  1:270700/149‑1 (MQ=255)                                                                              ccGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCa                                                                       <  2:249086/149‑1 (MQ=255)                                                                                        ccTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                              >  6:182468/1‑148 (MQ=255)                                                                                        ccTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                              <  5:182468/148‑1 (MQ=255)                                                                                          tGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCTGGCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATc                                                            >  1:264223/1‑149 (MQ=255)                                                                                           gAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCCCCGCCCGCGCGGCACACGAGATGGCAGGCGCTGCACGATTTGCGGCGTAAGCCCCGGCCCCTTATCATCACGCCACTGTCGGActctt                                                          >  6:568/1‑148 (MQ=255)                                                                                            cgCAAGGCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTg                                                          <  5:568/147‑1 (MQ=255)                                                                                             gccaGTCGGCCGGGGCGCATTCACGCTCTTGTCCGTGGCAGCGCCT‑GCCCGAGCGAGCCCCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGcc                                                        <  6:178866/146‑1 (MQ=255)                                                                                               aaGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                  >  7:148067/1‑137 (MQ=255)                                                                                               aaGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                  <  8:148067/137‑1 (MQ=255)                                                                                                 gTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGc                                                     >  7:73160/1‑148 (MQ=255)                                                                                                    ggCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGTTCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTc                                                                  <  5:312924/132‑1 (MQ=255)                                                                                                               aTTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                  >  6:289927/1‑137 (MQ=255)                                                                                                               aTTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCaa                                                  <  5:289927/137‑1 (MQ=255)                                                                                                                  aaCGCTCTTGTACGTTGCAGCGCCT‑GCCACATAGAGCACCGCCCCCGCGGCACACCATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATAACCCCACTGTCGTAATATGCCGGCCAACAGTATGCcacaag                                    >  6:201944/1‑146 (MQ=255)                                                                                                                      ctctTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCAcc                               >  6:326415/1‑149 (MQ=255)                                                                                                                         ttGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCCACTCTATGCCACAGGGCACCCgg                            >  8:320409/1‑149 (MQ=255)                                                                                                                                       gcCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc              <  5:326415/149‑1 (MQ=255)                                                                                                                                       gcCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc              >  2:269896/1‑149 (MQ=255)                                                                                                                                            gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTg               >  1:318993/1‑144 (MQ=37)                                                                                                                                            gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTg               <  2:318993/144‑1 (MQ=37)                                                                                                                                             ccAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc         >  2:173769/1‑149 (MQ=37)                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCTGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc         >  4:118115/1‑149 (MQ=255)                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc              <  4:79974/144‑1 (MQ=255)                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc              >  3:79974/1‑144 (MQ=255)                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt          >  4:298961/1‑148 (MQ=255)                                                                                                                                             ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGt          <  3:298961/148‑1 (MQ=255)                                                                                                                                              cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                        <  1:308667/117‑1 (MQ=255)                                                                                                                                              cAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCca                                        >  2:308667/1‑117 (MQ=255)                                                                                                                                                    gAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCTGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc  <  3:118115/149‑1 (MQ=255)                                                                                                                                                    gAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTtgctgc  >  3:31783/1‑149 (MQ=255)                                                                                                                                                     aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                           <  2:296685/107‑1 (MQ=255)                                                                                                                                                     aGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATg                                           >  1:296685/1‑107 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGGCATGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGC  >  NC_002947/4344155‑4344449

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