Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,344,305	A→G	7.2%	intergenic (+36/‑120)	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,305	0	A	G	7.2%	147.7 / 2.7	57	intergenic (+36/‑120)	PP_3817/PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (27/24);  new base G (2/2);  total (30/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 = 6.45e-02

TGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCT  >  NC_002947/4344160‑4344450                                                                                                                                                   |                                                                                                                                                  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                                                                                            >  2:44403/1‑149 (MQ=255)                                                     ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                            >  1:249086/1‑149 (MQ=255)                                                     ttGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGAGCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                            >  4:126473/1‑149 (MQ=255)                                                      tgaatatatcCCGAATGGGCCGAGCGCGACCTGACCAAGCCGGCGGGGGCGAATTAACGCTCTTGAAGGTTGCAGCGCGG‑GCCTTATCGCGCACCGCCCGCGCCCCACACTATCTGGCAGGCGCTGCACTATTGGAGGCGTACGcccc                                                                                            <  2:300877/142‑1 (MQ=255)                                                       gCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                       <  4:240771/120‑1 (MQ=255)                                                       gCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAgg                                                                                                                       >  3:240771/1‑120 (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                                                               <  5:182468/148‑1 (MQ=255)                                                                                   ccTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTaa                                                               >  6:182468/1‑148 (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                <  2:318993/144‑1 (MQ=37)                                                                                                                                       gCCAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTg                >  1:318993/1‑144 (MQ=37)                                                                                                                                        ccAGATCGTGTGCCGCGCGGGCGGTGCTCGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc          >  2:173769/1‑149 (MQ=37)                                                                                                                                        ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCTGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTc          >  4:118115/1‑149 (MQ=255)                                                                                                                                        ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc               >  3:79974/1‑144 (MQ=255)                                                                                                                                        ccAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc               <  4:79974/144‑1 (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)                                                                                                                                                 gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct  >  2:98609/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TGCCGATGCAGCCGGGCGACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCT  >  NC_002947/4344160‑4344450

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