Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	789,464	G→A	I172I (ATC→ATT)	galK ←	galactokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	789,464	0	G	A	100.0%	103.4 / NA	33	I172I (ATC→ATT)	galK	galactokinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/19);  total (14/19)

TACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAG  >  NC_000913/789331‑789583                                                                                                                                      |                                                                                                                        tACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATccc                                                                                                                    >  2:101078/1‑139 (MQ=255)  aCGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCg                                                                                                                   <  2:843898/139‑1 (MQ=255)         aaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                            <  1:104930/139‑1 (MQ=255)         aaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                            >  2:104930/1‑139 (MQ=255)                          gacgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGttt                                                                                           >  2:843456/1‑139 (MQ=255)                            cgacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTc                                                                                         >  2:755510/1‑139 (MQ=255)                                aGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCt                                                                                     >  2:137786/1‑139 (MQ=255)                                   cacaCCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTcc                                                                                  <  2:162310/139‑1 (MQ=255)                                           tGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGtt                                                                          >  2:77270/1‑139 (MQ=255)                                                  ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGcgcg                                                                   >  2:53823/1‑139 (MQ=255)                                                        tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACa                                                                                                    >  1:647613/1‑100 (MQ=255)                                                        tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACa                                                                                                    <  2:647613/100‑1 (MQ=255)                                                          cTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgc                                                           >  1:851273/1‑139 (MQ=255)                                                                tCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTc                                                     <  1:843456/139‑1 (MQ=255)                                                                     aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGtt                                                                                            >  1:178856/1‑95 (MQ=255)                                                                     aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGtt                                                                                            <  2:178856/95‑1 (MQ=255)                                                                     aGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCg                                                <  1:77270/139‑1 (MQ=255)                                                                        gAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCa                                             <  2:344859/139‑1 (MQ=255)                                                                                   gATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                  <  1:83988/139‑1 (MQ=255)                                                                                   gATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                  <  2:851273/139‑1 (MQ=255)                                                                                     tCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                <  2:309870/139‑1 (MQ=255)                                                                                        gCAAGGCATGATTTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGctgctg                             <  1:776707/139‑1 (MQ=255)                                                                                            ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                            <  1:493391/56‑1 (MQ=255)                                                                                            ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTAc                                                                                                            >  2:493391/1‑56 (MQ=255)                                                                                                       tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                >  2:887374/1‑121 (MQ=255)                                                                                                       tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctg                                <  1:887374/121‑1 (MQ=255)                                                                                                       tcttGCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGAc              >  2:394787/1‑139 (MQ=255)                                                                                                           gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg          <  1:665340/139‑1 (MQ=255)                                                                                                           gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg          <  1:53823/139‑1 (MQ=255)                                                                                                           gCCGAGCGCGGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgcg          <  1:137786/139‑1 (MQ=255)                                                                                                                 cgcgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                  >  2:733306/1‑109 (MQ=255)                                                                                                                 cgcgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgc                                  <  1:733306/109‑1 (MQ=255)                                                                                                                   cgGAAATTAGCTGATCCATAATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAg  <  1:141866/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                        TACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAG  >  NC_000913/789331‑789583

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