Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,950,559	(C)5→6	intergenic (‑37/‑273)	aspS ← / → yecD	aspartyl‑tRNA synthetase/isochorismatase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,950,554	1	.	C	100.0%	106.1 / NA	31	intergenic (‑32/‑278)	aspS/yecD	aspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (14/17);  total (14/17)

ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950416‑1950683                                                                                                                                            |                                                                                                                                  aTGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATg                                                                                                                                    <  1:97462/139‑1 (MQ=255)       aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                             >  1:26633/1‑138 (MQ=255)            gCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                         >  1:158913/1‑139 (MQ=255)            gCTACCAAGATCACGACGACGGTTAACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTACACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGTATGCCCCCctgctg                                                                                                                         >  2:489504/1‑139 (MQ=255)                          gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:723017/1‑139 (MQ=255)                          gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAg                                                                                                           >  1:67015/1‑139 (MQ=255)                               cgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCa                                                                                                      <  2:275801/139‑1 (MQ=255)                                ggTTGACCCAACCACGCAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTCAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAg                                                                                                     >  1:535058/1‑139 (MQ=255)                                         aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                            <  1:488688/139‑1 (MQ=255)                                              cacaGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGt                                                                                       >  1:160331/1‑139 (MQ=255)                                                   agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                  <  2:171091/139‑1 (MQ=255)                                                      cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtg                                                                               >  2:457127/1‑139 (MQ=255)                                                        ccTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtgaa                                                                             >  2:65011/1‑139 (MQ=255)                                                                        ggACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             <  2:26633/139‑1 (MQ=255)                                                                        ggACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                             <  1:88333/139‑1 (MQ=255)                                                                          aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtatt                                                           <  2:765972/139‑1 (MQ=255)                                                                          aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtatt                                                           <  2:419857/139‑1 (MQ=255)                                                                                   gCTGTCCACAATATTCTGTCCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGaa                                                  >  1:655167/1‑139 (MQ=255)                                                                                   gCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtattat                                                         >  2:319023/1‑132 (MQ=255)                                                                                   gCTGACCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtattat                                                         <  1:319023/132‑1 (MQ=255)                                                                                         cacaATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTg                                            <  1:453024/139‑1 (MQ=255)                                                                                                 tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGt                                                                                           <  1:419318/84‑1 (MQ=255)                                                                                                 tCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGt                                                                                           >  2:419318/1‑84 (MQ=255)                                                                                                   tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  2:103873/139‑1 (MQ=255)                                                                                                   tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                  <  1:560994/139‑1 (MQ=255)                                                                                                      aCGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAaga                               >  1:2345/1‑139 (MQ=255)                                                                                                          aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                           >  2:702472/1‑139 (MQ=255)                                                                                                             agaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGATAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTCGAAATTCTGCCGCACCGTAAGAGCCTGGc                        <  2:535058/139‑1 (MQ=255)                                                                                                                       tAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGa              >  1:645289/1‑139 (MQ=255)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  2:160331/139‑1 (MQ=255)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa   <  2:723017/139‑1 (MQ=255)                                                                                                                                   cggcggATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  <  2:2345/139‑1 (MQ=255)                                                                                                                                   cggcggATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAc  <  1:762028/139‑1 (MQ=255)                                                                                                                                            |                                                                                                                                  ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTAC  >  NC_000913/1950416‑1950683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑