Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,088	Δ1 bp	coding (766/780 nt)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,088	0	G	.	100.0%	133.4 / NA	33	coding (766/780 nt)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (14/19);  total (14/19)

GATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCG  >  NC_000913/4617953‑4618218                                                                                                                                        |                                                                                                                                   gATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagc                                                                                                                               <  1:73052/140‑1 (MQ=255)  aTGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagca                                                                                                                              <  2:401784/140‑1 (MQ=255)   tGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagcag                                                                                                                             <  2:452321/140‑1 (MQ=255)   tGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagcag                                                                                                                             <  2:425242/140‑1 (MQ=255)   tGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagcag                                                                                                                             <  1:7077/140‑1 (MQ=255)   tGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑Ccagcag                                                                                                                             >  1:452321/1‑140 (MQ=255)    gcgcAGAGATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCagcagc                                                                                                                            >  2:82634/1‑140 (MQ=255)      gcAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGcta                                                                                                                          <  1:416929/140‑1 (MQ=255)                  cGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATg                                                                                                              >  2:578757/1‑140 (MQ=255)                    ccATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCt                                                                                                            >  1:288600/1‑140 (MQ=255)                        tGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGACAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTAc                                                                                                        <  1:481145/140‑1 (MQ=255)                         gCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACg                                                                                                       >  2:510825/1‑140 (MQ=255)                             aTGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTg                                                                                                   <  1:494321/140‑1 (MQ=255)                                 aCTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgc                                                                                               >  1:462856/1‑140 (MQ=255)                                  cTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcg                                                                                              <  1:578757/140‑1 (MQ=255)                                    gTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                                                                            >  1:70254/1‑140 (MQ=255)                                    gTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGAtgcgct                                                                                            >  2:10705/1‑140 (MQ=255)                                               aCTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAg                                                                                 >  1:463333/1‑140 (MQ=255)                                                  gggCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGcc                                                                              <  1:510825/140‑1 (MQ=255)                                                  gggCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGcc                                                                              <  2:233654/140‑1 (MQ=255)                                                    gCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTa                                                                            >  2:664927/1‑140 (MQ=255)                                                    gCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTa                                                                            >  1:102996/1‑140 (MQ=255)                                                    gCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTa                                                                            >  1:247639/1‑140 (MQ=255)                                                         tGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACgaga                                                                       <  1:82634/140‑1 (MQ=255)                                                         tGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACgaga                                                                       <  2:463333/140‑1 (MQ=255)                                                                          ttGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGc                                                      <  2:2286/140‑1 (MQ=255)                                                                                          gctgGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGAc                                      <  1:10705/140‑1 (MQ=255)                                                                                           ctgGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACg                                     >  2:549745/1‑140 (MQ=255)                                                                                           ctgGCAAGCCTGCTGAAAGCGCTGGGCCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACg                                     >  2:24147/1‑140 (MQ=255)                                                                                              gCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACCCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCAt                                  <  2:397569/140‑1 (MQ=255)                                                                                                      gctgAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCaaa                          >  1:342639/1‑140 (MQ=255)                                                                                                         gAAAGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgc                       <  2:462856/140‑1 (MQ=255)                                                                                                          aaaGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgcc                      <  1:558042/140‑1 (MQ=255)                                                                                                           aaGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccg                     <  2:102996/140‑1 (MQ=255)                                                                                                           aaGCGCTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccg                     <  1:664927/140‑1 (MQ=255)                                                                                                             gcgcTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAgccgcc                   <  1:516006/140‑1 (MQ=255)                                                                                                                cTGGGTCACGGCGACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCa                <  1:24147/140‑1 (MQ=255)                                                                                                                            gACGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCt    <  2:342639/140‑1 (MQ=255)                                                                                                                              cGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCg  <  1:181473/140‑1 (MQ=255)                                                                                                                              cGGTAAGAGC‑CCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCg  <  1:549745/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                   GATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCG  >  NC_000913/4617953‑4618218

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