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%	110.9 / NA	32	intergenic (‑32/‑278)	aspS/yecD	aspartyl‑tRNA synthetase/isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (16/16);  total (16/16)

ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTA  >  NC_000913/1950416‑1950682                                                                                                                                            |                                                                                                                                 aTGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATg                                                                                                                                   <  2:601619/139‑1 (MQ=255)       aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                            >  2:965513/1‑138 (MQ=255)                          gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATg                                                                                                                                   >  1:585264/1‑114 (MQ=255)                          gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATg                                                                                                                                   <  2:585264/114‑1 (MQ=255)                                      cccAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACGTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAg                                                                                              <  1:3853/139‑1 (MQ=255)                                         aaCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                           <  1:965513/139‑1 (MQ=255)                                          aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                        >  2:812276/1‑125 (MQ=255)                                          aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTc                                                                                                        <  1:812276/125‑1 (MQ=255)                                           ccACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                         >  2:1283537/1‑139 (MQ=255)                                           ccACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                         <  1:1283537/139‑1 (MQ=255)                                                caGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTca                                                                                    <  2:153834/139‑1 (MQ=255)                                                      cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtg                                                                              >  2:436763/1‑139 (MQ=255)                                                             tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                        >  2:841969/1‑90 (MQ=255)                                                             tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                        <  1:841969/90‑1 (MQ=255)                                                                       gggACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCt                                                             >  2:604000/1‑139 (MQ=255)                                                                          aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCa                                                                                                                    <  1:1400046/81‑1 (MQ=255)                                                                          aCAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCa                                                                                                                    >  2:1400046/1‑81 (MQ=255)                                                                              aCGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaa                                                                       >  1:1297276/1‑122 (MQ=255)                                                                              aCGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaa                                                                       <  2:1297276/122‑1 (MQ=255)                                                                                              tatTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               <  2:1294440/49‑1 (MQ=255)                                                                                              tatTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               >  1:1294440/1‑49 (MQ=255)                                                                                                   tGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTa                                 >  2:17023/1‑139 (MQ=255)                                                                                                          aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                   <  1:538139/66‑1 (MQ=255)                                                                                                          aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                   >  2:538139/1‑66 (MQ=255)                                                                                                          aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCt                          >  2:368979/1‑139 (MQ=255)                                                                                                             agaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCa                                                                                                     <  1:320689/61‑1 (MQ=255)                                                                                                             agaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCa                                                                                                     >  2:320689/1‑61 (MQ=255)                                                                                                             agaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGc                       <  1:368979/139‑1 (MQ=255)                                                                                                               atatCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                        <  2:137760/40‑1 (MQ=255)                                                                                                               atatCCCTTAACTTAGCTGCCGGCGGATGCCCCCctgctg                                                                                                                        >  1:137760/1‑40 (MQ=255)                                                                                                                            tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAAc           >  2:419240/1‑136 (MQ=255)                                                                                                                            tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAAc           <  1:419240/136‑1 (MQ=255)                                                                                                                              gCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTc      >  2:1059633/1‑139 (MQ=255)                                                                                                                               cTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCg     >  1:158624/1‑139 (MQ=255)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                   <  2:211588/42‑1 (MQ=37)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGc                                                                                                   >  1:211588/1‑42 (MQ=37)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa  <  2:158624/139‑1 (MQ=255)                                                                                                                                  ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa  <  1:17023/139‑1 (MQ=255)                                                                                                                                            |                                                                                                                                 ATGAAGATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTA  >  NC_000913/1950416‑1950682

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