Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,306,087	A→G	intergenic (+22/+18)	tonB → / ← yciA	TonB system transport protein TonB/acyl‑CoA thioester hydrolase YciA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,306,087	0	A	G	100.0%	71.5 / NA	21	intergenic (+22/+18)	tonB/yciA	TonB system transport protein TonB/acyl‑CoA thioester hydrolase YciA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/14);  total (7/14)

AGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTC  >  NZ_CP009273/1305959‑1306174                                                                                                                                 |                                                                                        aGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCtt                                                                               <  1:248756/139‑1 (MQ=255)   gTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCtttt                                                                             <  2:406057/139‑1 (MQ=255)      aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                          <  2:605199/139‑1 (MQ=255)      aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                          <  1:784252/139‑1 (MQ=255)      aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                          <  2:264481/139‑1 (MQ=255)        aaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTa                                                                        >  1:848172/1‑139 (MQ=255)              gatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTc                                                                  <  1:583546/139‑1 (MQ=255)                tgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCaa                                                                <  1:135414/139‑1 (MQ=255)                               attaGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGa                                                  <  2:177264/135‑1 (MQ=255)                                   gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                             <  1:616980/139‑1 (MQ=255)                                       ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                         >  2:816451/1‑139 (MQ=255)                                       ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                         >  1:588209/1‑139 (MQ=255)                                       ggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCt                                         >  1:630323/1‑139 (MQ=255)                                                       ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                         <  2:848172/139‑1 (MQ=255)                                                       ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                         <  2:786027/139‑1 (MQ=255)                                                       ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                         <  1:145419/139‑1 (MQ=255)                                                                  aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaaa               >  1:608230/1‑138 (MQ=255)                                                                         tGTTTAAAATTAACGGCACCACCGAAATTCAGAAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt       >  1:749727/1‑139 (MQ=255)                                                                          gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                          >  2:366730/1‑71 (MQ=255)                                                                          gTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGAc                                                                          <  1:366730/71‑1 (MQ=255)                                                                               aaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc  <  1:111596/138‑1 (MQ=255)                                                                                                                                 |                                                                                        AGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTC  >  NZ_CP009273/1305959‑1306174

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