Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,073,600	T→C	8.1%	D219G (GAC→GGC)	PP_5327 ←	phosphate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,073,600	0	T	C	8.1%	138.5 / 2.8	49	D219G (GAC→GGC)	PP_5327	phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (22/23);  new base C (2/2);  total (24/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 = 9.86e-01

CACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTGGTTGGCACGCT  >  NC_002947/6073458‑6073747                                                                                                                                               |                                                                                                                                                      cACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGGGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑t                                                                                                                                                <  3:16549/149‑1 (MQ=255) cACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑t                                                                                                                                                >  5:145656/1‑149 (MQ=255)    tACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCgg                                                                                                                                             >  3:99969/1‑149 (MQ=255)             cAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGggcgg                                                                                                                                    >  1:113115/1‑149 (MQ=255)                gTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGCCACGGTCGGCGTCCATGATGGTCCTGggcggcg                                                                                                                                  >  1:111343/1‑149 (MQ=255)                    aTTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTc                                                                                                                             <  7:118412/149‑1 (MQ=255)                     tttCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTc                                                                                                                             >  4:196530/1‑148 (MQ=255)                              cACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGcc                                                                                                                    <  7:177812/148‑1 (MQ=255)                                cGGCCGTTTCCATCCCGGGCCCCCCGGCTGTGGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGGCGTGGATAGCCGTGTAGGGGGTGGTGAGTTGGGGCCGGTCGGCGTGCATG‑TCGGCCTGGGCGGCGGCGGCCCGCTTGCCGt‑c                                                                                                                >  4:233589/1‑149 (MQ=255)                                  gCCTTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCAATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGt‑cca                                                                                                              <  6:301711/149‑1 (MQ=255)                                         ccATCCCGGGCCACCAGGCTGGCGCGGGAAAAAGCCTGGTGCAGGCGGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGt‑ccagctccag                                                                                                       >  5:229362/1‑149 (MQ=255)                                           aTCCCGGGCCAGCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTc                                                                                                                             <  8:6341/126‑1 (MQ=255)                                           aTCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTc                                                                                                                             >  7:6341/1‑126 (MQ=255)                                           aTCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGCCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑ccagctccagct                                                                                                     >  4:57481/1‑149 (MQ=255)                                           aTCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGCCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑ccagctccagct                                                                                                     >  4:53133/1‑149 (MQ=255)                                                            tGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTgcaggcg                                                                                     >  6:9511/1‑148 (MQ=255)                                                                gcgGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTgcaggcgcaggc                                                                                >  2:147345/1‑149 (MQ=255)                                                                     aaaaaGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTc                                                                            >  2:244336/1‑148 (MQ=255)                                                                      aaaaGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCGCaa                                                                          <  6:145656/149‑1 (MQ=255)                                                                        aaGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAgg                                                                        <  2:113115/149‑1 (MQ=255)                                                                                 gCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGtt                                                                >  8:10977/1‑148 (MQ=255)                                                                                   aGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGAt                                                             <  1:147345/149‑1 (MQ=255)                                                                                     gCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATgg                                                           >  7:229261/1‑149 (MQ=255)                                                                                      ccGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGcc                                                                      >  5:95525/1‑137 (MQ=255)                                                                                      ccGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGcc                                                                      <  6:95525/137‑1 (MQ=255)                                                                                          tCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGt                                                                 <  2:49331/138‑1 (MQ=255)                                                                                          tCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGt                                                                 >  1:49331/1‑138 (MQ=255)                                                                                          tCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCg                                                      <  1:244336/149‑1 (MQ=255)                                                                                              gCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGt                                                  >  2:1348/1‑149 (MQ=255)                                                                                               cGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGcc                                                                                                                    >  5:141113/1‑83 (MQ=255)                                                                                               cGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGcc                                                                                                                    <  6:141113/83‑1 (MQ=255)                                                                                                ggTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTcc                                                <  3:196530/149‑1 (MQ=255)                                                                                                ggTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTcc                                                <  6:173414/149‑1 (MQ=255)                                                                                                  tCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCctt                                              <  5:9511/149‑1 (MQ=255)                                                                                                    ttCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATC‑TCGGACTGGGAGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCATGCGCTCAAGGCCATGGTTGATCGCACCGATGTCcttct                                            >  2:224650/1‑149 (MQ=255)                                                                                                      cGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCcttc                                             <  1:315135/146‑1 (MQ=255)                                                                                                      cGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGCCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCcttc                                             <  3:57481/146‑1 (MQ=255)                                                                                                      cGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGCCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCcttc                                             <  3:53133/146‑1 (MQ=255)                                                                                                            ccTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGTCCCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCAt                                                                    <  1:223633/118‑1 (MQ=255)                                                                                                            ccTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCAt                                                                    >  2:223633/1‑117 (MQ=255)                                                                                                            ccTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAgc                                    >  3:62184/1‑149 (MQ=255)                                                                                                                   gCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCtgctgc                              <  4:99969/148‑1 (MQ=255)                                                                                                                             aGTTCGCCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGtt                   <  2:111343/149‑1 (MQ=255)                                                                                                                              gTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTg                  <  7:10977/149‑1 (MQ=255)                                                                                                                                  ggCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCTCCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCt              <  1:224650/149‑1 (MQ=255)                                                                                                                                     aCGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTgca           <  3:233589/149‑3 (MQ=255)                                                                                                                                            gCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTGGTTGGCACg    >  2:23533/1‑149 (MQ=255)                                                                                                                                             cGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCTCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTGGTTGGCACGc   <  7:261972/149‑1 (MQ=255)                                                                                                                                              gTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTGGTTGGCACGCt  >  4:213157/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      CACTACCTTGCTCAGGTTGATTTCCACTTCACGGCCGTTGCCATCCCGGGCCACCAGGCTGTCGCGGGAAAAAGCCTGGTGCAGGCCGGTCAGGCGGTCTTCGATAGCCTTGTAGCGGCTGTTCAGTTCGGCACGGTCGGCGTCCATG‑TCGGCCTGGGCGGCGGCGTCCAGCTTGCCGT‑CCAGCTCCAGCTTGCGCGCTTGCAGGCGCAGGCGCTCAAGGCCATGGTTGATGGCACCGATGTCCTTCTTCTCCAGCTGCTGCAGCTCGCTGTTGAGCTGGTTGGCACGCT  >  NC_002947/6073458‑6073747

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