Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,233,597	G→A	*240* (TAG→TAA)	sanA →	DUF218 superfamily vancomycin high temperature exclusion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,233,597	0	G	A	100.0%	101.6 / NA	33	*240* (TAG→TAA)	sanA	DUF218 superfamily vancomycin high temperature exclusion protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

TTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCG  >  NC_000913/2233463‑2233725                                                                                                                                       |                                                                                                                                 tttAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTAt                                                                                                                              <  2:642807/139‑1 (MQ=255)  ttAAACGTGAACCGCG‑TTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATgg                                                                                                                            >  2:975320/1‑139 (MQ=255)           aaCCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcag                                                                                                                    >  2:911246/1‑139 (MQ=255)            aCCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcagc                                                                                                                   <  2:300982/139‑1 (MQ=255)                cgTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGtt                                                                                                               >  1:575711/1‑139 (MQ=255)                            ccGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGc                                                                                                   <  1:830302/139‑1 (MQ=255)                              gCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAAATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCg                                                                                                 <  1:147216/139‑1 (MQ=255)                                  gTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGtttt                                                                                             <  2:575711/139‑1 (MQ=255)                                   tCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGttttt                                                                                            <  2:655838/139‑1 (MQ=255)                                            cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATc                                                                                   <  1:373406/139‑1 (MQ=255)                                             ggCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATcc                                                                                  <  1:965091/139‑1 (MQ=255)                                             ggCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATcc                                                                                  <  1:656370/139‑1 (MQ=255)                                             ggCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAACAGATTTTTGTCGTTGCCGCTTTTTTCCTTATcc                                                                                  >  2:792569/1‑139 (MQ=255)                                              gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                 <  1:539399/139‑1 (MQ=255)                                                             ccGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcg                                                                  >  1:866212/1‑139 (MQ=255)                                                              cGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcgc                                                                 >  2:269472/1‑139 (MQ=255)                                                                 aaGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGaa                                                              >  1:889232/1‑139 (MQ=255)                                                                   gATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                            >  1:759516/1‑139 (MQ=255)                                                                       cgcAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATgg                                                        >  2:935678/1‑139 (MQ=255)                                                                                 tCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGgc                                              >  1:1053796/1‑139 (MQ=255)                                                                                 tCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGgc                                              >  1:1004257/1‑139 (MQ=255)                                                                                      ccgTCACACCCGAACCGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCCGGCTCGCa                                         >  2:363634/1‑139 (MQ=255)                                                                                          cacaCCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGgcg                                     <  1:419734/139‑1 (MQ=255)                                                                                            cacCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGcgcg                                   <  2:369961/139‑1 (MQ=255)                                                                                               ccGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGAtt                                <  1:975320/139‑1 (MQ=255)                                                                                                      gTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAAc                         <  2:759516/139‑1 (MQ=255)                                                                                                               aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTaaaa                >  1:159379/1‑139 (MQ=255)                                                                                                                    aaaaGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCa           <  2:1004257/139‑1 (MQ=255)                                                                                                                         aaaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                            <  2:913361/85‑1 (MQ=255)                                                                                                                         aaaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                            >  1:913361/1‑85 (MQ=255)                                                                                                                          aaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGaa     <  1:943742/139‑1 (MQ=255)                                                                                                                           aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAAc    >  1:128256/1‑139 (MQ=255)                                                                                                                             aaGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCg  >  2:192410/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                 TTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCG  >  NC_000913/2233463‑2233725

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