Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	322,258	A→G	9.1%	A13A (GCA→GCG)	lytM →	M23 family lysostaphin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	322,258	0	A	G	9.1%	100.5 / 2.0	41	A13A (GCA→GCG)	lytM	M23 family lysostaphin
Reads supporting (aligned to +/- strand):  ref base A (15/22);  new base G (2/2);  total (17/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

AATATTAAATAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACTTATTATACAATT  >  CP000730/322119‑322396                                                                                                                                            |                                                                                                                                           aaTATTAAATAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAAcagcagcag                                                                                                                                           >  1:22854/1‑141 (MQ=255) aaTATTAAATAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAAcagcagcag                                                                                                                                           >  1:299435/1‑141 (MQ=255)       aaaTAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTg                                                                                                                                     >  2:318787/1‑141 (MQ=255)          tAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCaa                                                                                                                                  <  2:234521/141‑1 (MQ=255)              tGGAAAATGAGAATTAAATCCTGTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTAAACATGAAAAAATTAACAGCAGCAGCGATTGCAACGAt                                                                                                                              <  2:306874/141‑1 (MQ=255)                 aaaaTGAGAATTAAATCCTTTACTCATTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATggg                                                                                                                           <  2:76993/141‑1 (MQ=255)                             aaaTCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCGGCGATTGCAACGATGGGCTTCGCTACAtt                                                                                                               >  1:326577/1‑141 (MQ=255)                             aaaTCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCGGCGATTGCAACGATGGGCTTCGCTACAtt                                                                                                               >  1:280058/1‑141 (MQ=255)                              aaTCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACAttt                                                                                                              >  2:335743/1‑141 (MQ=255)                                tCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACCATTAc                                                                                                            >  2:65484/1‑141 (MQ=255)                                  cTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACaa                                                                                                          <  2:334616/141‑1 (MQ=255)                                        tCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGgcgc                                                                                                    <  1:119030/141‑1 (MQ=255)                                                cTAATTCTTTGAGTCTGTGCAGTACAGTTAAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAACa                                                                                            >  2:55038/1‑141 (MQ=255)                                                     ttttttAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATgc                                                                                       <  2:107309/139‑1 (MQ=255)                                                       ttttAGTATGTGCAGTAAAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATgcag                                                                                     <  2:177092/141‑1 (MQ=255)                                                        tttAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATgcagc                                                                                    >  1:37374/1‑141 (MQ=255)                                                         ttAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATgcagca                                                                                   <  1:359813/141‑1 (MQ=255)                                                          tAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATgcagcaa                                                                                  >  1:126830/1‑140 (MQ=255)                                                           attATGTGCAGTGCAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTAAACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGa                                                                                 <  2:96255/139‑1 (MQ=255)                                                           aGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGa                                                                                 >  2:342398/1‑141 (MQ=255)                                                                  gCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACa                                                                          <  2:302494/141‑1 (MQ=255)                                                                        cAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACAcc                                                                    >  2:75251/1‑141 (MQ=255)                                                                            ttAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCcaac                                                                <  2:337320/141‑1 (MQ=255)                                                                                   aaaaCTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGcaca                                                         <  1:211806/141‑1 (MQ=255)                                                                                      actaactTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATac                                                      <  1:90503/141‑1 (MQ=255)                                                                                              taactttaaTGGAGGATCTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGt                                              <  2:183326/141‑1 (MQ=255)                                                                                                     aaTGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCGGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAAcaca                                       <  2:326577/141‑1 (MQ=255)                                                                                                     aaTGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCGGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAAcaca                                       <  2:280058/141‑1 (MQ=255)                                                                                                           ggaTGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATcaca                                 <  2:328685/141‑1 (MQ=255)                                                                                                           ggaTGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATcaca                                 <  2:328658/141‑1 (MQ=255)                                                                                                             aTGTTGTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAg                               <  2:105087/141‑1 (MQ=255)                                                                                                              tggtttATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGa                              <  2:71895/138‑1 (MQ=255)                                                                                                              tGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGa                              <  2:186901/141‑1 (MQ=255)                                                                                                                ttttATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACg                            <  1:157434/141‑1 (MQ=255)                                                                                                                 tttATACATGAAAAAATTAACAGCATCAGCGATTCCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGt                           <  1:198520/141‑1 (MQ=255)                                                                                                                     taCATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCt                       >  1:44487/1‑141 (MQ=255)                                                                                                                       cATGAAAAAATTAACAGCAGCACCGATTGCAACGATGGGCTTCGCTACATTTACAATAGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGCCAACACAATCACAAGACGTATCTTa                     >  2:214992/1‑141 (MQ=255)                                                                                                                        aTGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTAt                    <  1:374610/141‑1 (MQ=255)                                                                                                                           aaaaaaTTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGt                 >  2:226131/1‑141 (MQ=255)                                                                                                                                aTTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACtta            >  1:78348/1‑141 (MQ=255)                                                                                                                                aTTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACtta            >  1:80575/1‑141 (MQ=255)                                                                                                                                   aaCAGCAGCAGCGATTGCAACGATGGGCTTCGCTCCATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACttatta         >  1:332582/1‑141 (MQ=255)                                                                                                                                    aCAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACAtattat        >  1:174610/1‑141 (MQ=255)                                                                                                                                          gcagcGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACAACAAACACCCAACAAGCACATACACAAAGGTCAACACAATCACAAGACGTATCTAATGGTACTTATTATACAAtt  >  2:23315/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                           AATATTAAATAGTTGGAAAATGAGAATTAAATCCTTTACTCAGTTGTCTAATTCTTTTAGTATGTGCAGTACAGTTTAATTGAAAACTAACTTTAACTTTAATGGAGGATGTTTTATACATGAAAAAATTAACAGCAGCAGCGATTGCAACGATGGGCTTCGCTACATTTACAATGGCGCATCAAGCAGATGCAGCAGAAACGACAAACACCCAACAAGCACATACACAAATGTCAACACAATCACAAGACGTATCTTATGGTACTTATTATACAATT  >  CP000730/322119‑322396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑