Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,315,373	T→C	intergenic (‑401/+58)	yhbX ← / ← BW25113_RS16455	phosphoethanolamine transferase/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,315,373	0	T	C	100.0%	110.3 / NA	30	intergenic (‑401/+58)	yhbX/BW25113_RS16455	phosphoethanolamine transferase/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/12);  total (18/12)

AAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTA  >  NZ_CP009273/3315294‑3315459                                                                                |                                                                                       aaGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAg                                                                             >  1:6538/1‑91 (MQ=255)      gTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTc                                                                                    <  1:311123/79‑1 (MQ=255)      gTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTc                                                                                    >  2:311123/1‑79 (MQ=255)            gCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGa                                                                  <  1:697494/91‑1 (MQ=255)             cACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGCCCAGGACCAAAACGaa                                                                 >  1:366185/1‑91 (MQ=255)             cACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGaa                                                                 >  1:661305/1‑91 (MQ=255)             cACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGaa                                                                 >  2:429142/1‑91 (MQ=255)                 aTCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGaaaaaa                                                             >  1:439666/1‑91 (MQ=255)                    gTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGAc                                                          <  1:51988/91‑1 (MQ=255)                            gTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCa                                                  >  2:521422/1‑91 (MQ=255)                                 aGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                             <  2:763363/91‑1 (MQ=255)                                 aGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtc                                             <  1:429142/91‑1 (MQ=255)                                   cTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGtctc                                           <  1:782145/91‑1 (MQ=255)                                      cGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCttt                                        >  1:794799/1‑91 (MQ=255)                                      cGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCttt                                        >  2:538119/1‑91 (MQ=255)                                             aGCAACTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGaa                                 <  1:658664/91‑1 (MQ=255)                                                 aCTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTg                             >  2:89895/1‑91 (MQ=255)                                                  cTGGACGAGATACAGATACCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTgg                            <  1:538119/91‑1 (MQ=255)                                                         agatacagatacCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAg                     >  1:127620/1‑91 (MQ=255)                                                              cagatacCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                <  1:775269/91‑1 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                >  1:432874/1‑86 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACgg                <  2:432874/86‑1 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt           >  2:513348/1‑91 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt           >  1:343538/1‑91 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt           >  2:805065/1‑91 (MQ=255)                                                                   acCTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACtt           >  1:101929/1‑91 (MQ=255)                                                                     cTGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGa         >  1:346055/1‑91 (MQ=255)                                                                      tGACAACCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGaa        <  1:140912/91‑1 (MQ=255)                                                                           aCCATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGt   >  2:565051/1‑91 (MQ=255)                                                                            ccATCCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTa  <  2:127620/91‑1 (MQ=255)                                                                                |                                                                                       AAGGCGTACCGGCACCATCGTTTCAAGGTACCAGCTACGAGTAAAGCAACTGGACGAGATACAGATACCTGACAACCATTCCTCAGACCAGGACCAAAACGAAAAAAGACGCTTTTCAGCGTCTCTTTTCTGGAATTTGGTACCGAGGACGGGACTTGAACCCGTA  >  NZ_CP009273/3315294‑3315459

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