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%	95.3 / NA	30	*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/14);  total (16/14)

TTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGT  >  NC_000913/2233464‑2233727                                                                                                                                      |                                                                                                                                   ttAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATg                                                                                                                               >  2:581535/1‑139 (MQ=255)    aaCGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGAt                                                                                                                            >  1:631462/1‑139 (MQ=255)            ccGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcagca                                                                                                                    <  2:580686/139‑1 (MQ=255)                          gCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTg                                                                                                      <  1:59876/139‑1 (MQ=255)                              cTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGt                                                                                                  >  2:317879/1‑139 (MQ=255)                              cTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGt                                                                                                  >  2:550396/1‑139 (MQ=255)                                ggTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGttt                                                                                                <  2:367573/139‑1 (MQ=255)                                       aTTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCt                                                                                         <  1:581535/139‑1 (MQ=255)                                           cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATc                                                                                     >  2:129472/1‑139 (MQ=255)                                           cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATc                                                                                     <  2:426292/139‑1 (MQ=255)                                            ggCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATcc                                                                                    <  2:361711/139‑1 (MQ=255)                                             gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGAAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                   >  1:235733/1‑139 (MQ=255)                                                 tGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGata                                                                               >  2:259847/1‑139 (MQ=255)                                                            ccGGAAGATGCGCAGGGGTATTCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcg                                                                    <  2:469954/139‑1 (MQ=255)                                                            ccGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcg                                                                    >  1:797074/1‑139 (MQ=255)                                                              ggAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcgcg                                                                  <  2:724760/139‑1 (MQ=255)                                                                  gATGCGCAGGGGTATCCCGCCGTCACACCAGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                              >  1:589800/1‑139 (MQ=255)                                                                              tATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTg                                                  <  1:786188/139‑1 (MQ=255)                                                                                     ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCa                                           >  2:266371/1‑139 (MQ=255)                                                                                     ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCa                                           >  1:765/1‑139 (MQ=255)                                                                                     ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCa                                           >  1:376462/1‑139 (MQ=255)                                                                                          acacCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGgcgc                                      <  2:416923/139‑1 (MQ=255)                                                                                                              aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTaaaa                  >  1:725557/1‑139 (MQ=255)                                                                                                              aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTaaaa                  >  1:362343/1‑139 (MQ=255)                                                                                                                tACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATaaa                            >  1:699095/1‑127 (MQ=255)                                                                                                                tACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATaaa                            <  2:699095/127‑1 (MQ=255)                                                                                                                        aaaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGa        >  2:501156/1‑139 (MQ=255)                                                                                                                         aaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGaa       <  1:156827/139‑1 (MQ=255)                                                                                                                           cAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAAcc     <  2:797074/139‑1 (MQ=255)                                                                                                                              ggAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGgt  <  2:75245/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                   TTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGT  >  NC_000913/2233464‑2233727

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