Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	524,114	Δ1 bp	10.3%	coding (1181/1200 nt)	tyrS →	tyrosine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	524,114	0	T	.	10.3%	110.7 / 7.8	39	coding (1181/1200 nt)	tyrS	tyrosine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (16/19);  new base . (2/2);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CTGGTCTGGTCAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCG  >  NC_002947/523967‑524251                                                                                                                                                    |                                                                                                                                           ctggtctggtCAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGtt                                                                                                                                           <  3:112896/149‑1 (MQ=255)    gtctggtcAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAATTCCTTTGGCCGGGTTAc                                                                                                                                         >  7:179107/1‑148 (MQ=255)        cgtcAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGAGCGGTGAAAGTTGATGGTGCAGTGGTCGACAATGACTTCATGTTTCTCCTAGGTGCTACCCATCACTGCCAAGCGGGCAAGAAGACGATTGGCCGGCTTACCCTc                                                                                                                                     >  2:5703/2‑148 (MQ=255)         gtcAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTc                                                                                                                                     >  2:299955/1‑147 (MQ=255)          tcAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAg                                                                                                                                  >  5:98885/1‑149 (MQ=255)            aaGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGc                                                                                                                                >  1:266549/1‑149 (MQ=255)              gAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGG‑TACCCTCAAGGCt                                                                                                                               >  5:74560/1‑147 (MQ=255)              gAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGG‑TACCCTCAAGGCTGa                                                                                                                             >  5:69162/1‑149 (MQ=255)               accTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTAGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGt                                                                                                                                            <  7:289338/132‑1 (MQ=255)               aaCTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCt                                                                                                                               >  5:346365/1‑147 (MQ=255)               aaCTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCAAAGCTGTCAAGAAGTCGTTTGGCCGGGt                                                                                                                                            >  8:289338/1‑134 (MQ=255)                  tCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGt                                                                                                                           >  6:231348/1‑148 (MQ=255)                    ggCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGG‑TACCCTCAAGGCTGAGTGAtt                                                                                                                       <  6:69162/149‑1 (MQ=255)                    ggCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGG‑TACCCTCAAGGCTGAGTGAtt                                                                                                                       <  6:74560/149‑1 (MQ=255)                     gCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGAtt                                                                                                                       <  8:126012/149‑1 (MQ=255)                           gcgcGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTg                                                                                                                  >  3:352700/1‑148 (MQ=255)                                  aCTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGAGCTACCCATGTGTGCCACGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGAGTGCCTGGCGGgtg                                                                                                           >  6:384849/1‑148 (MQ=255)                                       cTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCt                                                                                                      <  6:346365/148‑1 (MQ=255)                                          aaTAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAg                                                                                                                            <  5:322784/123‑1 (MQ=255)                                          aaTAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAg                                                                                                                            >  6:322784/1‑123 (MQ=255)                                          aaTAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCtata                                                                                                   >  4:161398/1‑148 (MQ=255)                                          aaTAGCGGCGCGGGGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCtat                                                                                                    <  8:206358/147‑1 (MQ=255)                                                   gcgGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCggg                                                                                         <  6:101474/149‑1 (MQ=255)                                                     ggTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGa                                                                                       <  5:231348/149‑1 (MQ=255)                                                          aaGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCTGGGAGGCCg                                                                                  <  4:352700/149‑1 (MQ=255)                                                          aaGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCg                                                                                  <  3:161126/149‑1 (MQ=255)                                                                aTGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGc                                                                            <  3:237367/149‑1 (MQ=255)                                                                 tGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAg                                                                                                        <  3:370283/120‑1 (MQ=255)                                                                 tGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAg                                                                                                        >  4:370283/1‑120 (MQ=255)                                                                              gACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGtt                                                             >  5:392726/1‑149 (MQ=255)                                                                                cAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGtttt                                                           <  1:299955/149‑1 (MQ=255)                                                                                             tttGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTtggtgg                                              <  6:98885/149‑1 (MQ=255)                                                                                                gcccTTGGGGCTACCCATGTGTGCCAAGCGGCAAATAAGTCGTTTGGCCGGGTTACCCTCACGGCTGCGTGATTCCGTGGATGGTGTAGCTAAATACACGCGGGGAGGCCGGTAGGCC‑GCCCCGCTTTGTTTTTGGGGTTTGGTGGtt                                            <  4:184237/146‑1 (MQ=255)                                                                                                       gTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCtgttgt                                    <  3:161398/149‑1 (MQ=255)                                                                                                             cccATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTaggggg                              <  2:266549/149‑1 (MQ=255)                                                                                                                     tgCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTTATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCt                      <  3:249721/149‑1 (MQ=255)                                                                                                                      gCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGTGGCCGGCCCCCGGTTGTTTTT‑TGGTTTGGTGGGTGTTTTTGTttggggggggggggg                     >  5:24625/1‑146 (MQ=255)                                                                                                                           gCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTg                                           <  3:64195/122‑1 (MQ=255)                                                                                                                              ggCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTgaga             <  6:392726/149‑1 (MQ=255)                                                                                                                                        cGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGGCGGTACGCC‑ACCCCCCATTGTTTTTTTGGTTTTTTGGTTGGGGTTGTTTGGGGAGGGGGCGGTTTTTAGAAGGTGTTTGc   >  6:149671/1‑149 (MQ=255)                                                                                                                                         gTTTGGCGGGGTTACCCTCAAGGCTGAGTGATTGCCTGGTTGGTGTAGCTATATAGAAGGGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCg  <  1:5703/149‑1 (MQ=255)                                                                                                                                         gTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGGGGGCGGGTGGGC‑ACCCCCCATTTGTTTTGGGGTTTGGTGGGTGGGGGGGGTGGGGGGGGGGGGGGGGTTGAGAGGGTGtttggg  >  8:171014/1‑145 (MQ=255)                                                                                                                                                    |                                                                                                                                           CTGGTCTGGTCAAGAACTCGGCACAGGCGCGAGACTTGCTGAATAGCGGCGCGGTGAAAGTTGATGGTGCAGTGGTCGACAAGGACTTCATGTTTGTCCTTGGTGCTACCCATGTGTGCCAAGCGGGCAAGAAGTCGTTTGGCCGGGTTACCCTCAAGGCTGAGTGATTGCCTGGATGGTGTAGCTATATAGAAGCGGGGAGGCCGGTAGGCC‑TCCCCGCTTTGTTTTTGGGGTTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCG  >  NC_002947/523967‑524251

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