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%	182.7 / NA	56	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 (24/32);  total (24/32)

CAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATA  >  NC_000913/3360638‑3360896                                                                                                                                    |                                                                                                                                cAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAg                                                                                                                          <  2:131128/139‑1 (MQ=255) cAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAg                                                                                                                          <  1:409560/139‑1 (MQ=255)              gCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAAt                                                                                                             <  2:661756/139‑1 (MQ=255)                   tCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTtctc                                                                                                        >  1:752220/1‑139 (MQ=255)                       aTTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCAttt                                                                                                    >  2:703877/1‑139 (MQ=255)                              aCGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                             <  1:916781/139‑1 (MQ=255)                                      cgcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAAt                                                                                     <  1:311358/139‑1 (MQ=255)                                      cgcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAAt                                                                                     >  1:940029/1‑139 (MQ=255)                                      cgcgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAAt                                                                                     >  1:672134/1‑139 (MQ=255)                                        cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                                                   >  2:173099/1‑139 (MQ=255)                                        cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                                                   >  2:876769/1‑139 (MQ=255)                                        cgTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                                                   >  1:626543/1‑139 (MQ=255)                                           tgttTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACtt                                                                                <  1:217393/139‑1 (MQ=255)                                                     cAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACct                                                                      <  1:445728/139‑1 (MQ=255)                                                        aTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTg                                                                   >  2:285938/1‑139 (MQ=255)                                                                    cGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGAtt                                                       <  2:752220/139‑1 (MQ=255)                                                                           ccTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                <  2:60504/139‑1 (MQ=255)                                                                           ccTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                <  1:173099/139‑1 (MQ=255)                                                                           ccTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTCGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAg                                                <  1:706165/139‑1 (MQ=255)                                                                               tCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTAt                                            <  2:626543/139‑1 (MQ=255)                                                                                  tataATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCa                                         <  2:270383/139‑1 (MQ=255)                                                                                   ataataGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        <  2:672134/139‑1 (MQ=255)                                                                                   ataataGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        <  1:613270/139‑1 (MQ=255)                                                                                         gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        >  1:261194/1‑133 (MQ=255)                                                                                         gggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAg                                        <  2:261194/133‑1 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACtt                                                                                >  1:229396/1‑92 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACtt                                                                                <  2:229396/92‑1 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  2:662677/1‑139 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  2:152627/1‑139 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  1:587799/1‑139 (MQ=255)                                                                                          ggATAATTTTCAATAAATTGTGCTTTGTACACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTc                                 >  1:654620/1‑139 (MQ=255)                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                             >  1:195642/1‑78 (MQ=255)                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATcc                                                                                             <  2:195642/78‑1 (MQ=255)                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACATCTCTGAc                                                                 <  1:611853/106‑1 (MQ=255)                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCt                                                     <  1:714721/118‑1 (MQ=255)                                                                                           gATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCt                                                     >  2:714721/1‑118 (MQ=255)                                                                                           gATAATTTTCAATAAATTATGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACATCTGAc                                                                 >  2:611853/1‑106 (MQ=255)                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                           <  1:116437/95‑1 (MQ=255)                                                                                            aTAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTaa                                                                           >  2:116437/1‑95 (MQ=255)                                                                                             tAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCAtt                              <  1:876769/139‑1 (MQ=255)                                                                                              aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTc                             <  2:940029/139‑1 (MQ=255)                                                                                              aaTTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTc                             >  1:461652/1‑139 (MQ=255)                                                                                                 tttCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                          <  2:51548/139‑1 (MQ=255)                                                                                                 tttCAATAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCAt                          <  2:654620/139‑1 (MQ=255)                                                                                                     aaTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTAtt                                           <  2:115859/118‑1 (MQ=255)                                                                                                     aaTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTAtt                                           >  1:115859/1‑118 (MQ=255)                                                                                                      aTAAATTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTaaa                     <  1:703877/139‑1 (MQ=255)                                                                                                         aaTTGTGCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTg                  >  2:797496/1‑139 (MQ=255)                                                                                                               gCTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGa            <  1:662677/139‑1 (MQ=255)                                                                                                                cTTTGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGat           >  2:262145/1‑139 (MQ=255)                                                                                                                 tttGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGata          <  2:245763/139‑1 (MQ=255)                                                                                                                   tGTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATAtt        >  1:613604/1‑139 (MQ=255)                                                                                                                    gTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGtct               >  1:89873/1‑131 (MQ=255)                                                                                                                    gTTCACCAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGtct               <  2:89873/131‑1 (MQ=255)                                                                                                                         ccAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAAtata  <  2:587799/139‑1 (MQ=255)                                                                                                                         ccAAGTCTGATAAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAAtata  <  1:152627/139‑1 (MQ=255)                                                                                                                                    |                                                                                                                                CAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAAATTGTCTCGATATTAATATA  >  NC_000913/3360638‑3360896

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: ‑