Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,349,491	T→C	S470P (TCC→CCC)	gltB →	glutamate synthase large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,349,491	0	T	C	100.0%	181.4 / NA	51	S470P (TCC→CCC)	gltB	glutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (23/28);  total (23/28)

GTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACTCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGTTTGCC  >  NZ_CP009273/3349355‑3349616                                                                                                                                         |                                                                                                                              gTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGAcccc                                                                                                                             >  2:623124/1‑139 (MQ=255) gTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGAcccc                                                                                                                             >  2:541368/1‑139 (MQ=255) gTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGAcccc                                                                                                                             >  1:737744/1‑139 (MQ=255)  tggatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCg                                                                                                                            >  1:220395/1‑139 (MQ=255)  tggatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCg                                                                                                                            >  1:781906/1‑139 (MQ=255)    gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGt                                                                                                                           <  1:167399/138‑1 (MQ=255)    gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGt                                                                                                                           >  2:801031/1‑138 (MQ=255)    gatggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGt                                                                                                                           <  2:737744/138‑1 (MQ=255)     atggaGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTaa                                                                                                                         <  1:623124/139‑1 (MQ=255)              aaCGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTAc                                                                                                                <  2:819733/139‑1 (MQ=255)              aaCGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTAc                                                                                                                >  1:819733/1‑139 (MQ=255)                      ccgACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaa                                                                                                        <  2:225973/139‑1 (MQ=255)                        gACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      <  2:481096/139‑1 (MQ=255)                        gACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      <  1:681780/139‑1 (MQ=255)                        gACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      <  1:32302/139‑1 (MQ=255)                             gTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTc                                                                                                 <  1:716614/139‑1 (MQ=255)                                   tttGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAg                                                                                           >  2:617093/1‑139 (MQ=255)                                   tttGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAg                                                                                           >  1:184554/1‑139 (MQ=255)                                             tGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTc                                                                                 >  1:494789/1‑139 (MQ=255)                                               cccGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGa                                                                               >  2:752148/1‑139 (MQ=255)                                                   aTGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATggg                                                                           <  1:250818/139‑1 (MQ=255)                                                   aTGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATggg                                                                           <  1:649155/139‑1 (MQ=255)                                                   aTGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATggg                                                                           <  1:715189/139‑1 (MQ=255)                                                    tGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGc                                                                          <  1:406797/139‑1 (MQ=255)                                                       agaagTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCgat                                                                       >  2:719060/1‑139 (MQ=255)                                                          agTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCgatgat                                                                    <  2:673407/139‑1 (MQ=255)                                                              ggTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATAccc                                                                >  2:677052/1‑139 (MQ=255)                                                                   ccGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaa                                                                                                        <  1:693839/94‑1 (MQ=255)                                                                   ccGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaa                                                                                                        >  2:693839/1‑94 (MQ=255)                                                                   ccGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTc                                                           >  2:591292/1‑139 (MQ=255)                                                                    cgcgAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCg                                                          <  1:436621/139‑1 (MQ=255)                                                                     gcgAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGc                                                         >  2:314054/1‑139 (MQ=255)                                                                      cgAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGcc                                                        <  1:672363/139‑1 (MQ=255)                                                                        aaCTGGGCGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTCCTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGt                                                      <  2:103174/139‑1 (MQ=255)                                                                         aCTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTg                                                     <  2:510094/139‑1 (MQ=255)                                                                                  cgacACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGt                                            <  2:708359/139‑1 (MQ=255)                                                                                        gCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCcg                                      <  1:314054/139‑1 (MQ=255)                                                                                               aGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCatta                               >  2:424894/1‑139 (MQ=255)                                                                                                      aGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGa                        <  1:2512/139‑1 (MQ=255)                                                                                                              tttAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGcggt                                                                                       <  2:261280/68‑1 (MQ=255)                                                                                                              tttAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGcggt                                                                                       >  1:261280/1‑68 (MQ=255)                                                                                                               ttAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGa                                                                               >  2:784096/1‑75 (MQ=255)                                                                                                               ttAACTACAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGa                                                                               <  1:784096/75‑1 (MQ=255)                                                                                                                      cAGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCcagcag        >  2:565400/1‑139 (MQ=255)                                                                                                                       aGCGCGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGt       <  2:155221/139‑1 (MQ=255)                                                                                                                        gcgcGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      <  2:619204/43‑1 (MQ=255)                                                                                                                        gcgcGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      >  1:619204/1‑43 (MQ=255)                                                                                                                        gcgcGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      >  2:762957/1‑43 (MQ=255)                                                                                                                        gcgcGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGaaaa                                                                                                      <  1:762957/43‑1 (MQ=255)                                                                                                                        gcgcGGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGtt      >  1:159354/1‑139 (MQ=255)                                                                                                                         cgcgGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGttt     >  2:734418/1‑139 (MQ=255)                                                                                                                           cgGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGTTTGc   >  2:137676/1‑139 (MQ=255)                                                                                                                           cgGAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGGTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGTTTGc   >  2:51257/1‑139 (MQ=255)                                                                                                                            ggAAGAGCTGGACCCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGTTTGcc  <  1:591292/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                              GTGGATGGAGAAAAACGTCCGCCGACTGGTACCGTTTGAAGATCTGCCCGATGAAGAAGTGGGTAGCCGCGAACTGGACGACGACACGCTTGCCAGCTACCAGAAACAGTTTAACTACAGCGCGGAAGAGCTGGACTCCGTAATTCGCGTACTGGGCGAAAACGGTCAGGAAGCGGTCGGTTCGATGGGCGATGATACCCCATTCGCCGTGCTCTCCAGTCAGCCGCGCATTATTTACGACTACTTCCGCCAGCAGTTTGCC  >  NZ_CP009273/3349355‑3349616

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