Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	505,779	A→C	H289P (CAT→CCT)	ushA →	5'‑nucleotidase/UDP‑sugar hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	505,779	0	A	C	100.0%	91.5 / NA	30	H289P (CAT→CCT)	ushA	5'‑nucleotidase/UDP‑sugar hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (15/15);  total (15/15)

CCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCATGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGGTCG  >  NC_000913/505538‑506022                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ccAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTgggg                                                                                                                                                                                                                                            <  1:212419/251‑1 (MQ=255)               cgcgGCGACCCATCTGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGAc                                                                                                                                                                                                                              >  1:53095/1‑251 (MQ=255)                         cATATGGGGCATTACGCTAATGGTGAGCACGGCTCTAACGCACCGGGCGCTGTGGAGCTGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGCTTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGAttt                                                                                                                                                                                                                    >  2:78854/1‑251 (MQ=255)                              ggggCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGCTCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCCACCCAGATCAACAAAACGGCATCTGGATTGTGCAGGAGCCTGGGGGGGGAAAATACGTGGGACAGGCTGATTTTgggt                                                                                                                                                                                                               >  2:178386/1‑248 (MQ=255)                               gggCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGGCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGAGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGGTGAATTTGAGtt                                                                                                                                                                                                              >  1:60228/1‑251 (MQ=255)                                      aCGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGGGGGACGGGCTGATTTTGAGTTTCGTAAt                                                                                                                                                                                                       <  2:224866/251‑1 (MQ=255)                                         aTAATGGTGTGCCCGTCTCGAACGCACCGGGAGATTTGTAGATGGCACGCGCGCTGCCTGCCGGATCGCGGGCGATGATCGTCGGTGGTCACTCTCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTTCGCCATTCAAACCAGATCAACAAAACGGCAGCTGTATTGTGCAGGCGCCTGAGTGGGGCAAATACTTGGGAAGGGCTGATTTTGAGTTTCGTAATGGc                                                                                                                                                                                                    <  2:230424/251‑1 (MQ=255)                                              ggTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAAt                                                                                                                                                                                               <  2:53095/251‑1 (MQ=255)                                               gTGAGCACGGCTATAACGCACAGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATg                                                                                                                                                                                              <  2:127506/251‑1 (MQ=255)                                                      cGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATgg                                                                                                                                                                                       >  1:254400/1‑251 (MQ=255)                                                           ctAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGTAATGAAAATGGTTa                                                                                                                                                                                    <  1:79882/249‑1 (MQ=255)                                                                                            cgcTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTgaagaagaa                                                                                                                                                  >  1:198578/1‑250 (MQ=255)                                                                                                    gccGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTAGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTTACCTg                                                                                                                                         <  1:36703/251‑1 (MQ=255)                                                                                                        gATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATCACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGaa                                                                                                                                     <  1:178386/251‑1 (MQ=255)                                                                                                            gCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGACAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGCCGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGATAACTACCCGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACg                                                                                                                                 >  2:119288/1‑251 (MQ=255)                                                                                                                  gatgaaCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGaaaa                                                                                                                           <  2:257956/251‑1 (MQ=255)                                                                                                                   atgatCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAag                                                                                                                          <  1:78854/251‑1 (MQ=255)                                                                                                                   atgatCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAag                                                                                                                          <  2:131928/251‑1 (MQ=255)                                                                                                                            ggtggtCACTCGCAAGATCAGGTCTGCATGGCGGCAGAAAACAAAAACCAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCAACTGGATTGTGAAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGt                                                                                                                 >  1:225909/1‑251 (MQ=255)                                                                                                                                         aaGATCCGGTATGCATGGCGGCAGAAAACAAAAAACAGGTCGATAACGTGCCGGGTACGCCATGCAAACCAGATCAACAACACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATAAGTGGGACGGGCTGATTTTGAGTTTCGTAAAGGCGAAATGAACATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCAGGGAAGACGGGAAAagcgagcgagagattaacactcat                                                                                                    >  1:157017/1‑235 (MQ=255)                                                                                                                                          aGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAAAAAAACGGCATCTGGATTGTGCAGGCGCCTGAGGGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGAAATGGCGAAATGAAAATGGTTACCTACCAGCTGATTCCGGTGAAACTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGAGCGTGCTTTACACTCCTg                                                                                                   >  1:143656/1‑251 (MQ=255)                                                                                                                                                     gCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCt                                                                                           <  1:53807/248‑1 (MQ=255)                                                                                                                                                             gcAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGc                                                                                 >  1:39840/1‑250 (MQ=255)                                                                                                                                                              cAGAAAACAACAAACAGGCCGATTACGGGCCGGGAACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGCGAACCTGAAGAAGAAAGTGGCCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCa                                                                                >  2:51849/1‑250 (MQ=255)                                                                                                                                                                         aaaCAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCt                                                                    <  1:242206/251‑1 (MQ=255)                                                                                                                                                                                 cGATTACGAGCCGGGTACGCCATGCAAACCAGCTCAACAAACCGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGGCGGGCTGCTTTTGAGTTTCGGAATGGCGAAATGAAAATGGATAACTACCAGCTGATTCCGGTGACCCTGATGACGCAAGTGACCTGGGAAGAAGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAACTCGCTGCAAACCAGCAAATGATCTCGCTGTTATCAc                                                            >  2:253367/1‑251 (MQ=255)                                                                                                                                                                                     tACGAGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTCATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGCAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATAGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGtt                                                        >  2:215348/1‑251 (MQ=255)                                                                                                                                                                                                   gCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGc                                           >  1:205659/1‑250 (MQ=255)                                                                                                                                                                                                        gCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCCTGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAgc                                      <  2:73465/250‑1 (MQ=255)                                                                                                                                                                                                                                            aGGCGCCTGAGTGGGGCAACTACGTGGGACGGGCTGATTTTGAGCTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGCCCTGGGAAGACGGGAAAAGCGATCGCGTGCTTTGCACTCCTGAAATCGCTGAAAACCAGCAAATGAGCGCGCTGTTATCCCCGTTCCAGAGCAAAGGCAAAGCGCAGCTGGAAGTGAAGACAGGCGAATCCaaggggcg  <  1:207395/250‑7 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    CCAGACATTATTATCGCGGCGACCCATATGGGGCATTACGATAATGGTGAGCACGGCTCTAACGCACCGGGCGATGTGGAGATGGCACGCGCGCTGCCTGCCGGATCGCTGGCGATGATCGTCGGTGGTCACTCGCAAGATCCGGTCTGCATGGCGGCAGAAAACAAAAAACAGGTCGATTACGTGCCGGGTACGCCATGCAAACCAGATCAACAAAACGGCATCTGGATTGTGCAGGCGCATGAGTGGGGCAAATACGTGGGACGGGCTGATTTTGAGTTTCGTAATGGCGAAATGAAAATGGTTAACTACCAGCTGATTCCGGTGAACCTGAAGAAGAAAGTGACCTGGGAAGACGGGAAAAGCGAGCGCGTGCTTTACACTCCTGAAATCGCTGAAAACCAGCAAATGATCTCGCTGTTATCACCGTTCCAGAACAAAGGCAAAGCGCAGCTGGAAGTGAAAATAGGCGAAACCAATGGTCG  >  NC_000913/505538‑506022

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