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%	57.7 / NA	19	*240* (TAG→TAA)	sanA	DUF218 superfamily vancomycin high temperature exclusion protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/8);  total (11/8)

TTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACC  >  NC_000913/2233464‑2233724                                                                                                                                      |                                                                                                                                ttAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATg                                                                                                                            >  2:74840/1‑139 (MQ=255)                                     cTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTc                                                                                        <  1:72469/139‑1 (MQ=255)                                            ggCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATcc                                                                                 >  1:598237/1‑139 (MQ=255)                                             gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                >  2:632208/1‑139 (MQ=255)                                             gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                >  2:292998/1‑139 (MQ=255)                                             gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATAACAAAAGATACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                >  2:552857/1‑139 (MQ=255)                                             gCTATGCACCAGGTACCGGAAGATGCGCAGGGCTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATccc                                                                                >  2:458455/1‑139 (MQ=255)                                                              ggAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcgcg                                                               <  1:74840/139‑1 (MQ=255)                                                                                tCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATgg                                                                                                                           <  2:327733/61‑1 (MQ=255)                                                                                tCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATgg                                                                                                                           >  1:327733/1‑61 (MQ=255)                                                                                     ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcgcg                                                               >  2:386943/1‑116 (MQ=255)                                                                                     ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcgcg                                                               <  1:386943/116‑1 (MQ=255)                                                                                       gTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAgg                                      >  2:599017/1‑139 (MQ=255)                                                                                                    agctACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATaaa                         <  1:16727/136‑1 (MQ=255)                                                                                                      ttACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACg                       >  1:234443/1‑139 (MQ=255)                                                                                                              aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTaaaa               >  1:544072/1‑139 (MQ=255)                                                                                                                        aaaCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGa     <  2:690912/139‑1 (MQ=255)                                                                                                                          aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAAc   <  1:166825/139‑1 (MQ=255)                                                                                                                           cAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAAcc  <  1:552857/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                TTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACC  >  NC_000913/2233464‑2233724

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