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

ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTA  >  NC_000913/1950422‑1950682                                                                                                                                      |                                                                                                                                 aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                            >  1:847877/1‑138 (MQ=255) aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccc                                                                                                                            <  2:1282507/139‑2 (MQ=255) aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               <  2:1536980/136‑1 (MQ=255) aTCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               >  1:1536980/1‑136 (MQ=255)     ggCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctgct                                                                                                                        <  2:1648332/139‑6 (MQ=255)                    gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               <  2:1533396/117‑1 (MQ=255)                    gacgacgGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑ccc                                                                                                                               >  1:1533396/1‑117 (MQ=255)                                    aCCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGt                                                                                          <  1:968102/139‑1 (MQ=255)                                             agTCACCTGCTGCCCCACGTGGGAGAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                         >  1:1240697/1‑131 (MQ=255)                                             agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                 <  1:485999/139‑1 (MQ=255)                                             agTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                 <  2:1074031/139‑1 (MQ=255)                                             agTCACCTGCTGCCCCACGTGCGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTa                                                                                         <  2:1240697/131‑1 (MQ=255)                                                cACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACtg                                                                              >  2:969314/1‑139 (MQ=255)                                                 aCCTGCTACCACACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAATTCGCAGCGTTAGCTGTATGTCACAACtga                                                                             <  1:1531550/139‑1 (MQ=255)                                                       tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGaaa                                                                       >  1:762048/1‑139 (MQ=255)                                                       tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctg                                                                                                                          >  2:1003485/1‑84 (MQ=255)                                                       tgcCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑cccccctg                                                                                                                          <  1:1003485/87‑4 (MQ=255)                                                             aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCTGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                 <  2:1694067/139‑1 (MQ=255)                                                             aCGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                 <  1:969314/139‑1 (MQ=255)                                                              cGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                 >  2:163437/1‑106 (MQ=255)                                                              cGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGt                                                                                                 <  1:163437/106‑1 (MQ=255)                                                                tGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggcggc                                                              >  2:195054/1‑139 (MQ=255)                                                                   gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                           <  1:100203/107‑1 (MQ=255)                                                                   gACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                           >  2:100203/1‑107 (MQ=255)                                                                      aaaCGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCtattat                                                        <  1:195054/139‑1 (MQ=255)                                                                        aaggAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAg                                                                                                                   <  1:1314533/76‑1 (MQ=255)                                                                        aCGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAg                                                                                                                   >  2:1314533/1‑78 (MQ=255)                                                                            aGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGa                                                  >  1:1226676/1‑139 (MQ=255)                                                                                                  gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                      >  1:763127/1‑65 (MQ=255)                                                                                                  gCATGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGc                            >  2:136019/1‑139 (MQ=255)                                                                                                  gCATGAGATATCCCTAAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                      <  2:763127/65‑1 (MQ=255)                                                                                                    aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                 >  2:523573/1‑100 (MQ=255)                                                                                                    aTGAGATATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAggc                                                                 <  1:523573/100‑1 (MQ=255)                                                                                                        gaTATCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTACGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTACTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGcc                      <  2:1226676/139‑1 (MQ=255)                                                                                                           atCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGTTATTATACTGGAAATTCTGCCGCa                                      >  2:641949/1‑120 (MQ=255)                                                                                                           atCCCTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGTTATTATACTGGAAATTCTGCCGCa                                      <  1:641949/120‑1 (MQ=255)                                                                                                              ccTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                            <  1:759389/95‑1 (MQ=255)                                                                                                              ccTTAACTTAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCta                                                            >  2:759389/1‑95 (MQ=255)                                                                                                                      tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                      <  1:1390835/45‑1 (MQ=38)                                                                                                                      tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGc                                                                                                      >  2:1390835/1‑45 (MQ=38)                                                                                                                      tAGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGcc        >  1:388318/1‑139 (MQ=255)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                           >  2:1726021/1‑55 (MQ=39)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTg                                                                                           <  1:1726021/55‑1 (MQ=39)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                 <  2:168878/65‑1 (MQ=255)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACaa                                                                                 >  1:168878/1‑65 (MQ=255)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGaaa                                                >  2:1607683/1‑98 (MQ=255)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGaaa                                                <  1:1607683/98‑1 (MQ=255)                                                                                                                       aGCTGCCGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCt       >  1:844811/1‑139 (MQ=255)                                                                                                                            ccGGCGGATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTa  <  2:388318/139‑1 (MQ=255)                                                                                                                              ggcggATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCcg                    >  1:857819/1‑119 (MQ=255)                                                                                                                              ggcggATGCCCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCcg                    <  2:857819/119‑1 (MQ=255)                                                                                                                                      |                                                                                                                                 ATCAGGCTACCAAGATCACGACGACGGTTGACCCAACCACACAGAGTCACCTGCTGCCCCACGTGGGACAAACGGAGCTGTCCACAATATTCTGTACGCATGAGATATCCCTTAACTTAGCTGCCGGCGGATG‑CCCCCTGCTGCGCAGGTGACCAAGTCGCAGCGTTAGCTGTATGTCACAACTGAATGAAAAAAGGCGGCTATTATACTGGAAATTCTGCCGCACCGTAAGAGCCTGGCCCGCGCTGGAACGCCTCGTTA  >  NC_000913/1950422‑1950682

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