Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,618,148	G→A	intergenic (+46/‑81)	deoC → / → deoA	deoxyribose‑phosphate aldolase/thymidine phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,618,148	0	G	A	100.0%	101.5 / NA	30	intergenic (+46/‑81)	deoC/deoA	deoxyribose‑phosphate aldolase/thymidine phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/16);  total (14/16)

TTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTG  >  NC_000913/4617920‑4618391                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ttCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCgg                                                                                                                                                                                                                                <  2:47487/250‑1 (MQ=255)          gctggcggcGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTa                                                                                                                                                                                                                      <  2:275248/248‑1 (MQ=255)                ggcgTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCCCTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACACGAAGGCCGGATAAGTCGTAGACGc                                                                                                                                                                                                                 >  1:102632/1‑250 (MQ=255)                 gcgtgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCAt                                                                                                                                                                                                               >  2:27244/1‑251 (MQ=255)                   gtgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATcc                                                                                                                                                                                                             <  2:246889/251‑1 (MQ=255)                    tgcgtACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCg                                                                                                                                                                                                            >  2:189683/1‑251 (MQ=255)                            gCGGAAGAAGCGCAGAAATATCTCGCCACTGAAGATGAAATGGTCGGTGCTGACTGGGAAGATGCGCGTCACTACCGCTTTGGCGCTGCAAGGCTGCAGGCAAGCATGAAGAAAGCGATGGGTCAGGGAGACGGTAAGAGCTCCAGCAGCTACTAAGTAAGATGCTTTACGCCGGATGCGCTGAGCTTATCAGGCTTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGAAGCATGCGGCAAAAgc                                                                                                                                                                                                    <  2:176347/251‑1 (MQ=255)                                          gAAATATCTAGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGCTGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCCTCCGGCAAAAGCCGCCTCATACTCtt                                                                                                                                                                                      >  2:124964/1‑251 (MQ=255)                                                   cGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCggg                                                                                                                                                                             <  2:196606/251‑1 (MQ=255)                                                         tGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGtt                                                                                                                                                                        <  1:89543/250‑1 (MQ=255)                                                             gATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCtt                                                                                                                                                                   <  1:218560/251‑1 (MQ=255)                                                               tGAACTGTTCGGTGCTGACCGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGACGCCTCATACTCTTTTCCTCGGGAGGATACCTTGt                                                                                                                                                                 >  2:167787/1‑251 (MQ=255)                                                                               gACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAat                                                                                                                                                 <  2:210671/251‑1 (MQ=255)                                                                                 cTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAatt                                                                                                                                                >  2:232170/1‑250 (MQ=255)                                                                                         aTGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACACGTAGGCCGGATAAGGCGTAGACGCATCAGGCAAAAGCCGCCTCATACTCTTTTCATCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTaa                                                                                                                                        >  1:14638/1‑250 (MQ=255)                                                                                               gTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACgg                                                                                                                                 >  1:173774/1‑250 (MQ=255)                                                                                                       cgcgttGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTATGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCa                                                                                                                         <  1:250122/247‑1 (MQ=255)                                                                                                       cGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGt                                                                                                                                                                                                                       >  2:109151/1‑157 (MQ=255)                                                                                                       cGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGt                                                                                                                                                                                                                       <  1:109151/157‑1 (MQ=255)                                                                                                           ttGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCGCCAGCTACTAAGTAAGATGCTTTAACCCTGACGCGCTGCCCTTATCAGCCCCACGAAACGTACCTCAGAGCACCCCGCATCACGCGTCGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCAAAAGAAATTATTCGTAAAAAACGTGATGGTCATg                                                                                                                       <  2:158634/249‑1 (MQ=255)                                                                                                                              ggCAAGCCTTCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAAtt                                                                                                   <  1:53090/250‑1 (MQ=255)                                                                                                                                     ctgctgAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCtt                                                                                           >  2:238735/1‑251 (MQ=255)                                                                                                                                           aaaGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCaa                                                                                     <  1:102455/251‑1 (MQ=255)                                                                                                                                               cgcTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGt                                                                                 >  2:58293/1‑251 (MQ=255)                                                                                                                                                                      agCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAGAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCg                                                           >  1:63822/1‑250 (MQ=255)                                                                                                                                                                                   actaAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAAACGCATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGACGCCTCATACTCTTTTCCTAGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTgcc                                             >  1:180021/1‑251 (MQ=255)                                                                                                                                                                                                               gctgcgctTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGACGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACg                 >  1:251131/1‑251 (MQ=255)                                                                                                                                                                                                               gctgcgctTATCAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACg                 <  1:280576/251‑1 (MQ=255)                                                                                                                                                                                                                          cAGGCCTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAAt       <  1:189683/250‑1 (MQ=255)                                                                                                                                                                                                                              ccTACGAAACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTg  <  1:238434/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    TTCAAACCGGCGGGCGGCGTGCGTACTGCGGAAGATGCGCAGAAATATCTCGCCATTGCAGATGAACTGTTCGGTGCTGACTGGGCAGATGCGCGTCACTACCGCTTTGGCGCTTCCAGCCTGCTGGCAAGCCTGCTGAAAGCGCTGGGTCACGGCGACGGTAAGAGCGCCAGCAGCTACTAAGTAAGATGCTTTACGCCTGATGCGCTGCGCTTATCAGGCCTACGAGACGTATCTACCCGTAGGCCGGATAAGGCGTAGACGCATCCGGCAAAAGCCGCCTCATACTCTTTTCCTCGGGAGGTTACCTTGTTTCTCGCACAAGAAATTATTCGTAAAAAACGTGATGGTCATGCGCTGAGCGATGAAGAAATTCGTTTCTTTATCAACGGTATTCGCGACAACACTATCTCCGAAGGGCAGATTGCCGCCCTCGCGATGACCATTTTCTTCCACGATATGACAATGCCTG  >  NC_000913/4617920‑4618391

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