Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,581,105	A→G	V296A (GTC→GCC)	ydeQ ←	fimbrial protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,581,105	0	A	G	100.0%	123.3 / NA	32	V296A (GTC→GCC)	ydeQ	fimbrial protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (20/12);  total (20/12)

TGCGCGGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATT  >  NZ_CP009273/1580977‑1581235                                                                                                                                 |                                                                                                                                   tgCGCGGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGAc                                                                                                                          >  2:279579/1‑139 (MQ=255) tgCGCGGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGAc                                                                                                                          >  2:189217/1‑139 (MQ=255)   cgcgGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACgg                                                                                                                        >  1:848995/1‑139 (MQ=255)   cgcgGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACgg                                                                                                                        >  1:694502/1‑139 (MQ=255)   cgcgGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACgg                                                                                                                        >  1:568497/1‑139 (MQ=255)                aCAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGaa                                                                                                             <  1:189217/137‑1 (MQ=255)                aCAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGaa                                                                                                             <  2:624146/137‑1 (MQ=255)                  aGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCt                                                                                                         >  2:414482/1‑139 (MQ=255)                      gtAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGtt                                                                                                     >  2:816521/1‑139 (MQ=255)                       tAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGttt                                                                                                    >  1:255911/1‑139 (MQ=255)                         ggCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTc                                                                                                   <  2:694502/138‑1 (MQ=255)                              aTAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGttt                                                                                             >  1:560982/1‑139 (MQ=255)                              aTAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGttt                                                                                             >  1:573606/1‑139 (MQ=255)                              aTAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGttt                                                                                             >  1:775689/1‑139 (MQ=255)                                        caAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGt                                                                                   <  1:279579/139‑1 (MQ=255)                                         aaaaTAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGtt                                                                                  <  1:25825/139‑1 (MQ=255)                                               atTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCtt                                                                            <  2:736046/139‑1 (MQ=255)                                                                      tGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGt                                                                                               >  1:227586/1‑97 (MQ=255)                                                                      tGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGt                                                                                               <  2:227586/97‑1 (MQ=255)                                                                      tGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCttt                                                      >  2:68670/1‑138 (MQ=255)                                                                           cgcATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTtttattta                                                >  2:726018/1‑139 (MQ=255)                                                                             cATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACt                                              <  2:525556/139‑1 (MQ=255)                                                                             cATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACt                                              >  1:668552/1‑139 (MQ=255)                                                                                     aaGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAg                                      >  1:370922/1‑139 (MQ=255)                                                                                       ggTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAgg                                     >  2:225733/1‑138 (MQ=255)                                                                                       ggTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAgg                                     >  1:650284/1‑138 (MQ=255)                                                                                            aaaTCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTcc                                                                                                  >  1:307529/1‑72 (MQ=255)                                                                                            aaaTCACAATTATTCATAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTcc                                                                                                  <  2:307529/72‑1 (MQ=255)                                                                                                           aTAGATAAAAGTGACACCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGcc                <  2:370922/139‑1 (MQ=255)                                                                                                                        acacCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCAt   >  1:312726/1‑139 (MQ=255)                                                                                                                         cacCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCAtt  <  2:51591/139‑1 (MQ=255)                                                                                                                         cacCAATGGCTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCAtt  <  1:424432/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                   TGCGCGGATTTTCTTACAGGTGTAGGCTAATAGCCTTCACAAAATAATTATCCGCGCAATTTTCGTGAATGGAGCGCATAGACTAAGGTCAAAATCACAATTATTCATAGATAAAAGTGACACCAATGACTGACTGGACGGTTCCGGCTGAAACCTTGTTTCCTGTTTGTCCATAGGTTGCGCTTAAACCCAGAGGTACTTTGCTTTTATTTACTGTCCCCAGGGAGACATTTTCTCCCGTAGCCAGTATTTTTCCATT  >  NZ_CP009273/1580977‑1581235

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