Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	861,316	T→C	I169V (ATT→GTT)	moeA ←	molybdopterin molybdotransferase MoeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	861,316	0	T	C	100.0%	64.8 / NA	19	I169V (ATT→GTT)	moeA	molybdopterin molybdotransferase MoeA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/12);  total (7/12)

CTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAATCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCA  >  NZ_CP009273/861181‑861450                                                                                                                                        |                                                                                                                                       cTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAAcccc                                                                                                                                     >  2:346391/1‑139 (MQ=255)         aGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAg                                                                                                                             <  2:276013/139‑1 (MQ=255)              tACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAAt                                                                                                                         >  2:390210/1‑138 (MQ=255)                     agaCGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTgg                                                                                                                   >  1:390682/1‑137 (MQ=255)                              gTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCg                                                                                                         <  1:346391/138‑1 (MQ=255)                                 tCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTa                                                                                                      <  1:227823/138‑1 (MQ=255)                                                              aCCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGATTTCCCGCCGGGAAaacaac                                                                         <  1:425324/138‑1 (MQ=255)                                                                 gggCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCa                                                                     <  2:358434/139‑1 (MQ=255)                                                                       cTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCa                                                               <  1:105452/139‑1 (MQ=255)                                                                          gAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCa                                                               >  2:431711/1‑136 (MQ=255)                                                                              tcatcaCCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCAc                                                                    >  1:324942/1‑127 (MQ=255)                                                                              tcatcaCCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCAc                                                                    <  2:324942/127‑1 (MQ=255)                                                                                 tcaCCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCtt                                                     >  1:145603/1‑139 (MQ=255)                                                                                     cGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAg                                                                                                                             <  1:28898/63‑1 (MQ=255)                                                                                     cGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAACCCCCAGTGAAg                                                                                                                             >  2:28898/1‑63 (MQ=255)                                                                                                               ttACGAATCACCGGAACTTCGGCAACCCCCAGTGAAg                                                                                                                             >  2:58072/1‑37 (MQ=255)                                                                                                               ttACGAATCACCGGAACTTCGGCAACCCCCAGTGAAg                                                                                                                             <  1:58072/37‑1 (MQ=255)                                                                                                                 aCGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACg                      <  2:390682/138‑1 (MQ=255)                                                                                                                 aCGAATCACCGGAACTTCGGCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACg                      <  1:431711/138‑1 (MQ=255)                                                                                                                                    gCAACCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAacgca  <  2:145603/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                       CTAACATCAGGTGTACGGCGAGACGGTTGGTATCGTAGATTTGGCCGTCGCCCAGCGGCTGACCGGGCAACTGGAGTTCATCACCGGTAGAAAAAAGCGCTACACGCACTTTACGAATCACCGGAACTTCGGCAATCCCCAGTGAAGCAATTACTGGCAGCTCTGCGGTAGTCAGGCGAGTTCCCGCCGGGAAAACAACCGCACCTGCAGAGATATCTTCACCGCGACGGCGAATATTTTGCCCGCTACGCACTTCAGCAGTAAAACGCA  >  NZ_CP009273/861181‑861450

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