Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	500,543	C→G	10.0%	A364A (GCC→GCG)	PP_0411 →	polyamine ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	500,543	0	C	G	10.0%	147.9 / 6.9	50	A364A (GCC→GCG)	PP_0411	polyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base G (3/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

CGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAC  >  NC_002947/500396‑500688                                                                                                                                                    |                                                                                                                                                    cGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTCGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGccc                                                                                                                                                    >  5:233118/1‑149 (MQ=255)  gAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCcg                                                                                                                                                   >  7:300306/1‑149 (MQ=255)        aTTTACGTGGGCTAGCACGTGCGGGTGCGTCTGGAGGTTTGCGGCATAGGCGACTTCTTCGTGAAGCAGCGGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGAGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCt                                                                                                                                             <  8:199622/149‑1 (MQ=255)             ccTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑c                                                                                                                                       >  2:134126/1‑149 (MQ=255)                 ggCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTTGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGa                                                                                                                                    <  1:21315/148‑1 (MQ=255)                        aCGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCg                                                                                                                              >  3:344279/1‑147 (MQ=255)                                 tGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGa                                                                                                                  <  1:134126/149‑1 (MQ=255)                                  gCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGAt                                                                                                                 >  8:261084/1‑149 (MQ=255)                                         gAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGc                                                                                                          >  4:96956/1‑149 (MQ=255)                                           ggTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCtt                                                                                                        <  3:96956/149‑1 (MQ=255)                                            gTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCtt                                                                                                        >  3:81487/1‑148 (MQ=255)                                              ttGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGc                                                                                                                         >  6:99313/1‑129 (MQ=255)                                              ttGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGc                                                                                                                         <  5:99313/129‑1 (MQ=255)                                                gCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCAcc                                                                                                    <  8:300306/148‑1 (MQ=255)                                                    cAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAAccc                                                                                               >  8:22147/1‑149 (MQ=255)                                                         gCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTAc                                                                                          <  7:261084/149‑1 (MQ=255)                                                            aCTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGATCCCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACtg                                                                                        <  1:72790/149‑1 (MQ=255)                                                                    gTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgagaga                                                                                <  6:101843/148‑1 (MQ=255)                                                                    gTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaa                                                                               >  7:268923/1‑149 (MQ=255)                                                                     tGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaat                                                                              <  8:268923/149‑1 (MQ=255)                                                                          cagcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACtgt                                                                                       >  1:156110/1‑135 (MQ=255)                                                                          cagcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACtgt                                                                                       <  2:156110/135‑1 (MQ=255)                                                                            gcCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATg                                                                       <  4:81487/149‑1 (MQ=255)                                                                               gATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataataa                                                                         >  8:129616/1‑144 (MQ=255)                                                                               gATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataataa                                                                         <  7:129616/144‑1 (MQ=255)                                                                                 ttGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATgcgc                                                                    >  5:23649/1‑147 (MQ=255)                                                                                      gAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataa                                                                            >  2:242154/1‑134 (MQ=255)                                                                                      gAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGaataa                                                                            <  1:242154/134‑1 (MQ=255)                                                                                               ccGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGa                                                                                                                                    >  5:49035/1‑70 (MQ=255)                                                                                               ccGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGa                                                                                                                                    <  6:49035/70‑1 (MQ=255)                                                                                                cGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGAc                                                   >  6:144606/1‑149 (MQ=255)                                                                                                        gCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTTGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTg                                            <  7:22147/148‑1 (MQ=255)                                                                                                            tGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgag                                                                                   <  8:34598/105‑1 (MQ=255)                                                                                                            tGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGgag                                                                                   >  7:34598/1‑105 (MQ=255)                                                                                                              ggCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                      <  6:23649/148‑1 (MQ=255)                                                                                                               gCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTggg                                    >  1:206569/1‑149 (MQ=255)                                                                                                                     tagtaCCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGtt                              <  2:206569/149‑1 (MQ=255)                                                                                                                        taCCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCAc                           >  2:52962/1‑149 (MQ=255)                                                                                                                           cGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGc                        >  8:198282/1‑149 (MQ=255)                                                                                                                                ggctgggAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTggc                    <  6:233118/148‑1 (MQ=255)                                                                                                                                   tgggAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTggccatg                >  8:178392/1‑147 (MQ=255)                                                                                                                                         gTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATgg          <  4:344279/149‑1 (MQ=255)                                                                                                                                         atcgAGCGCGCGGGCGTGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTgggg                                   <  8:88294/121‑1 (MQ=255)                                                                                                                                         atcgAGCGCGCGGGCGTGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTgggg                                   >  7:88294/4‑124 (MQ=255)                                                                                                                                         atcgAGCGCGCGGGCGTGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                      >  3:43470/4‑121 (MQ=255)                                                                                                                                         atcgAGCGCGCGGGCGTGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTg                                      <  4:43470/118‑1 (MQ=255)                                                                                                                                           cgAGCGCGCGGGCGTGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGACCAGGCCATGGcc        >  1:119069/2‑149 (MQ=255)                                                                                                                                             aGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAtggcctggccag      <  7:178392/149‑3 (MQ=255)                                                                                                                                              gCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCg       >  2:169941/1‑147 (MQ=255)                                                                                                                                                 cGCCCGCGCGCTCGAT‑CCGATTGCCGATAGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGa   >  1:193045/1‑149 (MQ=255)                                                                                                                                                 cGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAc  <  1:169941/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGAATTCATTTACCTGGGCGACCACGTGCGGGTGCGTCTGGAGGTTTGCGGCAAAGGCGACTTCTTCGTGAAGCAGCCGATTGCCGAGCTCGACCCGGCGTTGGCCGTGGGCGATGTAGTACCGCTGGGCTGGGAGGTGGAGCACGCCCGCGCGCTCGAT‑CCGATTGCCGA‑AGCCCATTGATGGATTGCTTCACCAACCCTGTACTGTGGAGAGAATAATAATGCGCAAGCAGTTGAAACTGACCGCCCTGGCTTTGGGGCTGTTCACCGCTGGCCAGGCCATGGCCGCAGAC  >  NC_002947/500396‑500688

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