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

GAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTG  >  NC_000913/2233472‑2233728                                                                                                                              |                                                                                                                                    gAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTAt                                                                                                                                 <  2:186404/130‑1 (MQ=255) gAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTAt                                                                                                                                 >  1:186404/1‑130 (MQ=255)  aaCCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcag                                                                                                                       <  2:1599065/139‑1 (MQ=255)                 ggCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGtt                                                                                                        >  1:980705/1‑139 (MQ=255)                                 tCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTa                                                                                        <  1:262490/139‑1 (MQ=255)                                   cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATc                                                                                      >  1:874864/1‑139 (MQ=255)                                   cGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATc                                                                                      >  2:1586668/1‑139 (MQ=255)                                     gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCtt                                                                                         >  2:1500402/1‑134 (MQ=255)                                     gCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCtt                                                                                         <  1:1500402/134‑1 (MQ=255)                                                  tACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTc                                                                                            >  1:1043510/1‑118 (MQ=255)                                                  tACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTc                                                                                            <  2:1043510/118‑1 (MQ=255)                                                   aCCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTcc                                                                      >  1:161104/1‑139 (MQ=255)                                                       gAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGGATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGa                                                                  >  2:1579757/1‑139 (MQ=255)                                                          gATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                               >  1:29857/1‑139 (MQ=255)                                                            tGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCAt                                                             >  1:341459/1‑139 (MQ=255)                                                             gcgcAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATg                                                            >  1:555276/1‑139 (MQ=255)                                                                  ggggTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAg                                                       >  2:1557918/1‑139 (MQ=255)                                                                           cgccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGgcgcg                                              <  2:147912/139‑1 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTAt                                                                                                                                 >  1:774971/1‑54 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTAt                                                                                                                                 <  2:774971/54‑1 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAca                                                                                                                        >  1:642770/1‑63 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAca                                                                                                                        <  2:642770/63‑1 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGtttt                                                                                                <  1:848082/87‑1 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGtttt                                                                                                >  2:848082/1‑87 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcg                                                                     >  2:87929/1‑114 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCcg                                                                     <  1:87929/114‑1 (MQ=255)                                                                             ccgTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAG‑TTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAg                                           >  1:868413/1‑139 (MQ=255)                                                                                 cacaCCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGgcg                                        <  1:833307/139‑1 (MQ=255)                                                                                            aGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATaaa                             <  2:29857/139‑1 (MQ=255)                                                                                            aGTTACTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATaaa                             <  1:1399883/139‑1 (MQ=255)                                                                                                aCTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcag                                                                                                                       >  2:1151016/1‑45 (MQ=255)                                                                                                aCTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTAcag                                                                                                                       <  1:1151016/45‑1 (MQ=255)                                                                                                aCTTGAATTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACggg                         >  2:451080/1‑139 (MQ=255)                                                                                                      aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                               >  2:1298628/1‑95 (MQ=255)                                                                                                      aTTACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGc                                                               <  1:1298628/95‑1 (MQ=255)                                                                                                        tACAAAAGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAAt                 >  2:109662/1‑139 (MQ=255)                                                                                                            aaaGAAACAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCAc             <  1:1579757/139‑1 (MQ=255)                                                                                                                  aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTc                                                                                            >  2:270718/1‑54 (MQ=255)                                                                                                                  aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTc                                                                                            <  1:270718/54‑1 (MQ=255)                                                                                                                  aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGgcgc                                       <  2:1132472/107‑1 (MQ=255)                                                                                                                  aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGgcgc                                       >  1:1132472/1‑107 (MQ=255)                                                                                                                  aCAAGGAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAAc       <  2:341459/139‑1 (MQ=255)                                                                                                                       gAAAGTAATTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGgtg  <  1:109662/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                    GAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTG  >  NC_000913/2233472‑2233728

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