Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,625,832	T→C	Q21Q (CAA→CAG)	guaA ←	glutamine‑hydrolyzing GMP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,625,832	0	T	C	100.0%	135.6 / NA	35	Q21Q (CAA→CAG)	guaA	glutamine‑hydrolyzing GMP synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/17);  total (18/17)

CAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAG  >  NZ_CP009273/2625702‑2625966                                                                                                                                   |                                                                                                                                       cAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAgt                                                                                                                                   <  1:123200/136‑1 (MQ=255)  agtagtACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTa                                                                                                                                <  2:483651/138‑1 (MQ=255)      gtaCTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTgag                                                                                                                           <  2:606738/139‑1 (MQ=255)            tCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                      >  1:460479/1‑138 (MQ=255)              cGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAg                                                                                                                   >  2:244257/1‑139 (MQ=255)                                 atGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgat                                                                                                <  1:847993/139‑1 (MQ=255)                                  tGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcgatg                                                                                               >  1:660422/1‑139 (MQ=255)                                       cTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                      >  1:118654/1‑111 (MQ=255)                                       cTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCg                                                                                                                      <  2:118654/111‑1 (MQ=255)                                          ggATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATg                                                                                         <  2:660422/137‑1 (MQ=255)                                                  gTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAGCCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCAATGCTTATGAATGTTTTcc                                                                               <  1:398111/139‑1 (MQ=255)                                                     aCGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                  <  1:135889/117‑1 (MQ=255)                                                     aCGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGgatgcg                                                                                                  >  2:135889/1‑117 (MQ=255)                                                     aCGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGt                                                                             >  2:519203/1‑138 (MQ=255)                                                              tgCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCt                                                                   <  2:804216/139‑1 (MQ=255)                                                                             ccACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGt                                                    >  1:604733/1‑139 (MQ=255)                                                                               aCGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTg                                                   >  1:127653/1‑138 (MQ=255)                                                                                cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGa                                                  >  1:704015/1‑138 (MQ=255)                                                                                cGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGa                                                  >  1:295992/1‑138 (MQ=255)                                                                                    cacaGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCaa                                                      >  1:252482/1‑130 (MQ=255)                                                                                    cacaGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCaa                                                      <  2:252482/130‑1 (MQ=255)                                                                                       aGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGaa                                                 >  1:468552/1‑132 (MQ=255)                                                                                       aGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGaa                                                 <  2:468552/132‑1 (MQ=255)                                                                                           cGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGataata                                       >  1:654825/1‑138 (MQ=255)                                                                                                   acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                              <  1:573376/139‑1 (MQ=255)                                                                                                   acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                              >  2:573376/1‑139 (MQ=255)                                                                                                   acacCCAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGcc                              <  1:383599/139‑1 (MQ=255)                                                                                                        cAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACa                         >  1:733630/1‑139 (MQ=255)                                                                                                        cAGCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACa                         <  1:519203/139‑1 (MQ=255)                                                                                                          gCTCACGCACGCGGCGCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                        >  2:372102/1‑138 (MQ=255)                                                                                                                     cggcgCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg             >  2:399131/1‑138 (MQ=255)                                                                                                                     cggcgCGCAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGa            <  1:399131/139‑1 (MQ=255)                                                                                                                        cgcgcAACCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa           >  2:457354/1‑137 (MQ=255)                                                                                                                             aaCCAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATc    <  2:152605/139‑1 (MQ=255)                                                                                                                               ccAGCTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAg  <  1:661837/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       CAGTAGTACTTTCCGGGCCGCCGGAAAGAATAATGCCGCTTGGATTGAAGTCACGAATTTGTGCTTCTGTCACATCCCACGCCCACAGTTCGCAGTAAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAG  >  NZ_CP009273/2625702‑2625966

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