Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,515,380	T→C	Y648C (TAC→TGC)	yncD ←	TonB‑dependent receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,515,380	0	T	C	100.0%	169.1 / NA	51	Y648C (TAC→TGC)	yncD	TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (12/39);  total (12/39)

TTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCC  >  NZ_CP009273/1515293‑1515465                                                                                        |                                                                                      tttGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCaa                                                                                     <  2:160103/90‑1 (MQ=255)  ttGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAAt                                                                                    <  1:1135235/90‑1 (MQ=255)   tGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAAtt                                                                                   <  2:1537794/90‑1 (MQ=255)   tGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAAtt                                                                                   <  2:491918/90‑1 (MQ=255)    gACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTg                                                                                  <  2:1298444/90‑1 (MQ=255)    gACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTg                                                                                  <  2:749676/90‑1 (MQ=255)       tCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCAGAGTTAAATTGTGGCAATTGTAt                                                                               >  2:4162/1‑90 (MQ=255)          ttGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTta                                                                            >  1:1261081/1‑90 (MQ=255)               aaTGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATAccc                                                                       <  1:1282128/90‑1 (MQ=255)                 tGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCgg                                                                     <  2:1114866/90‑1 (MQ=255)                 tGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCgg                                                                     <  1:1543732/90‑1 (MQ=255)                      gAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAAt                                                                >  1:1271184/1‑90 (MQ=255)                      gAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAAt                                                                >  2:1180231/1‑90 (MQ=255)                      gAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAAt                                                                >  2:725541/1‑90 (MQ=255)                       aaCCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATa                                                               <  2:1418675/90‑1 (MQ=255)                           aaCGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGcc                                                           >  1:1207985/1‑90 (MQ=255)                           aaCGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGcc                                                           >  1:1120301/1‑90 (MQ=255)                           aaCGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGcc                                                           >  1:1111718/1‑90 (MQ=255)                               tATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCcgacg                                                       <  2:336780/90‑1 (MQ=255)                               tATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCcgacg                                                       <  1:1565283/90‑1 (MQ=255)                               tATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCcgacg                                                       <  1:189917/90‑1 (MQ=255)                                aTTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCcgacga                                                      >  1:364225/1‑90 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:1290793/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:98153/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:1366549/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:493577/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:670445/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  1:1001331/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  2:1185994/90‑1 (MQ=255)                                 ttCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACgag                                                     <  1:1180231/90‑1 (MQ=255)                                      tATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGtat                                                >  2:1237522/1‑90 (MQ=255)                                                   aTCGACACGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCttt                                   <  1:1031054/90‑1 (MQ=255)                                                        cacGACCAAATAAATCCAGAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCg                              <  1:4162/90‑1 (MQ=255)                                                        cacGACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCg                              <  1:1344482/90‑1 (MQ=255)                                                           gACCAAATAAATCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTAt                           >  1:1470567/1‑90 (MQ=255)                                                                 ataaatCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTtcat                     <  1:1237522/90‑1 (MQ=255)                                                                    aatCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTcatcat                  >  1:846761/1‑90 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  2:146891/90‑1 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  2:1510021/90‑1 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  2:529439/90‑1 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  2:260133/90‑1 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  2:342743/90‑1 (MQ=255)                                                                      tCCACAGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCt                <  1:1366816/90‑1 (MQ=255)                                                                           aGTTAAATTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCat           <  1:1507641/90‑1 (MQ=255)                                                                               aaaTTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCataata       <  2:969262/90‑1 (MQ=255)                                                                                 aTTGTGGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATc     <  1:1198447/90‑1 (MQ=255)                                                                                   tgtgGCATTTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGc   <  2:48089/90‑1 (MQ=255)                                                                                   tgtgGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGc   <  1:254377/90‑1 (MQ=255)                                                                                   tgtgGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGc   <  1:960765/90‑1 (MQ=255)                                                                                    gtgGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGcc  <  1:499923/90‑1 (MQ=255)                                                                                    gtgGCAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGcc  <  2:1261081/90‑1 (MQ=255)                                                                                        |                                                                                      TTTGACTCATTGACAATGACAGAACCAACGTATTCTTTATCGAATAAATTATCGACACGACCAAATAAATCCACAGTTAAATTGTGGTAATTGTATTTATACCCGGTGAATAAGCCGACGAGAGTATAAGACGGCGCTTTTGCCGTATTTTCATCATCTGCCATAATATCGCC  >  NZ_CP009273/1515293‑1515465

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