Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,679,189	A→G	N225S (AAC→AGC)	tus →	DNA replication terminus site‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,679,189	0	A	G	100.0%	140.1 / NA	37	N225S (AAC→AGC)	tus	DNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/16);  total (21/16)

TACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAA  >  NZ_CP009273/1679057‑1679315                                                                                                                                     |                                                                                                                               tACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGcgcg                                                                                                                            >  1:171087/1‑137 (MQ=255)        tgatgaAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTaaa                                                                                                                    >  1:548097/1‑138 (MQ=255)              aGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGTCACAGAGCGCGAAGTTAAAAATCaaa                                                                                                             <  1:477293/139‑1 (MQ=255)                      cacaGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGt                                                                                                     <  1:635003/139‑1 (MQ=255)                                   aaGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGa                                                                                                   <  2:237662/128‑1 (MQ=255)                                   aaGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGa                                                                                                   >  1:237662/1‑128 (MQ=255)                                               accacGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTc                                                                            <  2:171087/139‑1 (MQ=255)                                                         gtCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAg                                                                  >  1:184230/1‑139 (MQ=255)                                                          tCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCGGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAgg                                                                 <  2:548097/139‑1 (MQ=255)                                                             cACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAg                                                                  >  1:296567/1‑135 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGt                                                                                                     <  2:457206/95‑1 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGt                                                                                                     >  1:457206/1‑95 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTc                                                                            >  1:154676/1‑120 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTc                                                                            <  2:154676/120‑1 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCaa                                                          >  1:372314/1‑138 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCaa                                                          >  2:411386/1‑138 (MQ=255)                                                                  tGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCaa                                                          >  2:392669/1‑138 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGcc                                        >  2:13158/1‑139 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGcc                                        >  2:532549/1‑139 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGcc                                        <  1:182098/139‑1 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGcc                                        >  2:266753/1‑139 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGc                                         <  2:601121/138‑1 (MQ=255)                                                                                   gCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGc                                         <  2:184230/138‑1 (MQ=255)                                                                                       aGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTc                                                                            <  1:309947/99‑1 (MQ=255)                                                                                       aGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTc                                                                            >  2:309947/1‑99 (MQ=255)                                                                                         aaaaCTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGccc                                   >  1:366986/1‑138 (MQ=255)                                                                                          aaaCTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTa                                 <  1:392669/139‑1 (MQ=255)                                                                                            aCTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTaca                               <  1:411386/139‑1 (MQ=255)                                                                                                                 tatCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAAt          <  1:13158/139‑1 (MQ=255)                                                                                                                   tCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCcgc                                                                               >  2:624153/1‑68 (MQ=255)                                                                                                                   tCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCcgc                                                                               <  1:624153/68‑1 (MQ=255)                                                                                                                   tCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCg        >  2:284875/1‑139 (MQ=255)                                                                                                                   tCGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCg        >  2:239385/1‑139 (MQ=255)                                                                                                                    cGCTGCCCTGCCACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCgg       >  2:129454/1‑139 (MQ=255)                                                                                                                         ccctgccACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGc                                         >  2:170748/1‑100 (MQ=255)                                                                                                                         ccctgccACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGc                                         <  1:170748/100‑1 (MQ=255)                                                                                                                         ccctgccACAGAGCGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATaa  >  2:326831/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                               TACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAGCCGATTGCCCGCGTCTGGTACAAAGGAGATCAAAAACAAGTCCAACACGCCTGCCCTACACCACTGATTGCACTGATTAATCGGGATAA  >  NZ_CP009273/1679057‑1679315

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