Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,344,318	C→T	7.2%	intergenic (+49/‑107)	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,318	0	C	T	7.2%	151.8 / 2.0	55	intergenic (+49/‑107)	PP_3817/PP_3818	polyamine ABC transporter ATP‑binding protein/phosphate‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (25/26);  new base T (2/2);  total (27/28)
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.97e-01

GACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCT  >  NC_002947/4344177‑4344459                                                                                                                                               |                                                                                                                                              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                        >  2:269896/1‑149 (MQ=255)                                                                                                                 gcCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGc                        <  5:326415/149‑1 (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)                                                                                                                                gCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTgctgct           >  2:98609/1‑148 (MQ=255)                                                                                                                                   ccgcccgcGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc       <  8:73160/149‑1 (MQ=255)                                                                                                                                        cgcgcgGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCt  <  2:264223/149‑1 (MQ=255)                                                                                                                                              gCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc       <  2:271025/138‑1 (MQ=255)                                                                                                                                              gCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCAcc       >  1:271025/1‑138 (MQ=255)                                                                                                                                               |                                                                                                                                              GACAGCGTGAAGGTGCGTTGGGCAGTGGCAGACGCTTGCATTTATACCGAATGGGCCGAGAGCGACCTGAGCAAGTCGGCCGGGGCGCATTAACGCTCTTGTACGTTGCAGCGCCT‑GCCAGATCGAGCACCGCCCGCGCGGCACACGATCTGGCAGGCGCTGCACTATTGGAGGCGTACGCCCCGTCACATTATCATCACGCCACTGTCGTAATCTGCCGGCCAACTCTATGCCACAGGTCACCCGGCATGCTCCGCCTGCTTGTCTTGCTGCTCACCCTGCT  >  NC_002947/4344177‑4344459

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