Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,521,249	A→C	K315N (AAA→AAC)	yfeH →	bile acid:sodium symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,521,249	0	A	C	100.0%	114.9 / NA	32	K315N (AAA→AAC)	yfeH	bile acid:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (17/15);  total (17/15)

AAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCA  >  NZ_CP009273/2521113‑2521378                                                                                                                                         |                                                                                                                                  aaaaGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCg                                                                                                                                 <  2:92213/139‑1 (MQ=255)      agTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGa                                                                                                                            <  2:86839/139‑1 (MQ=255)            gCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGaa                                                                                                                       <  1:187973/138‑1 (MQ=255)                        ccGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGgcgc                                                                                                          <  2:427898/139‑1 (MQ=255)                                             cccACATCGGTTATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGAt                                                                                      <  1:486327/138‑1 (MQ=255)                                                            ggTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  2:439611/138‑1 (MQ=255)                                                              tATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGc                                                                    >  1:479486/1‑139 (MQ=255)                                                                      tgctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAgggg                                                            >  2:126968/1‑139 (MQ=255)                                                                       gctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  1:318854/127‑1 (MQ=255)                                                                       gctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       >  2:318854/1‑127 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCg                                                                                                                                 <  2:639923/68‑1 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCg                                                                                                                                 >  1:639923/1‑68 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  2:741483/126‑1 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       >  1:741483/1‑126 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCt                                                          >  1:796734/1‑139 (MQ=255)                                                                        ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCt                                                          >  2:339139/1‑139 (MQ=255)                                                                               tgatgaTTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAg                                                   >  1:554597/1‑139 (MQ=255)                                                                                gatgatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGt                                                  >  2:307447/1‑139 (MQ=255)                                                                                 atgatTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAg                                                   <  1:339139/137‑1 (MQ=255)                                                                                                aTCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCgg                                  <  1:685198/139‑1 (MQ=255)                                                                                                          tGGTCTGTGCGGTGCTGGCGCGGCGATACAGCCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATcag                        <  1:397094/139‑1 (MQ=255)                                                                                                           ggTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  1:9091/91‑1 (MQ=255)                                                                                                           ggTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       >  2:9091/1‑91 (MQ=255)                                                                                                                      tGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTg            >  2:33303/1‑139 (MQ=255)                                                                                                                        cTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCa          >  1:770049/1‑139 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             >  1:694300/1‑37 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             <  2:694300/37‑1 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt         >  1:193743/1‑137 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt         >  2:730553/1‑137 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt         >  2:575127/1‑137 (MQ=255)                                                                                                                           gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCACCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt         >  2:719422/1‑137 (MQ=255)                                                                                                                                tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  <  2:193743/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                  AAAAGAGTCTGGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCA  >  NZ_CP009273/2521113‑2521378

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