Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,386,861	A→G	L443P (CTG→CCG)	aaeB ←	p‑hydroxybenzoic acid efflux system component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,386,861	0	A	G	100.0%	180.8 / NA	50	L443P (CTG→CCG)	aaeB	p‑hydroxybenzoic acid efflux system component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/23);  total (27/23)

CCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGCGTGGATTCGGCAAAC  >  NC_000913/3386714‑3387002                                                                                                                                                    |                                                                                                                                              ccGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCa                                                                                                                                            <  1:495179/151‑1 (MQ=255)  cgacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCAc                                                                                                                                           <  1:154167/151‑1 (MQ=255)   gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACt                                                                                                                                          <  1:200180/151‑1 (MQ=255)   gacgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACt                                                                                                                                          <  1:713690/151‑1 (MQ=255)     cgaTTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGc                                                                                                                                        >  2:718982/1‑151 (MQ=255)         ttGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAgg                                                                                                                                    >  1:632668/1‑151 (MQ=255)              cTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAAt                                                                                                                               >  1:359647/1‑151 (MQ=255)               tAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATg                                                                                                                              >  2:256527/1‑151 (MQ=255)                   gcgcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGcac                                                                                                                           >  1:437107/1‑150 (MQ=255)                     gcTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCAcagc                                                                                                                        <  1:3972/151‑1 (MQ=255)                                 ctgactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGctc                                                                                                            >  2:658336/1‑151 (MQ=255)                                  tgactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGctct                                                                                                           >  2:352743/1‑151 (MQ=255)                                   gactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGctct                                                                                                           <  2:128778/150‑1 (MQ=255)                                   gactgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGctct                                                                                                           <  1:526888/150‑1 (MQ=255)                                      tgaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTg                                                                                                       >  1:588021/1‑151 (MQ=255)                                       gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGtt                                                                                                        <  2:205629/149‑1 (MQ=255)                                       gaAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTgg                                                                                                      >  1:719512/1‑151 (MQ=255)                                             gAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTa                                                                                                <  1:81909/151‑1 (MQ=255)                                               aaGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAgg                                                                                              >  1:187364/1‑151 (MQ=255)                                                aGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                                                             >  1:649780/1‑151 (MQ=255)                                                      cGGGTTATCCAGCCCGATAATATTTATGGTGCTGGCCAGTGCCCACATCGAGCCCAGTCGCCGTTTCTGTCCTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAggg                                                                                             <  2:473506/145‑1 (MQ=255)                                                        ggTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCAc                                                                                     <  1:677244/151‑1 (MQ=255)                                                            aTCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCaa                                                                                  >  1:95424/1‑150 (MQ=255)                                                                   aCGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTaga                                                                          <  2:376450/151‑1 (MQ=255)                                                                   aCGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTag                                                                           <  2:58489/150‑1 (MQ=255)                                                                       taataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagc                                                                        >  2:168480/1‑149 (MQ=255)                                                                        aataTTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGagcagc                                                                     >  2:78715/1‑151 (MQ=255)                                                                           atTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGccc                                                                   <  1:240571/150‑1 (MQ=255)                                                                               aTGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCGAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTaa                                                                <  1:720616/149‑1 (MQ=255)                                                                                  gTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGc                                                            >  2:268862/1‑150 (MQ=255)                                                                                       ggCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGcgcg                                                      <  1:573748/151‑1 (MQ=255)                                                                                        gCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAgcgc                                                       <  2:514125/149‑1 (MQ=255)                                                                                         ccAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCGACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGc                                                    >  2:605096/1‑151 (MQ=255)                                                                                          cAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGcc                                                   >  1:704492/1‑151 (MQ=255)                                                                                           aGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGcc                                                   >  2:249730/1‑150 (MQ=255)                                                                                            gTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAg                                                 <  1:611796/151‑1 (MQ=255)                                                                                              gCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCg                                               >  1:398302/1‑151 (MQ=255)                                                                                                  ccATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTccc                                           >  1:468031/1‑151 (MQ=255)                                                                                                   cATCGAGCCCAGTCGCCGTTTCTGTACTTCTCTACCGAGGAAGAATCCCGGCACTGCCAGGCCAGTGCACAGCAGCAACCTGCTCTGTTGGGCATTAGGGATAATCACCACAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCg                                          >  2:681605/1‑151 (MQ=255)                                                                                                    aTCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCg                                          >  2:410382/1‑150 (MQ=255)                                                                                                       gAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGt                                         >  2:31437/1‑148 (MQ=255)                                                                                                       gAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGt                                         >  2:327370/1‑148 (MQ=255)                                                                                                              gTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTc                               <  2:382004/151‑1 (MQ=255)                                                                                                                cGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCg                              >  1:593376/1‑150 (MQ=255)                                                                                                                cGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCg                              <  1:636797/150‑1 (MQ=255)                                                                                                                  ccGTTTCTGTACTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATc                           <  1:137412/151‑1 (MQ=255)                                                                                                                            aCTTCTATACCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGc                 <  1:719026/151‑1 (MQ=255)                                                                                                                                  ataCCGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGCGTGGAt           >  2:386731/1‑151 (MQ=255)                                                                                                                                     ccGAGGAAGAATCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGCGTGGATTCg        <  1:754897/151‑1 (MQ=255)                                                                                                                                            agaaTCCCGGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGCGTGGATTCGGCAAAc  >  2:308555/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                              CCGACGATTTGCCCTAATGCGCTGTCGAGAAACTGACTGAAATGGAAAGTCATCGGGTTATCCAGCACGATAATATTTATGGTGCTGGCCAGTGCCCCCATCGAGCCCAGTCGCCGTTTCTGTACTTCTATACCGAGGAAGAATCCCAGCACTGCCAGGCTAATGCACAGCAGCAACATGCTCTGTTGGGTATTAGGGATAATCACCAAAAAGTAGAGCAGCCCTAACGGCAGCGCGGCCAGCGTCCCGTAGATAAAGTCGATCGCCACCATGCGTGGATTCGGCAAAC  >  NC_000913/3386714‑3387002

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