Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,171,616	Δ1 bp	9.3%	coding (54/828 nt)	PP_1026 ←	M23/M37 family peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,171,616	0	G	.	9.3%	132.1 / 7.4	43	coding (54/828 nt)	PP_1026	M23/M37 family peptidase
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base . (2/2);  total (21/22)
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.81e-01

TTGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAGGCGGGCGGTCAGG  >  NC_002947/1171473‑1171760                                                                                                                                                 |                                                                                                                                                 ttGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGAGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGc                                                                                                                                             <  1:179469/149‑1 (MQ=255)      cTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCcag                                                                                                                                          >  2:110128/1‑147 (MQ=255)         cTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCcagcagca                                                                                                                                     <  5:176853/149‑1 (MQ=255)                          ccGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAg                                                                                                                     >  6:114639/1‑148 (MQ=255)                          ccGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAg                                                                                                                     >  3:322362/1‑148 (MQ=255)                                  ccGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAgg                                                                                                                          <  1:339409/135‑1 (MQ=255)                                  ccGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAgg                                                                                                                          >  2:339409/1‑135 (MQ=255)                                    gTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAg                                                                                                                           >  2:24453/1‑132 (MQ=255)                                       gAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGACTG‑GGCGCCGCCGGCCAGAAGCAGAAGGGAATGGGCAAGCATGGGAGCGAGCa                                                                                                       >  6:122777/1‑149 (MQ=255)                                           cGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCTGGTGATGTAGCTGGCCTT‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCAcgcg                                                                                                   <  5:122777/149‑1 (MQ=255)                                             ggCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCggg                                                                                                 <  5:114639/149‑1 (MQ=255)                                             ggCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCggg                                                                                                 <  4:322362/149‑1 (MQ=255)                                                gcTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCa                                                                                              <  7:119348/149‑1 (MQ=255)                                                             gcggcgGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCAt                                                                                             >  2:41504/1‑137 (MQ=255)                                                             gcggcgGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCAt                                                                                             <  1:41504/137‑1 (MQ=255)                                                              cggcggGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCagt                                                                                <  1:264655/149‑1 (MQ=255)                                                                   gggCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑CGCGCCGCCGGCCAGCAGCAGCAGGGAAAGTGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTga                                                                           <  1:110128/149‑1 (MQ=255)                                                                       ccAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                                    <  2:7517/136‑1 (MQ=255)                                                                       ccAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                                    >  1:7517/1‑136 (MQ=255)                                                                       ccAGGTCGACGACTGCCACACCGCCGGGCACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAgg                                                                       >  7:142727/1‑149 (MQ=255)                                                                                     gCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGa                                                                                                          >  4:193930/1‑100 (MQ=255)                                                                                            cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCaa                                                                                                                      >  2:221018/1‑81 (MQ=255)                                                                                            cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCaa                                                                                                                      <  1:221018/81‑1 (MQ=255)                                                                                            cgccgGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTg                                                  >  3:259924/1‑149 (MQ=255)                                                                                                 ggcaCCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTc                                             <  8:142727/146‑1 (MQ=255)                                                                                                      cGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGAccc                                        >  6:47993/1‑149 (MQ=255)                                                                                                             ttCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAg                                 >  5:270109/1‑149 (MQ=255)                                                                                                               cAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGc                                <  6:270109/148‑1 (MQ=255)                                                                                                                aGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCaa                              >  2:315970/1‑149 (MQ=255)                                                                                                                  tgtgCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                                    <  6:13375/93‑1 (MQ=255)                                                                                                                  tgtgCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTc                                                                                    >  5:13375/1‑93 (MQ=255)                                                                                                                       gggTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCt                       <  5:47993/149‑1 (MQ=255)                                                                                                                              gTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAggg                                                                                                                <  4:101531/53‑1 (MQ=255)                                                                                                                              gTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAggg                                                                                                                >  3:101531/1‑53 (MQ=255)                                                                                                                               ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc                     <  4:8897/141‑1 (MQ=255)                                                                                                                               ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc                     >  7:239645/3‑143 (MQ=255)                                                                                                                               ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc                     <  8:239645/141‑1 (MQ=255)                                                                                                                               ctgCTGGCCGGCGGCGCC‑CAGGCCAGCAGCAGCAGGGAAAGCGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTc                     >  3:8897/3‑143 (MQ=255)                                                                                                                                 gCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCACGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAggc             >  4:205833/1‑149 (MQ=255)                                                                                                                                  cTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAggcg            <  4:259924/149‑1 (MQ=255)                                                                                                                                   tGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCggg                                                              >  2:296831/1‑98 (MQ=255)                                                                                                                                   tGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCggg                                                              <  1:296831/98‑1 (MQ=255)                                                                                                                                           ggCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCACGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAGGCGGGCGGt      <  3:205833/147‑1 (MQ=255)                                                                                                                                             cgccgcCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCCCCAGCAACTCGCCTTCGTTCAGGCGGGCGGTCAgg  >  3:42206/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 TTGTCCTGCTCCTTCACCACCAGCACCGGCTTGCCGTCGAAACGGGCGCTGGGTGCATTGGCGGCGGGGCCCAGGTCGACGACTGCCACACCGCCGGGTACCGGTTTGTTCAGTGTGCGGGTGATGTAGCTGGCCTG‑GGCGCCGCCGGCCAGCAGCAGCAGGGAAAGGGCAAGCAGGGGCGCGAGCACGCGGGGCATGTGTCAGTCCAGTAGTGAGAGGGTGACCGGGGTCAGGTGGTTGTCTTCGACCCGCACCAGCAACTCGCCTTCGTTCAGGCGGGCGGTCAGG  >  NC_002947/1171473‑1171760

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