Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,406,008	A→G	D22D (GAT→GAC)	yfbU ←	YfbU family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,406,008	0	A	G	100.0%	122.4 / NA	35	D22D (GAT→GAC)	yfbU	YfbU family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/20);  total (15/20)

CATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATA  >  NZ_CP009273/2405876‑2406137                                                                                                                                     |                                                                                                                                  cATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGa                                                                                                                               >  2:371229/1‑137 (MQ=255) cATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGa                                                                                                                               <  1:371229/137‑1 (MQ=255)   tAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCAt                                                                                                                           <  1:486079/139‑1 (MQ=255)      tGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGt                                                                                                                        >  1:198038/1‑139 (MQ=255)        tCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTTTTTGCCTTCTACGGTTACGTTCAGCATTTGCCGGGTCGATCATAGtct                                                                                                                      <  1:518089/139‑2 (MQ=255)        tCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtca                                                                                                                      >  1:38405/1‑139 (MQ=255)        tCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtca                                                                                                                      <  2:206003/139‑1 (MQ=255)          gatgatAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGtcatc                                                                                                                    <  1:325864/139‑1 (MQ=255)                  tGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTa                                                                                                            >  2:412650/1‑139 (MQ=255)                  tGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTa                                                                                                            >  1:522728/1‑139 (MQ=255)                                                    gcgATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGt                                                                           <  1:278345/138‑1 (MQ=255)                                                      gATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAt                                                                        <  2:198038/139‑1 (MQ=255)                                                             ttCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCatcatc                                                                                                                 <  2:183303/91‑1 (MQ=255)                                                             ttCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCatcatc                                                                                                                 >  1:183303/1‑91 (MQ=255)                                                                aCGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgtt                                                              <  1:224789/139‑1 (MQ=255)                                                                aCGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCAttgtt                                                              <  1:425678/139‑1 (MQ=255)                                                                   cATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGAc                                                           <  2:575190/139‑1 (MQ=255)                                                                   cATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGAc                                                           <  1:412650/139‑1 (MQ=255)                                                                           aTCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTAt                                                   <  1:376101/139‑1 (MQ=255)                                                                           aTCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTAt                                                   <  1:327033/139‑1 (MQ=255)                                                                              cGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAt                                                                        >  1:286115/1‑115 (MQ=255)                                                                              cGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCAt                                                                        <  2:286115/115‑1 (MQ=255)                                                                                tAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGt                                                                           >  1:29977/1‑110 (MQ=255)                                                                                tAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGt                                                                           <  2:29977/110‑1 (MQ=255)                                                                                tAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCActc                                              >  2:149422/1‑139 (MQ=255)                                                                                tAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCActc                                              >  2:323912/1‑139 (MQ=255)                                                                                      cGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGa                                                                                                   >  1:165947/1‑80 (MQ=255)                                                                                      cGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGa                                                                                                   <  2:165947/80‑1 (MQ=255)                                                                                       gCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACg                                       >  2:282508/1‑139 (MQ=255)                                                                                                     tGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAaag                           >  2:281577/1‑137 (MQ=255)                                                                                                     tGCAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAaag                           >  2:142366/1‑137 (MQ=255)                                                                                                       cAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAgg                       <  2:310439/139‑1 (MQ=255)                                                                                                       cAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAgg                       >  1:294749/1‑139 (MQ=255)                                                                                                       cAGGCGACGGTAACGTTCAGCATTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAgg                       <  1:142366/139‑1 (MQ=255)                                                                                                                            aTTTGCCGGGTCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTtata  <  2:522728/139‑1 (MQ=255)                                                                                                                                     |                                                                                                                                  CATAATGTCGATGATAGTGCGGCAGGTTTCTTCTTTCAGCTCGCCAAACTCGCGATCCAATTCACGCATCTGTAATCCGTAACCACGCTCAATAATTGTTTGCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATA  >  NZ_CP009273/2405876‑2406137

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