Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,873,583	A→G	S333S (AGT→AGC)	recF ←	DNA replication/repair protein RecF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,873,583	0	A	G	100.0%	129.8 / NA	37	S333S (AGT→AGC)	recF	DNA replication/repair protein RecF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/21);  total (16/21)

TTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGT  >  NZ_CP009273/3873447‑3873714                                                                                                                                         |                                                                                                                                     tttGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCt                                                                                                                                    >  2:173441/1‑139 (MQ=255)   tGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGa                                                                                                                                  >  1:113187/1‑139 (MQ=255)   tGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGa                                                                                                                                  >  1:107356/1‑139 (MQ=255)   tGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGa                                                                                                                                  >  1:76905/1‑139 (MQ=255)        cTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATcgcg                                                                                                                             <  1:173441/139‑1 (MQ=255)        cTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATcgcg                                                                                                                             <  2:113187/139‑1 (MQ=255)        cTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATcgcg                                                                                                                             <  1:425718/139‑1 (MQ=255)        cTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATcgcg                                                                                                                             <  1:96752/139‑1 (MQ=255)           gaggagTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCgcgc                                                                                                                            <  2:354204/137‑1 (MQ=255)               agTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACaa                                                                                                                      <  2:316258/139‑1 (MQ=255)                           tCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGCTCAGCGCTGATCGCGCTGACAAAGACCTGTGAtt                                                                                                          >  1:324444/1‑139 (MQ=255)                               cATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACAGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGt                                                                                                      <  2:58987/139‑1 (MQ=255)                                                           gggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCCGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTc                                                                                                     >  1:437324/1‑112 (MQ=255)                                                           gggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTc                                                                                                     <  2:437324/112‑1 (MQ=255)                                                           gggTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑Ccg                                                                         >  2:154572/1‑139 (MQ=255)                                                                    cGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCtca                                                                <  2:147263/139‑1 (MQ=255)                                                                      ttATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCtcat                                                               >  2:78371/1‑138 (MQ=255)                                                                                cTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCa                                                    <  1:201456/139‑1 (MQ=255)                                                                                                     cGAATTTTCGTCCGTCATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                   <  2:28410/135‑1 (MQ=255)                                                                                                     cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                   <  2:76905/135‑1 (MQ=255)                                                                                                     cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTAt                                   <  2:45110/135‑1 (MQ=255)                                                                                                     cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAgg                               >  2:75702/1‑139 (MQ=255)                                                                                                     cGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAgg                               >  1:214081/1‑139 (MQ=255)                                                                                                                      tGTCTATAACGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGccgccc              <  2:242356/139‑1 (MQ=255)                                                                                                                        tCTGTAACGTGTTCAGCGCTTATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCGCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGccgccc              <  2:324444/138‑1 (MQ=255)                                                                                                                              aCGTGTTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCAc      <  2:214081/139‑1 (MQ=255)                                                                                                                                 tgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAatc                                         >  1:400635/1‑101 (MQ=255)                                                                                                                                 tgtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAatc                                         <  2:400635/101‑1 (MQ=255)                                                                                                                                  gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCtcatcatc                                                           >  2:224361/1‑82 (MQ=255)                                                                                                                                  gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCtcatcatc                                                           <  1:224361/82‑1 (MQ=255)                                                                                                                                  gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCa       >  2:312912/1‑134 (MQ=255)                                                                                                                                  gtTCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCa       <  1:312912/134‑1 (MQ=255)                                                                                                                                   ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCt                                     <  2:292885/103‑1 (MQ=255)                                                                                                                                   ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCt                                     >  1:292885/1‑103 (MQ=255)                                                                                                                                   ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCAc      >  1:242994/1‑134 (MQ=255)                                                                                                                                   ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCAc      <  2:242994/134‑1 (MQ=255)                                                                                                                                   ttCAGCGCTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGt  >  2:143804/1‑138 (MQ=255)                                                                                                                                         |                                                                                                                                     TTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGC‑CCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGT  >  NZ_CP009273/3873447‑3873714

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