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%	130.4 / NA	41	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 (17/24);  total (17/24)

GCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATA  >  NC_000913/3360651‑3360896                                                                                                                       |                                                                                                                                gCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAg                                                                                                                               <  2:852666/121‑1 (MQ=255)                gACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATc                                                                                              <  2:1052349/139‑1 (MQ=255)                        gcgcGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGt                                                                                           <  2:642401/134‑1 (MQ=255)                        gcgcGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGt                                                                                           >  1:642401/1‑134 (MQ=255)                                 tttCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGtt                                                                             >  2:243417/1‑139 (MQ=255)                                       gCAAATGCGGCATGATCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                       <  1:243417/139‑1 (MQ=255)                                       gCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                       <  1:864743/139‑1 (MQ=255)                                                    gAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCg                                                          <  2:33965/139‑1 (MQ=255)                                                     aGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGa                                                         >  2:187424/1‑139 (MQ=255)                                                         tttGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCt                                                     <  2:667330/139‑1 (MQ=255)                                                                    cTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTc                                          <  1:41451/139‑1 (MQ=255)                                                                      ataataGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        <  1:187424/139‑1 (MQ=255)                                                                            gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                            >  1:750411/1‑97 (MQ=255)                                                                            gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                            <  2:750411/97‑1 (MQ=255)                                                                            gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                       >  1:519629/1‑102 (MQ=255)                                                                            gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACAcc                                                                       <  2:519629/102‑1 (MQ=255)                                                                             ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  1:51610/1‑139 (MQ=255)                                                                             ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  2:27526/1‑139 (MQ=255)                                                                              gATAATTTTCAATAACTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCCATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCa                                <  1:243337/139‑1 (MQ=255)                                                                              gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCa                                <  1:929347/139‑1 (MQ=255)                                                                              gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCa                                >  2:929347/1‑139 (MQ=255)                                                                                tAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        >  2:482144/1‑129 (MQ=255)                                                                                tAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        <  1:482144/129‑1 (MQ=255)                                                                                 aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTa                                                                                                   >  2:99377/1‑69 (MQ=255)                                                                                 aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTa                                                                                                   <  1:99377/69‑1 (MQ=255)                                                                                 aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTc                             >  2:234529/1‑139 (MQ=255)                                                                                 aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTc                             <  1:234529/139‑1 (MQ=255)                                                                                        aaTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                            >  2:297350/1‑85 (MQ=255)                                                                                        aaTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTa                                                                            <  1:297350/85‑1 (MQ=255)                                                                                             aTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACtt                                                                                <  1:607697/76‑1 (MQ=255)                                                                                             aTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACtt                                                                                >  2:607697/1‑76 (MQ=255)                                                                                               tgtgCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGtct               <  1:99871/139‑1 (MQ=255)                                                                                                 tgCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCg             >  1:586273/1‑139 (MQ=255)                                                                                                  gCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGa            <  2:586273/139‑1 (MQ=255)                                                                                                        ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                <  1:363368/97‑1 (MQ=255)                                                                                                        ttCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                >  2:363368/1‑97 (MQ=255)                                                                                                           aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTAc                                                                                  >  2:18289/1‑60 (MQ=255)                                                                                                           aCCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTAc                                                                                  <  1:18289/60‑1 (MQ=255)                                                                                                            ccAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAAtata  <  1:27526/139‑1 (MQ=255)                                                                                                                gTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAaca                                                                         >  1:246261/1‑64 (MQ=255)                                                                                                                gTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAaca                                                                         <  2:246261/64‑1 (MQ=255)                                                                                                                   tgataaGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                          <  1:312444/108‑1 (MQ=255)                                                                                                                   tgataaGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                          >  2:312444/1‑108 (MQ=255)                                                                                                                       |                                                                                                                                GCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATA  >  NC_000913/3360651‑3360896

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