Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_vch_tpiA	482,230	G→A	P725L (CCG→CTG)	acrB ←	multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_vch_tpiA	482,230	0	G	A	100.0%	116.4 / NA	39	P725L (CCG→CTG)	acrB	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/18);  total (21/18)

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCC  >  NC_000913_3_vch_tpiA/482084‑482369                                                                                                                                                   |                                                                                                                                            tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                           <  2:420133/149‑1 (MQ=255)  gacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCcggg                                                                                                                                          >  1:738831/1‑149 (MQ=255)       aaaCTTTCTTCACACGACCGCGGTCGCTAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATCTCGTTGATAGAAACACCAAGCGCCTGCGCTTTATCCTGGTCGATATCAATCTTAAACTGCAGGGTATCt                                                                                                                                     >  1:553921/1‑149 (MQ=255)           tttCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCa                                                                                                                                 <  2:66087/149‑1 (MQ=255)                  acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                          >  1:15262/1‑149 (MQ=255)                          gcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTa                                                                                                                  <  1:248527/149‑1 (MQ=255)                               cGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTg                                                                                                             <  1:273317/149‑1 (MQ=255)                                 atattggcGTGCACATTGCTGCCGCCCCATGCAGCGTCCAGATTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                           <  2:193621/142‑1 (MQ=255)                                  tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                          <  1:291581/149‑1 (MQ=255)                                     aGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATAGCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGCCGTACGCTGGTCAAc                                                                                                       >  1:658507/1‑149 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGc                                                                                                               >  2:763653/1‑137 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCCACatat                                                                                                   >  1:430825/1‑149 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                   >  1:492082/1‑149 (MQ=255)                                         gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTGGGACGTACGc                                                                                                               <  1:763653/137‑1 (MQ=255)                                          ttCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATc                                                                                                  >  1:666968/1‑149 (MQ=255)                                          tacacaTACCGGCCGCCCCATGCAGCGCCCAGTGTGGTGTTAATATCGTTGATAGAAACACCCAGCGCCTGCGGTTTGTCCTGGTCTTTATCACTCTTAAACTGCAGGGGATCTTCCAGATCGTTTGGACGTACGCTGGTCAACATATc                                                                                                  <  1:268073/147‑1 (MQ=255)                                           tCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATCTCACTCTTAAACTGCAGGGTATATCACCTATCGTCTGGACCTACGCTGGTCAACATATCa                                                                                                 <  2:581972/149‑1 (MQ=255)                                             acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTAGATATCAATCTTAAACTGCAGGGTATATTACAGACCGTTTGGACGTACTCTCGTCCACATATGAgc                                                                                               >  2:37183/1‑148 (MQ=255)                                                     gccgccCCATCCACCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCg                                                                                       <  1:392266/149‑1 (MQ=255)                                                       cgccCCAGGCACGGCCCTGTGTGGCTTTAAGGTCCTTGATAGAACCACCCAGCGCCTGCGCGTTTTGCTGGTCGAAATCAATCGGAAACTGGAGGGTATATGCCAGACCGTTTGGACGTACGCTGGTCAACATAGCAGGGTGCTTCgct                                                                                     <  2:669335/149‑1 (MQ=255)                                                           ccATGCCGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATTCCAGGGTGCTTCGCTGCtt                                                                                 >  1:49629/1‑149 (MQ=255)                                                              tGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTg                                                                              >  1:29496/1‑149 (MQ=255)                                                                  gcgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCCGGGTGCTTCGCTGCTCCTgcaag                                                                          >  2:326447/1‑149 (MQ=255)                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                             <  2:492082/149‑1 (MQ=255)                                                                               tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                             <  2:18007/149‑1 (MQ=255)                                                                                   aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCt                                                         <  2:15262/149‑1 (MQ=255)                                                                                          ttGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAg                                                  >  1:402152/1‑149 (MQ=255)                                                                                           tGATAGAAACACCCGGCGCCTGCGCTTTTTCCTCGTTGATATCAATCTTAAACTGCAGGGTATCTTCCACACCGTTTGGACGTACGCTGGTCAACATGTCAGGGTGCTTCGCTGCTTCTGCACGCACCGGGTTCCGCGCCCGAGTCacg                                                 >  2:169341/1‑147 (MQ=255)                                                                                              tAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCGGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCGGGGTCAGtttt                                              >  2:649778/1‑149 (MQ=255)                                                                                               aGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACCGCAGGAGAGCTGCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTAACCCGCCTGAGTCAGttttt                                             >  1:152628/1‑149 (MQ=255)                                                                                                gAAACACCCAGCTCCTGAGCTTTTTTCTGTTTTTTATCAATCTAAAACTGCAGGGTATCTTCCAGACTGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCACGCAACTGGTTACGCCCCTGAGTCAGTTTTTc                                            >  2:565908/1‑149 (MQ=255)                                                                                                     acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCCGTATGTCTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGAc                                       >  1:542966/1‑149 (MQ=255)                                                                                                       ccAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCa                                     >  1:5245/1‑149 (MQ=255)                                                                                                            gcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGcc                                >  1:760399/1‑149 (MQ=255)                                                                                                                    tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                        <  2:760399/149‑1 (MQ=255)                                                                                                                        tCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                    <  1:738116/149‑1 (MQ=255)                                                                                                                                 atatCAATCTTAAACTGCATGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACCGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa           <  1:123878/149‑1 (MQ=255)                                                                                                                                 atatCAATCTTAAACGGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa           <  2:658507/149‑1 (MQ=255)                                                                                                                                          ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                              <  2:471637/105‑1 (MQ=255)                                                                                                                                          ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtttt                                              >  1:471637/1‑105 (MQ=255)                                                                                                                                          ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc  >  1:25045/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                            TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCC  >  NC_000913_3_vch_tpiA/482084‑482369

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑