Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,360,769	G→A	intergenic (+153/‑407)	gltD → / → gltF	glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,360,769	0	G	A	100.0%	116.2 / NA	37	intergenic (+153/‑407)	gltD/gltF	glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/18);  total (19/18)

TCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAA  >  NC_000913/3360656‑3360892                                                                                                                  |                                                                                                                            tCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTtctc                                                                                                    >  1:438106/1‑139 (MQ=255)                    cgcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCGATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAg                                                                                                             <  1:1491614/111‑1 (MQ=255)                    cgcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCGATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAg                                                                                                             >  2:1491614/1‑111 (MQ=255)                      cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                                               >  1:219811/1‑139 (MQ=255)                            tttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACaa                                                                                            >  1:1155450/1‑120 (MQ=255)                            tttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACaa                                                                                            <  2:1155450/120‑1 (MQ=255)                             ttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                        <  2:219811/139‑1 (MQ=255)                              tCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                       <  2:438106/139‑1 (MQ=255)                                     aaTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAATCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                       >  2:1012070/1‑132 (MQ=255)                                     aaTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                       <  1:1012070/132‑1 (MQ=255)                                            cATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATa                                                         >  2:373180/1‑139 (MQ=255)                                                    tttGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTtct                                                                                                     <  1:1195332/87‑1 (MQ=255)                                                    tttGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTtct                                                                                                     >  2:1195332/1‑87 (MQ=255)                                                                     taGGGATAATTTTCAGTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAAttt                                                                                                       <  1:1017602/68‑1 (MQ=255)                                                                     taGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAAttt                                                                                                       >  2:1017602/1‑68 (MQ=255)                                                                       gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACa                                                                                             >  1:1325756/1‑76 (MQ=255)                                                                       gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACa                                                                                             <  2:1325756/76‑1 (MQ=255)                                                                        ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                        <  2:840355/96‑1 (MQ=255)                                                                        ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                        >  1:840355/1‑96 (MQ=255)                                                                        ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                     <  2:346511/115‑1 (MQ=255)                                                                        ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                     >  1:346511/1‑115 (MQ=255)                                                                        ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                             >  2:538355/1‑139 (MQ=255)                                                                          aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                    >  1:519842/1‑130 (MQ=255)                                                                          aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                    <  2:519842/130‑1 (MQ=255)                                                                          aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCAt                           <  1:538355/139‑1 (MQ=255)                                                                                  cAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTa                   >  2:467594/1‑139 (MQ=255)                                                                                        aTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGt             <  1:373180/139‑1 (MQ=255)                                                                                          tgtgCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGtct           <  1:467594/139‑1 (MQ=255)                                                                                                 tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                     <  1:1116705/90‑1 (MQ=255)                                                                                                 tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                     >  2:1116705/1‑90 (MQ=255)                                                                                                   ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTaa  >  2:1409636/1‑139 (MQ=255)                                                                                                      aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAAt                                                                                           <  1:394641/47‑1 (MQ=255)                                                                                                      aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAAt                                                                                           >  2:394641/1‑47 (MQ=255)                                                                                                        cAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                      >  2:1354149/1‑114 (MQ=255)                                                                                                        cAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                      <  1:1354149/114‑1 (MQ=255)                                                                                                           gTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAAt                                                                                 <  2:698449/52‑1 (MQ=255)                                                                                                           gTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAAt                                                                                 >  1:698449/1‑52 (MQ=255)                                                                                                                  |                                                                                                                            TCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAA  >  NC_000913/3360656‑3360892

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