Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,632,046	G→T	P164T (CCG→ACG)	guaA ←	GMP synthetase (glutamine aminotransferase)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,632,046	0	G	T	100.0%	113.4 / NA	33	P164T (CCG→ACG)	guaA	GMP synthetase (glutamine aminotransferase)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (15/18);  total (15/18)

CTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTG  >  NC_000913/2631905‑2632177                                                                                                                                              |                                                                                                                                    cTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCtt                                                                                                                              >  1:126969/1‑149 (MQ=255)    gCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCg                                                                                                                           <  8:129578/149‑1 (MQ=255)            aGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGc                                                                                                                   >  1:50569/1‑149 (MQ=255)                    acacTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtga                                                                                                           <  5:77388/149‑1 (MQ=255)                      acTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGtgatg                                                                                                         <  2:126969/149‑1 (MQ=255)                                   aCGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGc                                                                                                                   <  5:10237/126‑1 (MQ=255)                                   aCGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGc                                                                                                                   >  6:10237/1‑126 (MQ=255)                                   aCGCACAAAACGCTCCAGCATGCGCATACCCAGCCGGGTATGAGACACTACCGGGAGCAACACTACGCCATAGAAGCCTTTTTCTTCGATAGCCATAATGCCAAACGTGCAGATTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGaa                                                                                            >  4:6427/1‑149 (MQ=255)                                         aaaaCGCTCCAGCATGCGCATACCCTGCCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCgg                                                                                      >  1:8286/1‑149 (MQ=255)                                          aaaCGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGt                                                                                     <  2:50569/149‑1 (MQ=255)                                                   aGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATc                                                                            >  8:64813/1‑149 (MQ=255)                                                                  tGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                     >  1:119600/1‑141 (MQ=255)                                                                  tGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTgg                                                                     <  2:119600/141‑1 (MQ=255)                                                                    gCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGa                                                           <  1:97305/149‑1 (MQ=255)                                                                        gTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGGTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATc                                                       <  7:64813/149‑1 (MQ=255)                                                                         tATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATc                                                       >  1:174216/1‑148 (MQ=255)                                                                            gAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTa                                                                                                                 <  6:141982/87‑1 (MQ=255)                                                                            gAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTa                                                                                                                 >  5:141982/1‑87 (MQ=255)                                                                                      cGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCg                                               <  3:33698/143‑1 (MQ=255)                                                                                      cGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCg                                               >  4:33698/1‑143 (MQ=255)                                                                                            gAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTg                                   >  8:62212/1‑149 (MQ=255)                                                                                            gAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTg                                   <  2:8286/149‑1 (MQ=255)                                                                                                tGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGaa                                                                                                       >  3:85216/1‑77 (MQ=255)                                                                                                tGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGaa                                                                                                       <  4:85216/77‑1 (MQ=255)                                                                                                    cGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAg                            <  2:174216/148‑1 (MQ=255)                                                                                                           gAAGCGTTTTTCTTCGTTAGCCTTAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGGCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCGAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCtt                    <  3:6427/149‑1 (MQ=255)                                                                                                           gAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCTGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTTGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCtt                    >  4:101168/1‑149 (MQ=255)                                                                                                                gTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTc                                     <  3:1367/127‑1 (MQ=255)                                                                                                                gTTTTTCTTCGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCAGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTc                                     >  4:1367/1‑127 (MQ=255)                                                                                                                     tcttcGTTAGCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACcgc           <  6:162754/148‑1 (MQ=255)                                                                                                                       cacGTTTGCCATTATGGCTAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTa                                          <  2:120759/113‑1 (MQ=255)                                                                                                                       cacGTTTGCCATTATGGCTAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTa                                          >  1:120759/3‑115 (MQ=255)                                                                                                                              gCCATAATGGCAAACGTGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTg  <  7:62212/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                    CTGGCGTCCACAGGGCTTCACACTGGCAGATATCACGCACAAAACGCTCCAGCATGCGCATACCCTGGCGGGTATGAGTCACTTCCGGGTGGAACTGTACGCCATAGAAGCGTTTTTCTTCGTTAGCCATAATGGCAAACGGGCAGCTTTCGGTGCTGGCTACGGTGATGAAGTCGGACGGAATAGCGGTAACTTTATCGCCGTGGCTCATCCAGACATCGAGCAGCGGTTTACCGTCTGCGGTCAGCGCATCTTCGATACCGCGAACCAGTG  >  NC_000913/2631905‑2632177

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