Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,278,701	+C	coding (1060/1191 nt)	bcr ←	multidrug efflux pump Bcr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,278,697	1	.	C	100.0%	44.8 / NA	14	A355G (GCA→GGA)	bcr	multidrug efflux pump Bcr
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (12/2);  total (12/2)

TAGATCACCGTTTTTTCGGCCGACTGGCGTACAGACAGAAGAGAATGGAGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑G‑CCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGATGACACCATCGACACGCAGCCCACAAAC  >  NC_000913/2278566‑2278816                                                                                                                                      |                                                                                                                           tAGATCACCGTTTTTTCGGCCGACTGGCGTACAGACAGAAGAGAATGGAGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑G‑c                                                                                                                           >  2:12712/1‑133 (MQ=255)                                     agaagaGAATGGAGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAg                                                                                     >  1:9937/1‑133 (MQ=255)                                            aaTGGAGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGt                                                                              >  2:40904/1‑133 (MQ=255)                                             aTGGCGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGCCGCAAGAGAAAGCAATGCGCCAACAATT‑G‑CCCGCTAC‑GCC‑AAAACCGAAAGTTCCTGCAAGCGGACATGCCGtt                                                                             >  2:200227/1‑133 (MQ=255)                                                      ttggcgcAGAATGCAATTGACCAAATCATCGGCCATGCAGATTAAAGGGTCGCAAGAAAAAGCAATGCGCAAACAATTCC‑CCC‑CTAT‑GCA‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCata                                                                    <  2:365480/130‑1 (MQ=255)                                                      gTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCata                                                                    <  1:40904/133‑1 (MQ=255)                                                      gTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCC‑CCatat                                                                   >  1:6146/1‑133 (MQ=255)                                                                aTGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACt                                                          >  1:278929/1‑133 (MQ=255)                                                                 tGCAATTGACCAAATCATCGGCCATGCAGAGTTTAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACtc                                                         >  1:409709/1‑133 (MQ=255)                                                                         aCCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGaa                                                 >  1:304367/1‑133 (MQ=255)                                                                           cAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATg                                               >  1:338698/1‑133 (MQ=255)                                                                                  tcGGCCATGCAGATTAAAAGGTCGAAAGAAAAAGCAAGGCCCCAACAATTCC‑CCC‑CTAT‑CCA‑AAAACGGAAGGTTCCGGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCa                                        <  2:77801/133‑1 (MQ=255)                                                                                   cGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCAt                                       >  1:264587/1‑133 (MQ=255)                                                                                   cGGCAATGAAGAGTAAAAGGCCGAAAAAAAAAGAAAGGCCCAAACAATT‑GCCCC‑CTAT‑GCA‑AAAAGGAAAGTTCCCTCCAAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCAt                                       <  2:304367/133‑1 (MQ=255)                                                                                     gcaaTGAAGAGTAAAAGGTCGCAAGAGAAACAAATGCCCCAACATTTCC‑CCC‑CTGT‑ACA‑AAAGCGGAAGGTTCCGGCCAGCTAAAGGGCCTTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCg                                     <  2:480193/130‑1 (MQ=255)                                                                                             gagTAAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATTCC‑CCC‑CTAT‑CCA‑AAAACGGAAGGTCCCTGCCAGCGAAGATCCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGAt                             <  2:264587/133‑1 (MQ=255)                                                                                             gagTAAAAGGTCGCAAGAAAAACCAAGGC‑CCCACAACT‑G‑CCC‑CCCTAGCCAAAAACGGAAGTTTCCCCCCAGCGAAGAGGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGAt                             <  2:67760/133‑1 (MQ=255)                                                                                                        cGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGATGACACCATCGa                  >  2:309190/1‑133 (MQ=255)                                                                                                        cGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGATGACACCATCGa                  >  1:445985/1‑133 (MQ=255)                                                                                                        cGCAAGAGAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAAATCATCAAGAATGTCCGTCATCGCATCGGATCACACCATCGa                  >  2:147969/1‑133 (MQ=255)                                                                                                            agagAAAGCAATGCGCCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCC‑CCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGATGACACCATCGACACGc             <  2:6146/133‑1 (MQ=255)                                                                                                                        gcgcCAACAATT‑GCCCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTAGATGACAACATAGACACGCAGACCACAAAc  >  1:221593/1‑133 (MQ=255)                                                                                                                                      |                                                                                                                           TAGATCACCGTTTTTTCGGCCGACTGGCGTACAGACAGAAGAGAATGGAGCTGGTTGCGCAGAATGCAATTGACCAAATCATCGGCCATGCAGAGTTAAAGGTCGCAAGAGAAAGCAATGCGCCAACAATT‑G‑CCC‑CTAT‑GCC‑AAAACGGAAGGTTCCTGCCAGCGAAGATGCCGTTCCCGCCATATGGGGAAACTCATCAAGAATGACCGCCATCGCATTGGATGACACCATCGACACGCAGCCCACAAAC  >  NC_000913/2278566‑2278816

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