Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,373,213	A→G	Y122H (TAT→CAT)	nanA ←	N‑acetylneuraminate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,373,213	0	A	G	100.0%	122.6 / NA	34	Y122H (TAT→CAT)	nanA	N‑acetylneuraminate lyase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/18);  total (16/18)

GTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGA  >  NC_000913/3373081‑3373347                                                                                                                                     |                                                                                                                                       gTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCg                                                                                                                                 <  2:508068/140‑1 (MQ=255)   ttCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCa                                                                                                                               <  1:248340/140‑1 (MQ=255)       gcgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCa                                                                                                                               <  2:52839/136‑1 (MQ=255)       gcgcACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCa                                                                                                                               >  1:52839/1‑136 (MQ=255)                cGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGa                                                                                                                                    >  1:153245/1‑122 (MQ=255)                cGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGa                                                                                                                                    <  2:153245/122‑1 (MQ=255)                  ccAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCt                                                                                                                <  2:293772/140‑1 (MQ=255)                       cAATGTAACAAGTGTGTTGATCTGATCCAGGTTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAg                                                                                                           <  1:116187/140‑1 (MQ=255)                       cAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAg                                                                                                           <  1:127112/140‑1 (MQ=255)                       cAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAg                                                                                                           <  1:460145/140‑1 (MQ=255)                                    gtgtTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACgg                                                                                              >  1:104392/1‑140 (MQ=255)                                     tgtTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACggg                                                                                             >  2:76954/1‑139 (MQ=255)                                     tgtTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCACGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGc                                                                                             >  1:85084/1‑140 (MQ=255)                                          atctgatcCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGAc                                                                                        >  1:175568/1‑140 (MQ=255)                                            ctgatcCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGAcgg                                                                                      >  1:490891/1‑140 (MQ=255)                                                    gggTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACg                                                                              >  2:56337/1‑140 (MQ=255)                                                                ccccACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCa                                                                  <  2:104392/140‑1 (MQ=255)                                                                     cTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACg                                                             >  2:250861/1‑140 (MQ=255)                                                                          ggCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAg                                                        >  2:502390/1‑140 (MQ=255)                                                                           gCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGc                                                       >  2:312633/1‑140 (MQ=255)                                                                                aaTGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATg                                                  <  2:175568/140‑1 (MQ=255)                                                                                         caccacCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAgtt                                         >  2:251434/1‑140 (MQ=255)                                                                                            caccaTCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCg                                               <  1:423751/131‑1 (MQ=255)                                                                                            caccaTCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCg                                               >  2:423751/1‑131 (MQ=255)                                                                                              ccaTCGGCAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTgg                                    >  1:224957/1‑140 (MQ=255)                                                                                                     cAAACCATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTcgg                             <  1:439956/140‑1 (MQ=255)                                                                                                          cATCCGCCGAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTg                        <  2:224957/140‑1 (MQ=255)                                                                                                              cgccgAATCAATAATTGCCCGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGa                    <  1:227687/140‑1 (MQ=255)                                                                                                                               cccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACg         <  2:206810/134‑1 (MQ=255)                                                                                                                               cccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACg         >  1:206810/1‑134 (MQ=255)                                                                                                                               cccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCg   >  1:66470/1‑140 (MQ=255)                                                                                                                                ccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGa  <  1:251434/140‑1 (MQ=255)                                                                                                                                ccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGa  <  1:502390/140‑1 (MQ=255)                                                                                                                                ccGATGGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGa  <  2:66470/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       GTTTCAGCGCACCTACGCCAGGCAATGTAACAAGTGTGTTGATCTGATCCAGGGTCAGTTTTACCCCACTCAGGGCTGGAATGTTGTACACCACCATCGGCAAACCATCCGCCGAATCAATAATTGCCCGATAGTGATCGCAGTGTTCTTCAAAGCTGAAAGGATAGTAGAACGGCGTGACGGCGGAGACGGCATCGAAGCCATAACGTTTAGCCGATGCCGCAAGTTGTTGGCTTTCGGCGGTGCTGACGCAACCGACGTGGGCGA  >  NC_000913/3373081‑3373347

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